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1.0  INTRODUCTION 


The  description  of  nonequilibrium  phenomena  in  gases  is  of  impor¬ 
tance  to  the  aerospace  sciences  because  these  phenomena  are  present 
in  such  diverse  applications  as  hypersonic  flow  fields,  rocket  exhaust 
plumes,  and  low  density,  high  enthalpy  wind  tunnels.  Experimental 
and  calculational  techniques  to  study  the  physicad  mechanisms  contri¬ 
buting  to  the  nonequilibrium  situation  are  at  best  difficialt  because  of 
the  complexity  of  the  necessary  mathematical  system  as  well  as  a 
lack  of  knowledge  of  many  of  the  fundamental  constants.  Although 
still  present,  these  drawbacks  are  reduced  in  decaying  low  density 
plasmas  and  such  plasmas  can  be  used  for  studying  the  fundamental 
manifestations  of  the  nonequilibrium  phenomena.  This  study  usually 
proceeds  through  interpretation  of  measured  spectral  line  intensities, 
electron  densities,  etc.  in  terms  of  recombination  rates  and  other 
useful  properties  of  the  plasma.  Proper  and  complete  interpretation 
of  these  data  requires  knowledge  of  the  distribution  of  population 
densities  among  the  different  energy  modes  of  the  plasma,  and  in 
the  nonequilibrium  environment,  this  is  generally  not  describable  by 
a  Boltzmann  distribution. 

The  majority  of  the  work  in  recent  years  has  been  concerned  with 
electron- ion  recombination  and  the  distribution  of  the  population  den¬ 
sities  of  electronic  states  in  atomic  plasmas.  In  most  of  this  work, 
the  assumption  that  the  population  densities  of  excited  states  are 
quiescent  was  required  in  order  to  make  the  calculations  tractable. 

This  assumption,  called  the  quasi-steady- state  (QSS)  assumption, 
although  apparently  reasonable,  results  in  the  masking  of  the  transient 
physical  mechanisms  that  contribute  to  establishing  the  final  distribu¬ 
tion.  Implicit  to  the  QSS  is  the  assumption  that  the  occurrence  of 
atomic  collisional  processes  is  sufficiently  rapid  that  macroscopic 
phenomena,  such  as  gas  dynamic  effects,  are  negligible. 

The  intent  of  this  study  is  to  remove  the  QSS  assumption  so  that  the 
detailed  processes  establishing  the  transient  behaviour  of  the  popula¬ 
tion  densities  can  be  studied  during  the  establishment  of  the  final 
distribution.  The  present  study  is  confined  to  helium  plasmas  because 
the  helium  atom  has  been  extensively  studied  both  experimentally  and 
theoretically  so  that  the  fundamental  quantum  mechanical  structure  is 
well  known.  The  study  further  ignores  the  possible  effects  introduced 
by  gas  dynamical  and  other  macroscopic  effects.  Including  these 
phenomena  at  this  point  would  only  complicate  the  problem  without 


9 


AEDC-TR-76-5 


materially  adding  to  the  fimdamental  description  of  the  necessary  mi- 
crosopic  processes. 


1.1  BACKGROUND 

The  following  review  of  the  wealth  of  literature  available  on  the 
subject  of  electron -ion  recombination  in  atomic  hydrogen  and  helium 
plasmas  is  not  an  exhaustive  one,  and  the  bibliography  is  not  complete 
by  any  means.  There  have  been  a  large  number  of  investigations  of 
recombination  in  helium,  and  enumeration  of  each  woiold  not  illuminate 
the  subject  materially  beyond  those  which  are  reviewed.  Rather,  the 
events  and  highlights  leading  to  the  present  study  are  indicated. 

Early  experiments  on  recombination  in  ionized  gases  provided 
recombination  rates  which  were  about  two  orders  of  magnitude  greater 
than  theoretical  values  obtained  by  considering  only  two-body  encoun¬ 
ters  (i.  e. ,  radiative  recombination)  (Ref.  1).  Early  in  1961,  D'Angelo 
(Ref.  1)  examined  the  recombination  rates  obtained  by  considering 
three -body  encounters  (two  electrons  and  an  ion),  resulting  in  a  neu¬ 
tral  atom  and  a  free  electron  for  hydrogen  plasmas.  The  captured 
electron  was  considered  to  be  in  a  bound  energy  state,  and  the  excited 
atom  subsequently  decayed  radiatively  to  the  ground  state.  The  results 
of  D'Angelo's  calculations  showed  that  considerations  of  three-body 
collisions  yield  recombination  rates  appreciably  larger  than  for  radia¬ 
tive  recombination.  Thus,  it  was  believed  that  three-body  collisions 
would  explain  the  large  recombination  rates  observed  in  plasmas.  This 
appears  to  have  been  the  first  work  in  which  the  three -body  collision 
theory  was  applied  to  recombination  studies. 

Byron,  Stabler,  and  Bortz  (Ref.  2)  subsequently  (1962)  further 
refined  the  calculative  procedure.  They  considered  three -body 
recombination  along  with  coUisional  de-excitation  and  radiative  de¬ 
excitation  and  calculated  recombination  rates  for  hydrogen  plasmas. 
Their  work  was  based  on  the  simplification  that  at  equilibrium  there 
is  a  minimum  in  the  total  rate  of  de -excitation  of  excited  atoms  with 
increasing  quantum  number.  This  minimum  results  because  the 
coUisional  de-excitation  rate  increases  with  increasing  quantum  num¬ 
ber,  caused  by  the  lower  energy  difference  between  states,  while  the 
radiative  de-excitation  rate  decreases  with  increasing  quantum  number 
due  to  the  smaller  transition  probabilities  for  states  of  larger  quan¬ 
tum  number.  They  then  solved  the  problem  for  the  total  recombi¬ 
nation  rate  by  obtaining  the  quantum  level  at  which  the  minimum 
occurs  and  finding  its  net  rate  of  de -excitation. 
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During  this  time.  Bates,  Kingston,  and  McWhirter  published  sev¬ 
eral  papers  concerning  recombination.  Their  work,  based  upon  a 
statistical  approach  to  the  problem,  culminated  in  two  comprehensive 
papers  in  1962  on  collisional -radiative  recombination  (Refs,  3  and  4). 
They  included  all  the  salient  processes  which  contribute  to  level  pop¬ 
ulations  of  a  hydrogenic  plasma  in  a  set  of  rate  equations.  The  result 
of  these  studies  was  a  model  which  could  be  generally  applied  to  decay¬ 
ing  plasmas;  the  model  is  usually  termed  the  "collisional-radiative 
recombination"  (CRR)  model.  In  this  model  a  Boltzmann  distribution 
of  population  of  upper  excited  states  is  maintained  which  is  in  equili¬ 
brium  with  the  free  electron  density. 

The  aforementioned  studies  form  the  basis  of  the  present  day  theo¬ 
ries  of  electron- ion  recombination  in  atomic  plasmas.  The  primary 
emphasis  in  each  has  been  to  obtain  recombination  rates  which  were 
more  consistent  with  experiment.  In  each  it  was  necessary  to  invoke 
the  QSS  assumption,  and  thus  information  as  to  the  mechanisms  for  the 
establishment  of  a  final  distribution  of  excited  state  population  densities 
was  lost. 

Before  the  establishment  of  the  CRR  model,  it  was  generally  nec¬ 
essary  to  revert  to  two  basic  models  for  the  analytic  study  of  excited 
state  population  densities.  One  of  these,  the  loc^  thermodynamic 
equilibrium  (LTE)  model  (Ref.  5)  is  useful  when  the  densities  are 
sufficiently  high  to  maintain  strong  collision  dominance  and  Boltzmann 
distributions  of  excited  state  densities  are  maintained.  The  other 
model,  the  corona  model,  applies  to  the  low  density  situation  (Ref.  5) 
where  radiative  rates  are  competitive  with  collisional  rates  and  the 
distributions  are  non-Boltzmann.  McWhirter  and  Hearn  (Ref.  6)  cal¬ 
culated  the  instantaneous  population  densities  for  the  excited  levels  of 
hydrogenic  plasmas  based  upon  the  CRR  model  for  a  variety  of  plasma 
conditions.  Their  results  are  in  basic  agreement  with  the  LTE  and 
corona  models  and  provided  validity  for  the  QSS  calculation  at  interme¬ 
diate  conditions  when  the  simpler  models  were  not  valid.  There  was  no 
way  in  using  these  models  for  excited  state  density  determinations  that 
the  validity  of  the  QSS  assumption  applied  to  a  real  plasma  could  be 
checked.  Rather,  all  that  was  known  was  that  the  various  theories, 
based  upon  the  QSS  assumption,  agreed. 

The  feeling  that  apparently  has  prevailed  since  about  1964  is  that 
the  recombination  rate  problem  has  been  solved  for  hydrogenic  plasmas. 
There  are  occasional  papers  in  which  the  CRR  ionization  and  recombi¬ 
nation  coefficients  are  computed  using  later  and  more  accurate  cross- 
section  calculations  (Ref.  7)  or  introducing  another  class  of  collisions 
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(Ref.  7),  or  including  estimated  relaxation  times  to  the  steady  state 
(Ref.  9).  This  work  all  depends  on  the  QSS,  however. 

The  difficulties  encountered  in  studies  of  recombination  in  helium 
plasmas  largely  paralleled  those  for  hydrogen  plasmas.  A  wealth  of 
experimental  data  from  microwave  and  spectroscopic  measurements  in 
decaying  afterglows  was  available,  which  indicated  significant  differences 
in  the  calculated  and  measured  recombination  rates.  Indeed  it  was  just 
such  discrepancies  that  prompted  the  analytic  work  in  the  hydrogenic 
plasmas.  Shortly  after  D'Angelo  (Ref.  1)  published  his  work,  Hinnov 
and  Hirschberg  (Ref.  10)  applied  the  three -body  theory  to  helium  recom¬ 
bination  experiments  in  the  B-1  Stellerator  at  Princeton.  By  also  includ¬ 
ing  inelastic  electron  collisions  for  excitation  and  superelastic  collisions 
between  bound  atomic  states  they  obtained  good  correlation  between 
theory  and  experiment.  In  later  work  Hinnov  and  Hirschberg  (Ref.  11) 
extended  the  interpretation  to  locating  a  ’’critical  level' '  which  is  defined 
as  the  lowest  principal  quantum  level  in  equilibrium  with  the  free  elec¬ 
trons.  This  critical  level  serves  as  an  upper  bound  on  the  quantum 
levels  affecting  the  electron  density  decay.  They  refer  to  favorable 
comparisons  of  their  work  to  the  preliminary  studies  which  led  to  the 
CRR  model  of  Bates,  Kingston,  and  Mc\^'hirter  (Refs.  3  and  4). 

Robben,  Kunkel,  and  Talbot  (Ref.  12)  made  a  spectroscopic  study 
of  the  freely  expanding  recombining  plume  of  a  helium  arc  jet  in  1963. 
They  determined  population  densities  of  upp>er  excited  states  including 
substates  for  each  principal  quantum  level  greater  than  two.  Their 
measurements  showed  that  the  sublevels  within  a  principal  quantum 
state  exhibited  Boltzmann  distributions  at  the  same  excitation  tempera¬ 
ture.  Recombination  rates  based  upon  the  QSS  were  approximately 
determined  and  were  in  reasonable  agreement  with  the  CRR  models  of 
Bates,  et  al.  (Refs.  3  and  4)  and  Hinnov  and  Hirschberg  (Ref.  11). 

Their  population  density  measurements  were  not  compared  to 
calculation. 

One  of  the  recurring  problems  of  experimental  studies  of  helium 
recombination  was  the  role  of  the  molecular  helium  ion  He"^.  In  1964 
Collins  and  Robertson  (Refs.  13,  14,  and  15),  with  a  series  of  selective 
excitation  experiments  on  a  flowing  afterglow,  and  Niles  and  Robertson 
(Refs.  16  and  17),  using  a  series  of  interference  filters  to  study  the 
si>ectral  emission,  foimd  that  He"^  recombination  dominated  the  after¬ 
glow  at  low  pressures  (less  than  about  20  torr)  and  above  this  the  mo¬ 
lecular  recombination  quickly  dominated  the  afterglow.  Both  species 
were  found  to  follow  CRR,  although  at  different  rates,  even  in  the 
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presence  of  the  other.  At  about  the  same  time,  Rogers  and  Biondi  (Ref. 
18)  published  the  results  of  experiments  on  the  helium  afterglow  in 
which  the  attempt  was  made  to  reconcile  the  observed  large  electron- 
ion  recombination  losses.  Although  the  goal  was  not  fully  realized,  they 
argue  that  dissociative  recombination  appeared  to  be  the  most  likely 
candidate  for  the  line  shapes  observed  in  the  late  afterglow. 

Collins  and  Hurt  (Ref.  19)  report  on  helium  afterglow  spectral  line 
measurements  made  in  1967  from  which  they  determined  population 
densities  of  eigenlevels  through  principal  quantum  number  25.  They 
found  that  there  were  three  separate  groups  of  quantum  levels  charac¬ 
terized  by  decay  times.  A  lower  group  and  an  upper  group  exist  whose 
quantum  levels  within  the  group  decay  with  the  same  lifetimes,  although 
the  lifetime  is  different  between  the  two  groups.  The  other  group  is  an 
intermediate  one  in  which  the  quantum  levels  within  the  group  all  decay 
with  different  lifetimes.  They  foxind  excitation  temperatures  of  300°K 
for  the  upper  group  and  2,  000*^  for  the  lower  group.  The  upper  group 
was  determined  to  be  in  equilibrium  with  the  free  electrons  by  compar¬ 
ing  exponential  decay  constants  of  the  levels,  the  free  electron  density, 
and  a  calculated  decay  constant  for  He”^.  By  a  similar  comparison  it 
was  .determined  that  the  decay  rate  of  the  lower  group  is  proportional 
to  the  square  of  the  electron  density  and  the  first  power  of  the  concen¬ 
tration  of  He’*'.  In  a  follow-up  analysis,  Collins  (Ref.  20)  computed  the 
logarithmic  derivative  of  excited  state  populations  as  functions  of 
logarithmic  electron  density.  These  calculations,  based  upon  the  quasi¬ 
steady-state  approximation,  yielded  recombination  rates  in  semiquanti- 
tative  agreement  with  experiment  (Ref.  19). 

Recent  studies  of  helium  recombination  of  some  significance  are 
those  of  Chen  (Ref.  21)  and  Johnson  and  Hinnov  (Ref.  22),  both  published 
in  1969.  Chen  appears  to  have  made  the  first  direct  calculation  of  helium 
recombination  rates,  as  well  as  those  of  other  atoms,  based  on  the  QSS 
assumption.  Chen's  approach  was  the  same  as  that  of  Byron,  et  al. 

(Ref.  2)  and  Bates,  et  al,  (Ref.  23)  for  hydrogen,  in  which  the  minimum 
in  the  total  rate  of  de-excitation  of  atoms  as  a  function  of  the  energy 
level  of  the  excited  states  is  found.  Chen  also  made  experimental  mea¬ 
surements  of  recombination  rates  in  a  discharge  tube  and  found  more  or 
less  satisfactory  agreement  with  calculation. 

Johnson  and  Hinnov  (Ref.  22)  made  spectroscopic  measurements  of 
excited  state  population  densities  in  the  C-Stellerator  at  Princeton  and 
compared  them  to  quasi-steady-state  calculations  of  the  densities  using 
various  assumed  formtilations  for  collisional  cross  sections  and  at 
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various  excitation  temperatures.  These  measurements  were  at  quite 
low  pressures  so  that  the  presence  of  molecular  helium  could  be  discount¬ 
ed.  The  comparison  between  their  calculated  and  measured  recombina¬ 
tion  rates  was  quite  good  for  moderate  times  in  the  afterglow.  At  early 
times,  however,  the  calculated  rates  were  noticeably  smaller  than  the 
measured  rates. 

The  aforementioned  studies  of  helium  recombination  have  assumed 
the  QSS  a  physical  as  well  as  a  mathematical  reality,  just  as  in  the 
analytic  studies  of  hydrogenic  recombination.  Little  of  that  previously 
cited  work  was  concerned  with  comparison  of  measured  and  calculated 
excited  state  densities  in  the  region  of  nonequilibrium.  The  work  of 
Johnson  and  Hinnov  (Ref.  22)  is  an  exception  to  this  but  they  also  use 
the  QSS  and  this  leads  to  discrepancies.  This  is  discussed  further  in 
Section  7.  0. 

The  solution  to  the  original  system  of  rate  equations  from  which 
Bates,  et  al.  (Refs.  3  and  4)  obtained  the  CRR  model  does  not  demand 
the  QSS  assumption,  although  the  calculations  are  made  much  easier  by 
its  use.  A  solution  of  the  complete  set  of  Eigenstate  Rate  Equations 
(ERE)  is  possible  using  modern  computational  methods,  and  makes 
possible  a  more  detailed  examination  of  the  approach  to  the  QSS  and 
permits  estimates  to  be  made  of  the  time  to  reach  the  QSS  solution  to 
the  ERE. 

Since  1967,  two  groups  of  investigators  (Refs.  24,  25,  26,  and  27) 
have  attacked  the  problem  of  removing  the  QSS  assumption  for  hydrogen 
plasmas.  Limbaugh,  Carstens,  McGregor,  and  Mason  (Refs.  24,  26, 
and  27)  were  primarily  concerned  with  the  pure  afterglow  problem  in 
which  the  electron  temperature  remains  constant.  The  purpose  of  that 
investigation  was  to  examine  the  fundamental  consistency  of  the  collision- 
al- radiative  recombination  model  to  the  full  transient  solution  and  to 
determine  quantitatively  the  differences  in  the  decay  of  plasmas  in  ^ 
which  a  radiative  metastability  existed.  Among  the  findings  was  that 
the  metastable  state  definitely  affected  the  decay  and  that  under  reason¬ 
able  conditions,  there  could  be  a  physically  significant  time  before  the 
QSS  is  established.  Gordiets,  Gudsenko,  and  Shelepin  (Ref.  25)  were 
primarily  concerned  with  the  effects  of  sudden  changes  in  the  plasma 
parameters,  such  as  excitation  temperature  or  free  electron  density 
upon  the  plasma  decay.  They  show  that  in  such  cases  the  parametric 
changes  could  have  an  appreciable  effect  upon  the  plasma  decay.  Al¬ 
though  detailed  quantitative  comparisons  between  the  results  of  the  two 
studies  are  difficult  because  of  the  different  problems  being  solved,  they 
do  show  qualitative  agreement  with  each  other. 
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1.2  APPROACH  AND  GOALS 

The  relative  success  of  those  earlier  studies  illustrated  the  potential 
of  the  ERE  to  provide  a  detailed  accounting  of  the  interaction  between 
free  electrons  and  the  various  energy  states  of  the  afterglow  plasma. 

This  success,  along  with  the  realization  that  only  through  use  of  the  ERE 
rather  than  CRR  theory  could  the  acquisition  of  the  QSS  be  studied,  and  a 
curiosity  about  the  processes  which  bring  a  plasma  to  the  QSS,  provided 
the  underlying  motivation  for  the  study  reported  here. 

For  these  purposes,  it  was  decided  to  extend  the  earlier  work, 
summarized  in  Ref.  26,  to  the  helium  atom.  Helium  is  next  to  hy¬ 
drogen  in  simplicity  and  has  been  investigated  extensively  experimental¬ 
ly.  In  addition,  its  electronic  structure  seems  to  promise  more  interest¬ 
ing  study  since  the  low-lying  energy  levels  are  farther  apart  thzin  those 
of  hydrogen  and  contain  both  singlet  and  triplet  states.  The  larger  energy 
spacing  causes  the  collisional  effects  to  be  less  effective  than  for  hydro¬ 
gen.  To  the  author's  knowledge,  this  is  the  first  application  of  the  ERE 
to  helium  and  the  present  computer  program  is  the  only  one  in  existence 
which  can  obtain  the  transient  solution  to  the  ERE. 

The  purpose  of  this  study  is  to  investigate  the  basic  mechanisms  by 
which  the  QSS  is  established  and  to  attempt  to  assess  the  validity  of  CRR 
theory  as  applied  to  physically  realizable  plasmas.  To  achieve  this,  the 
transient  solution  to  the  ERE  with  the  restrictions  and  assumptions  to  be 
discussed  below  will  be  obtained  for  several  plasma  conditions.  These 
studies  break  into  two  groups:  <1)  pure  afterglows  and  (2)  perturbations. 
The  pure  afterglows  are  studies  in  which  the  plasma  is  allowed  to  relax 
from  an  initial  LiTE  distribution  toward  the  QSS  distribution.  The  pertur¬ 
bations  are  studies  in  which  certain  plasma  parameters  are  perturbed 
from  the  acquired  QSS  and  the  relaxation  of  the  plasma  is  observed. 

In  the  realization  of  the  final  goal  of  this  study  there  are  several 
significant  objectives: 

1.  To  determine  the  mechanisms  by  which  the  QSS  is 
established  for  helium  plasmas. 

2.  To  observe  the  decay  characteristics  of  both  the 
excited  state  distribution  and  the  free  electron 
density  in  He  plasmas  and  to  determine  the 
applicability  of  the  results  of  CRR  theory  to 
helium  in  the  early  portion  of  plasma  decay. 
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3.  To  observe  the  variation  in  the  time  for  QSS  to  be 
established  in  He  plasmas  for  different  values  of 
electron  temperature,  electron  density,  and 

total  number  density,  and  to  observe  the  systematic 
trends. 

4.  To  gain  insight  into  the  total  coupled  processes 
which  plasmas  decay. 

In  the  following  pages.  Section  2.  0  contains  the  theoretical  model  and 
indicates  those  physical  processes  which  are  ignored  and  the  physical 
assumptions  which  are  made  in  setting  up  the  theoretical  description  of 
the  plasma.  Section  3.  0  includes  a  brief  examination  of  the  various  rate 
coefficients  and  transition  probabilities  used  in  the  ERE  and  indicates  the 
accuracy  of  the  specific  technique  used  to  evaluate  the  coUisional  cross 
sections.  In  Section  4.  0  the  numerical  techniques  used  to  effect  the 
solution  to  the  ERE  are  described  and  the  limitations  of  the  present 
computer  program  are  indicated.  In  Section  5.  0  and  6.  0  the  results  of 
the  computations  are  reported,  and  in  Section  7.  0  the  data  are  inter¬ 
preted  in  terms  of  effects  on  measurements  and  application  to  decaying 
helium  afterglow  plasmas.  The  findings  of  this  study  and  comments 
upon  their  impact  on  the  interpretation  of  existing  measurements  and 
the  design  of  future  experiments  are  summarized  in  Section  8.  0.  Also 
included  in  Section  8.  0  is  a  discussion  of  recommended  future  research 
using  the  ERE,  and  certain  simplifications  which  should  be  incorpor¬ 
ated  in  further  work  are  pointed  out. 


2.0  TRANSIENT  PLASMA  BEHAVIOR 


This  section  examines  the  development  of  the  system  of  equations 
which  will  describe  the  transient  properties  of  the  popvilation  density  of 
the  various  internal  energy  states  available  to  a  singly  ionized  mona¬ 
tomic  plasma.  In  the  following  subsections  a  generalized  description  of 
the  physical  problem  shows  the  role  of  the  separate  processes  contri¬ 
buting  to  the  excitation  or  de- excitation.  Subsequent  to  this  discussion 
each  of  the  phenomena  will  be  discussed  separately  leading  to  a  detailed 
mathematical  description  of  each  process.  Included  in  this  discussion 
will  be  the  simplifications  and  approximations  which  are  necessary  to 
make  the  mathematical  system  tenable  to  the  physical  problem  at  hand. 
Finally  the  detailed  descriptions  are  combined  to  illustrate  the  full 
description  of  the  physical  problem  and  the  development  of  the  CRR 
model  is  described. 
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2.1  THE  GENERALIZED  PROBLEM 

To  provide  a  detailed  description  of  a  plasma  which  may  not  be  in 
equilibrium  it  is  necessary  to  include  each  energy  state  available  to 
each  plasma  constituent.  By  considering  some  arbitrary  control  volume 
in  space  and  examining  each  energy  state  within  the  control  volume  it 
is  possible  to  arrive  at  a  simple  model  for  the  plasma.  The  dimensions 
of  the  control  volume  will  be  assumed  to  be  small  with  respect  to  the 
total  extent  of  the  plasma  yet  large  enough  that  there  are  a  sufficient 
number  of  particles  contained  within  it  for  statistical  properties  to  be 
meaningful.  Any  property  measured  or  ascribed  to  the  constituents  of 
the  plasma  within  the  control  volume  is  assumed  constant  throughout 
the  control  volume.  If  the  plasma  exists  in  field-free  space,  then  the 
population  density  of  a  given  energy  state  will  be  affected  by  local 
interactions  between  the  constituents  of  the  plasma  contained  within  the 
control  volume  and  by  a  net  flow  of  the  plasma  constituents  across  the 
boundaries  of  the  control  volume. 

The  time  rate  of  change  of  the  population  density  of  some  specific 
energy  level  within  the  control  volume  can  be  written  symbolically: 


dn(p) 

dt 


micro 


space 


(1) 


where  n  is  population  density,  p  symbolically  represents  the  set  of 
quantum  numbers  necessary  to  describe  the  particiilar  energy  state, 
t  is  time,  and  the  subscripts  "micro"  and  "space"  refer  to  the  micro¬ 
scopic  processes  and  spatial  variation  described  above.  An  equation  of 
this  general  type  must  be  written  for  each  of  the  energy  states  available 
to  each  plasma  constituent  and  the  solution  of  the  resultant  system  of 
equations  will  yield  a  time -dependent  function  for  the  population 
density  of  each  state. 

Including  that  portion  of  Equation  (1)  which  refers  to  the  spatial 
variation  of  n  (p)  will  generally  necessitate  including  the  statements  of 
conservation  of  energy  and  momentum  as  well  as'  the  statement  of  con¬ 
servation  of  mass.  As  discussed  earlier,  the  effort  here  is  an  examina¬ 
tion  of  the  mechanisms  by  which  the  QSS  is  obtained  and  an  assessment 
of  the  importance  of  including  the  transient  coupling  in  analysis  of  helium 
plasmas.  The  present  work  is  to  be  considered  a  first  step  leading  to 
eventual  modeling  of  plasma  phenomena,  and  inclusion  of  the  spatial 
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variations,  although  necessary  for  final  realization  of  the  goal,  is 
inappropriate  at  this  stage  of  development.  Hence,  subsequent  dis¬ 
cussion  will  consider  all  of  space  as  being  filled  with  the  plasma,  with 
no  sources  or  sinks,  and  with  no  gradients  present. 

2.2  THE  EIGENSTATE  RATE  EQUATION  (ERE) 

The  physical  phenomena  comprising  the  contributions  to  the  tran¬ 
sient  description  of  the  population  density  on  the  microscopic  level  can 
be  treated  by  examining  the  plasma  within  the  control  volume  as  if 
there  were  no  convective  variations.  Within  this  control  volume,  the 
monatomic  plasma  consists  of  neutral  atoms  in  each  of  the  possible 
electronic  energy  states,  ions,  and  free  electrons.  For  this  work  which 
will  be  concerned  with  the  helium  atom  single  ionization  will  be  assumed 
so  that  the  concentration  of  ions  and  free  electrons  will  be  identical. 
Further,  the  electrons  will  be  assumed  to  exhibit  a  Maxwell-Boltzmann 
velocity  distribution  independent  of  the  kinetic  temperature  of  the  atoms 
or  other  heavy  bodies. 

Each  of  the  mechanisms  causing  a  change  in  the  population  density 
of  a  particuleir  eigenstate  can  be  characterized  by  a  rate  coefficient. 

The  instantaneous  rate  by  which  the  level  is  being  populated  or  depopu¬ 
lated  by  a  given  mechanism  is  the  product  of  the  rate  coefficient  and 
the  concentrations  of  the  reaction  partners.  The  microscopic  processes, 
the  instantaneous  rates,  and  the  rate  coefficients  included  in  this 
study  can  be  listed  as  follows: 

1.  Inelastic  and  super -elastic  electron- atom  collisions. 

For  inelastic  excitation,  the  necessary  energy  to 
effect  the  transition  is  supplied  by  the  kinetic  energy 
of  the  free  electron.  For  the  super-elastic  de-exci¬ 
tation,  the  excess  energy  increases  the  kinetic  energy 
of  the  free  electron.  The  rate  coefficient  for  this  reac¬ 
tion  is  K<p;q)  where  p  and  q  represent  initial  and  final 
energy  states  of  the  bound  electron.  The  instaintaneous 
rate  for  this  process  is  given  by 

^§4^  =  -n^  n(p)  K(p;q)  (2) 

for  depopulating  of  state  p  to  state  q.  In  the  above,  n 
represents  population  densities  and  the  subscript  e 
refers  to  the  free  electron  density.  Obviously  the  same 
expression  with  a  plus  sign  will  describe  populat¬ 
ing  of  State  q  from  state  p.  Units  of  K(p;q)  are  cm^/sec. 
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2.  Electron-ion  encounters  in  which  a  single  ion  captures 
a  free  electron  into  some  bound  energy  state  with  an 
attendant  emission  of  radiation.  This  is  more  commonly 
called  radiative  recombination  and  is  described  by 

the  radiative  recombination  coefficient,  j3  {p).  The 
instantaneous  rate  of  filling  of  state  p  is  given  by 

dn  (p)  +  „  ,  . 

where  n"*"  is  the  ion  concentration.  Units  of  ^  (p)  are 
cm3 /sec. 

3.  Electron- electron- ion  encounters  resulting  in  a 
neutral  atom  with  one  electron  in  a  boimd  energy  state 
and  the  excess  energy  contributing  to  the  kinetic  energy 
of  the  free  electron.  This  three -body  recombination 
reaction  is  characterized  by  K  (c;p)  where  c  symbolizes 
the  electron  continuum  as  the  initial  state  of  the  reaction 
partners.  The  instantaneous  rate  is  given  by 

n+  K(c;p)  (4) 

Units  of  K  (c;p)  are  cm^/sec. 

4.  Ionization  collisions  in  which  a  free  electron  impinges  upon 
a  neutral  atom  with  sufficient  energy  to  ionize  the  atom 
resulting  in  a  singly  charged  ion  and  two  free  electrons. 

This  reaction  might  be  thought  of  as  the  equilibrium 
partner  of  the  three -body  recombination  described  above  and 
is  characterized  by  K  (p;c).  The  instantaneous  rate  is 

"  -He  n(p)  K(p;c)  (5) 

Units  of  K  (pjc)  are  cm^/sec. 

5.  Spontaneous  radiative  transitions  in  which  a  bound 
electron  relaxes  spontaneously  to  a  lower  bound  state 
with  an  attendant  emission  of  radiation.  This  mechanism 
can  contribute  to  both  populating  and  depopulating  a  level 
and  is  characterized  by  the  Einstein  transition  probability 
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A  (p;q).  The  instantaneous  rate  for  depopulating  state  p 
to  some  lower  state  q  is 

"  -n(p)  A(p;q)  (6) 

while  the  instantaneous  rate  for  populating  state  p  from 
some  higher  state  m  is 

=  n(m)  ‘A(m;p) 

Units  of  A  (p;q)  are  1/sec. 

In  addition  to  the  above  processes  other  phenomena  which  must  be 
considered  for  their  effect  upon  the  population  densities  are:  atom-atom 
and  atom- ion  collisions,  photon  absorption,  induced  emission,  and  the 
molecular  ion  effects.  This  study  will  involve  pressures  sufficiently 
low  that  the  molecular  ion  can  be  ignored  (Ref.  14). 

Heavy  body  collisions  can  be  ignored  for  this  application  also 
because  of  the  low  densities  and  heavy  body  temperature  examined  in  this 
study.  For  example,  for  a  total  particle  density  of  10^^  l/cm^  and  a 
heavy  body  temperature  of  3,  000“K,  an  elevated  state  population  density 
typical  for  an  electron  velocity  distribution  at  10,  000°K  is  of  order  10^ 

1/  cm^.  Considering  collisions  between  excited  state  atoms  and  assuming 
an  atomic  radius  of  1  A,  one  finds  the  collision  frequency  to  be  of  order  102 
encounters /sec/ cm3  foj.  these  atom-atom  collisions.  The  electron-atom 
collision  frequency,  assuming  an  electron  temperature  of  10,000°K  and 
1 -percent  ionization,  is  of  order  10 encounters /sec/ cm^.  Obviously 
the  atom- atom  encounters  are  insignificant  in  this  study  compared  to 
electron- atom  collisions. 

Absorption  and  stimulated  emission  can  also  be  discounted  for  the 
purposes  of  this  study.  Since  the  stimulated  emission  and  absorption 
coefficients  differ  only  by  the  ratio  of  statistical  weights,  only  absorp¬ 
tion  of  resonance  radiation  can  be  significant.  Accurate  accounting  of 
absorption  requires  including  geometry  considerations,  and  this  is 
beyond  the  scope  of  this  study.  An  estimate  of  absorption  effects  is 
obtained  from  Robben,  Kunkel,  and  Talbot  (Ref.  12)  who  found  that  for 
a  heavy  body  temperature  of  1,000‘’K  and  particle  density  of  order  10 
l/cm3  the  mean  free  path  of  the  helium  l^S  -  2^?  resonance  radiation 
was  about  0.  2  mm.  This  distance,  although  small  enough  to  make  it 
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necessary  to  consider  resonance  radiation  absorption  in  many  physical 
experiments,  is  svifficiently  large  that  there  will  not  be  appreciable  ab¬ 
sorption  in  the  small  control  volume  on  which  the  present  study  is  based. 
Further,  as  will  be  seen,  the  particle  density  used  in  the  estimate  (Ref. 
12)  is  near  the  upper  bound  of  the  total  densities  of  the  plasmas  studied 
here.  Thus,  absorption  can  be  neglected  for  the  purposes  of  this  study 
and  the  plasmas  are  considered  to  be  optically  thin.  For  one  set  of 
plasma  conditions,  comparative  calculations  are  made  at  different 
optical  thicknesses  to  examine  the  effect  upon  the  mechanisms  which 
QSS  is  acquired. 


The  total  time  rate  of  change  of  the  population  density  of  state  p  is 
just  the  algebraic  sum  of  the  five  basic  processes  described  above: 


3n  (p) 
St 


-n  n(p)[  Z  K(p;q) 
q=l 


+  K(p;c)  ] -n  (p)  Z  A(p;q) 
q<p 


+n  Z  n(q)  K  (q;p) 

q=i 


+  Z  n(q)  A(q;p) 
q>p 


n"*"  K(c;p)  +  n^  n'*'  3  (p) 


(7) 


In  the  formulation  of  Equation  (7)  all  K(p;p)  are  zero  and  any  A(p;q) 
is  zero  if  the  dipole  selection  rules  are  violated. 

In  Eq.  (7)  the  first  three  terms  express  the  rate  with  which  state  p 
is  depopulated  due  to  collisional  internal  transitions  (Eq.  (2)),  collision- 
al  ionization  (Eq.  (5)),  and  spontaneous  radiative  transitions  (Eq.  (6)), 
respectively.  The  remaining  terms  express  populating  rates  due  to  col¬ 
lisional  internal  transitions  (Eq.  (2)),  radiative  transitions  (Eq.  (6)), 
three-body  recombination  (Eq.  (4)),  and  radiative  recombination 
(Eq.  (3)). 

There  will  be  a  similar  equation  for  each  eigenstate  avilable  to  the 
atom.  In  addition,  the  continuum  density  n0  enters  into  each  equation 
as  a  product  with  each  of  the  other  bound  state  densities.  Hence,  the 
resultant  system  of  equations  is  a  rectangular  non-linear  system  of  in¬ 
finite  extent.  However,  by  requiring  that  there  by  no  loss  or  gain  of 
atoms  hy  diffusion  or  convection,  the  totsd  concentration  will  be  constant. 
Thus,  taking  advantage  of  single  ionization  so  that  the  electron  and  ion 
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densities  are  identical,  the  conservation  of  heavy  particles  is  expressed 
by  the  equation 

00 

In  (p)  +  n  =  n  (8) 

p=l  ®  ° 

where  no  is  the  total  heavy  particle  concentration.  Equation  (8)  added 
to  the  system  causes  the  system  to  become  square  and  in  princijde  de¬ 
terminate. 

The  system  of  equations  has  the  capability  of  yielding  a  complete 
description  of  the  instantaneous  population  density  distribution  for  those 
monatomic  plasmas  in  which  there  are  electronic  transitions  only  be¬ 
cause  of  reactions  between  the  microscopic  constituent  particles. 
Additional  processes,  such  as  those  ignored  in  this  study,  can  be  account¬ 
ed  for  by  expressing  them  functionally  the  same  way  as  processes  1 
through  5. 

Although  the  system  is  infinite  in  extent,  thus  making  the  solution 
impossible,  it  may  be  truncated  at  some  finite  energy  level  as  a  reason¬ 
able  approximation.  This  is  because  for  progressively  higher  energy 
levels  the  population  densities  decrease  monotonically  with  energy  for 
an  equilibrium  plasma.  Hence,  there  will  be  an  upper  energy  level 
above  which  all  levels  have  population  densities  insignificant  with  respect 
to  the  lower  energy  levels  and  the  continuum.  Also,  reactions  between 
plasma  constituents  causing  a  transfer  to  bound  electrons  between  the 
low-lying  states  and  upper  states  will  be  insignificant  compared  to 
transfer  rates  between  energy  states  lying  close  together.  Finally, 
because  of  the  close-lying  energy  levels  of  the  upper  states  and  their 
close  proximity  to  the  ionization  potential,  reactions  with  the  free 
electrons  will  estdblish  an-equilibrium  configuration  among  these  upper 
states  as  well  as  with  the  continuum.  Hence,  because  of  the  relative 
insignificance  of  the  densities  of  the  high-lying  levels,  the  relative  un¬ 
importance  of  reaction  rates  between  the  low -lying  and  high-lying  energy 
levels,  and  because  there  will  be  some  upper  energy  level  which  is 
assuredly  in  equilibrium  with  the  continuum,  the  system  of  equations 
can  be  truncated  at  this  upper  level  with  minimal  loss  of  accuracy. 

The  minimum  energy  level  for  which  this  is  true,  of  course,  depends 
upon  the  kinetic  temperature  exhibited  by  the  free  electrons.  It  should 
also  be  noted  that  this  is  in  effect  postvilating  the  existence  of  the  so- 
called  critical  level,  which  has  been  described  as  the  lowest-lying 
energy  level  in  equilibrium  with  the  continuum.  For  the  purposes  of 
this  study,  the  upper  level  considered  was  dictated  more  by  available 
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computer  size  than  by  physical  considerations  and  thus  the  critical 
level  will  be  indicative  of  limits  on  plasma  conditions  which  can  be 
studied  with  the  present  program. 

2.3  THE  COLLISIONAL-RADIATIVE  RECOMBINATION  (CRR)  MODEL 

The  system  of  Eqs,  (7)  and  (8),  with  the  truncation  considerations 
discussed  above,  provides  the  starting  point  for  the  development  of  the 
CRR  model  {Refs.  3  and  4)  which  previous  investigators  have  used  to 
study  plasma  decays.  This  model  for  optically  thin  plasmas  utilizes 
the  QSS  approximation  which  assumes  that  the  time  rate  of  change  is 
zero  for  all  bound  states  other  than  the  ground  state.  The  CRR  model 
is  applied  to  those  plasmas  in  which  the  ground  state  density  and  free 
electron  density  are  much  larger  than  excited  state  densities  and  in 
which  the  electronic  kinetic  temperature  is  much  less  than  the  first 
excitation  energy  so  that  collisional  transfers  from  the  ground  state  are 
small.  These  conditions  lead  to  the  assumption  that  a  quasi-equilibrium 
distribution  is  established  "almost  instantaneously"  without  appreciable 
effect  on  the  continuum  density  and  that  collisional  excitation  and  de¬ 
excitation  rates  will  so  largely  dominate  the  decay  rate  for  the  plasma 
that  the  QSS  configuration  is  maintained.  Consequently,  the  deriva¬ 
tives  for  the  time  rate  of  change  of  these  excited  levels  can  be  set  to 
zero.  In  such  a  situation  the  time  rate  of  change  of  the  ground  state  is 
the  algebraic  negative  of  the  time  rate  of  change  of  the  free  electrons, 
li  (1)  =  -lig.  The  excited  state  densities  are  constant  and  the  system, 
Eqs.  (7)  and  (8),  becomes  a  single  differential  equation  and  a  system 
of  algebraic  equations.  The  only  variables  are  the  ground  state 
density  and  the  continuum  density  so  that  Eq,  (7)  for  the  ground  state 
can  be  written 

=  an  ^  -  S  n(l)  n  (9) 

dt  e  e 

where  a  and  S  are  defined  as  the  collisonal-radiative  recombination  and 
collisional-radiative  ionization  coefficients,  respectively.  The  a  and  S 
can  be  expressed  as  functions  of  ng,  electron  temperature,  and  atomic 
parameters  by  reduction  of  the  algebraic  portion  of  the  system.  The 
determination  of  a  and  S  under  these  approximations  was  the  problem  to 
which  Bates,  Kingston,  and  McWhirter  (Refs.  3  and  4)  addressed  them¬ 
selves. 

Although  the  values  of  a  and  S  for  various  plasmas  are  periodically 
improved  and  updated  as  knowledge  progresses,  the  basic  model  utiliz¬ 
ing  the  QSS  assumption  and  depicted  by  Eq.  (9),  or  some  variations  of 
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Eq.  (9),  has  become  the  standard  model  by  which  plasma  decay  and 
growth  phenomena  are  analyzed.  Even  in  applications  where  the  physi¬ 
cal  problem  supports  transient  density  or  electron  temperature  changes 
the  QSS  condition  is  assumed  to  establish  itself  in  infinitesimal  times 
compared  to  other  characteristic  times  for  the  plasma. 

In  most  applications,  as  in  the  case  here,  the  collisional -radiative 
ionization  coefficient,  S,  is  an  ignorable  quantity.  In  this  case  the 
recombination  of  the  free  electrons  depends  only  upon  the  CRR  coeffi¬ 
cient,  a. 


The  very  assumption  of  QSS  is  alien  to  a  detailed  study  of  the 
mechanisms  by  which  a  plasma  decays.  The  only  way  to  examine  this 
is  to  actually  solve  the  transient  problem  with  all  the  non-linearities 
and  coupling  between  states  for  several  plasma  conditions  and  this  is 
the  subject  of  this  study. 


2.4  THE  CRITICAL  LEVEL  AND  EQUILIBRIUM  CONFIGURATION 


The  concept  of  the  so-called  critical  level  is  very  important  to 
simplified  methods  of  obtaining  recombination  rates  (Refs.  11,  21,  and 
22).  The  concept  is  used  here  only  to  provide  a  bound  on  the  plasma 
conditions  for  which  the  present  computer  program  is  valid. 


Physically  the  critical  level  is  that  lowest  energy  level  which  is 
maintained  in  equilibrium  with  the  free  electron  density.  This  equili¬ 
brium  density  is  given  by  the  modified  Saha  equation: 


g  u- 

n(p)  =  n.  n  = — r-sr 

^  1  e  2g.  27rinkT 

^1  e 


3/2 


exp  (I  /kT  ) 
P  ® 


where  n^  and  n^  are  the  groimd  state  ion  and  electron  densities,  respec¬ 
tively,  and  Ip  is  the  ionization  potential  of  quantum  level  p.  All  quantum 
states  above  the  critical  level  will  also  be  in  equilibrium  with  the  free 
electrons.  Thus  the  net  rate  with  which  botind  electrons  above  and 
including  the  critical  level  fall  to  energy  levels  below  the  critical  level 
is  immediately  taken  up  by  recombination  of  the  free  electrons  into 
these  upper  levels  and  this  is  the  recombination  rate. 
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The  critical  state  in  other  studies  (Ref.  2)  is  found  by  formally 
minimizing  the  recombination  coefficient  with  respect  to  energy  state. 
Because  of  the  form  of  the  collisional  rate  coefficients  Section  3.0, 
evaluation  can  become  quite  tedious. 

Rather  than  attempt  a  formal  analysis  for  the  critical  level  in  this 
study,  the  approach  will  be  to  observe  the  poptilation  densities  of  the 
various  quantum  levels  in  relation  to  the  value  which  would  be  in  theo¬ 
retical  equilibrium  with  the  free  electron  density  calculated  from  the 
modified  Saha  equation.  In  this  manner  the  critical  level  can  be  empi¬ 
rically  observed  independently  of  problems  attendant  with  evaluation 
of  complicated  expressions. 

A  term  which  is  used  in  this  study  which  merits  explanation  is 
"equilibrium  configuration.  "  This  term  applies  to  any  two  or  more 
excited  states  whose  densities  satisfy  the  Boltzmann  relationship 

^  =  i  exp  t-(Ep-Eg)AT^l 

The  densities  of  states  which  are  in  an  equilibrium  configuration 
satisfy  the  above  equation  regardless  of  their  relationship  to  the  rest 
of  the  distribution. 

s 

3.0  TRANSITION  PROBABILITIES,  RATE  COEFFICIENTS, 

AND  CROSS  SECTIONS 


Although  the  detailed  rate  coefficients  and  transition  probabilities 
are  not  in  themselves  the  object  of  this  study,  they  are  critical  compo¬ 
nents  of  the  full  ERE  and  CRR  models.  The  magnitude  of  each  of  the  rate 
coefficients  and  transition  probabilities  indicates  to  what  extent  a  parti¬ 
cular  type  of  transition  will  occur  and  the  relative  importance  of  the 
various  mechanisms.  This  dependence  of  the  results  of  either  of  the  ERE 
or  CRR  models  upon  the  various  coefficients  and  probabilities  makes  it 
desirable  to  look  at  the  detailed  evaluation  of  these  parameters. 


3.1  EINSTEIN  TRANSITION  PROBABILITY 

The  Einstein  spontaneous  transition  probability  is  the  rate  with 
which  atoms  radiatively  decay  from  an  excited  electronic  state  to  a 
lower  electronic  state.  The  value  can  be  obtained  from  time -dependent 
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perturbation  theory  if  one  has  at  hand  the  wave  functions  of  the  various 
bound  states.  The  transition  probabilities  presently  available  have 
been  compiled  in  several  publications  (e.  g. ,  Refs.  28  and  29)  but  neither 
these  nor  any  other  reported  results  found  in  the  literature  include  much 
detail  concerning  transitions  involving  high  principal  quantum  numbers 
(P  ^  5).  However,  for  these  higher  energy  states  the  wave  functions 
will  be  sufficiently  similar  to  hydrogen  wave  functions  that  transition 
probabilities  for  hydrogen  and  helium  will  be  approximately  equal.  Care 
must  be  exercised  as  this  approximation  will  only  be  valid  for  interac¬ 
tions  among  upper  states  with  some  unspecified  lower  bound. 

Niles  (Ref.  28)  has  tabulated  most  of  the  radiative  transition  proba¬ 
bilities  available  at  that  time  (1967).  The  compilation  includes  data 
for  all  transitions  from  levels  with  principal  quantum  number  p  <  15 
to  states  with  principal  quantum  number  p  <  3.  In  addition  to  these,  the 
compilation  includes  transitions  terminating  on  quantum  level  4  except 
for  the  4^'  ^  F  -  p^>  ^  G  transitions.  The  transitions  listed  include  the 
interactions  of  those  lower  eigenstates  where  strong  departure  from  the 
hydrogenic  approximations  should  be  observed.  Hence,  for  these  transi¬ 
tions  to  the  lower  energy  states,  the  values  reported  by  Niles  (Ref.  28) 
were  used. 

Transition  probabilities  for  the  4^»^  F  -  p^'^  G  and  all  higher 
transitions  were  computed  using  the  hydrogenic  approximation.  These 
values  for  the  transition  from  state  (p,  f)  to  state  (p'  £.')  are  given  by 
(Ref.  30) 


4 

A(p,il;p  Z  )  =  - 

3h  c 

where  v  is  the  frequency  of  emitted  radiation  and  ao  is  the  Bohr  radius. 
The  matrix  elements. 


^  max  {%,%')  _2_  2  L  -'n'  _  J 
3  -JTT-I ®  r  P  ^  r  pi  >1 


I  r|  pi  >  = 

o 

where  R(p,  i)  represents  the  hydrogen  radial  wave  function,  are  listed 
by  Green,  Rush,  and  Chandler  (Ref.  30)  for  angular  momentum  values 
through  a  principal  quantum  number  of  20.  Detailed  results  of  these 
calculations  will  not  be  given  because  of  the  large  bulk  of  the  data  and 
the  relative  availability  of  the  values  in  the  literature.  The  detailed 
transitions  among  the  upper  levels  are  relatively  unimportant  for 
this  investigation  because  the  conditions  will  be  restricted  to  those 


/ 


R(p',£'')  r  R(p,£)  r  dr 


26 


AEDC-TR-76-5 


cases  for  which  the  plasma  is  strongly  collision  dominated  among  the 
upper  levels.  Under  these  conditions,  the  population  density  of  the 
angular  momentum  states  within  a  principal  quantum  level  can  be  char¬ 
acterized  by  a  Maxwell -Boltzmann  distribution  and  the  entire  set  of 
states  is  thus  adequately  described  by  a  single  state  which  has  proper¬ 
ties  represented  by  an  average  of  the  detailed  properties  of  the  indi¬ 
vidual  substates.  In  low  density  applications  where  radiative  terms 
are  important,  the  approximations  made  here  would  not  necessarily  be 
valid,  and  care  must  be  exercised  in  applying  the  hydrogenic 
approximation. 

3.2  RATE  COEFFICIENTS 


As  mentioned  earlier,  the  rate  coefficients  are  related  to  the 
probability  that  a  particular  transition  will  occur  upon  collision.  It  is 
in  fact  the  probability  that  a  coUisional  transition  will  occur  when  an 
atom  encounters  an  electron  with  kinetic  energy  between  E  and  E+dE 
integrated  over  the  kinetic  pnergy  distribution  of  the  electrons. 


If  fe(E^)  is  the  fractional  number  of  electrons  with  energy  in  the 
range  E  to  E+dE,  Vg  is  the  kinetic  velocity  of  these  electrons,  and  n^ 
is  the  concentration  pf  these  electrons  on  a  volume  basis,  then  the 
product  ngVgfg(E)  is  the  number  of  encovuiters  between  the  electrons 
and  a  hypothetical  unit  area  ppr  unit  tinie  along  a  velocity  vector 
traveled  by  the  electron.  Now  if  Q(E)  is  the  effective  area  the  atom 
presents  to  the  electrons  in  this  energy  range  for  a  particular  transi¬ 
tion,  the  number  of  encpimters  with  thie  area  per  unit  time  is  given 
simply  by  the  product 

Q(E)  ngVgf^(E) 

Hence,  if  the  concentration  of  atoms  with  cross  section  Q(E)  is  n, 
the  total  time  rate  of  change  of  atoms  due  to  collisions  with  electrons 
with  all  energies  is 


an  T 

at  =  “  "e  / 

o 


Q(E)  f  (E)  V  d  E 
e  e 


Defining  the  rate  coefficient  as  the  integral 


K  = 


} 


Q(E)  fg(E)  Vg  d  E 
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noting  that 


E 


2 


and  requiring  the  electron  distribution  to  be  Maxwell-Boltzmann, 


%  (ki 


3/2 


^  exp  (-E/kT) 


yields 


K  = 


Stt 

e 


3/2 


2irkT 


Q(E)  exp  (-E/kT)  E  dE 


for  the  rate  coefficient. 


The  development  to  this  point  has  been  general,  requiring  only  the 
assumption  of  Maxwell-Boltzmann  statistics  for  the  electron  energy 
distribution  and  a  delta  function  distribution  in  energy  space  for  the  atoms. 
Hence,  if  one  has  Q(E)  for  each  of  the  various  processes,  the  evaluation 
of  each  of  the  rate  coefficients  will  proceed  straightforwardly.  It  is  this 
evaluation  of  Q(E)  to  which  subsequent  discussion  is  turned. 


3.2.1  Two-Body  Internal  Transition  Rate  Coefficient,  K(p,x|) 


The  two-body  internal  transition  rate  coefficient  is  given  by 

m 

3/2 


K{p;q)  = 


M8tt 


1 

l27rkT 


I 


Q  (E)  exp  (-E/kT)EdE  (9) 


Ie  -e 

I  p  q 


In  the  above  equation,  p  and  q  symbolically  represent  sets  of  quantum 
numbers  describing  initial  (p)  and  final  (q)  states,  M  is  a  parameter 
giving  the  number  of  electrons  available  for  excitation  (M  =  2  for  p  = 
ground  state,  m  =  1  for  p  >  ground  state),  Qpq  (E)  is  the  energy- 
dependent  cross  section  for  collisional  transitions  from  the  p^^^  to 
the  q^^  bound  states,  E  is  the  kinetic  energy  of  the  bombarding 
electron,  and  Ep  and  Eq  are  kinetic  energies  of  the  bound  electron 
in  the  p  and  q  states. 
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There  are  many  techniques  for  the  evaluation  of  Qpq  (E)  for  tran¬ 
sitions  between  bound  states  utilizing  various  approximations.  Moisei- 
witsch  and  Smith  (Ref.  31)  give  a  comprehensive  review  of  the  various 
techniques  and  results  from  both  theoretical  calculations  and  experimen¬ 
tal  measurements  available  in  1968.  The  most  extensive  tabulation  they 
give  for  helium  excitation  was  first  calculated  by  Ochkur  and  Bratsev 
(Ref.  32)  and  is  based  upon  quantum  mechanical  considerations. 

Rather  than  use  the  more  satisfying  quantum  mechanical  calculations, 
however,  it  was  decided  to  use  the  simple  classical  Gryzinski  method 
(Ref.  33)  for  these  determinations.  This  choice  is  motivated  primarily 
by  the  ease  with  which  these  cross  sections  may  be  evaluated,  the 
observation  that  most  previous  investigators  of  CRR  also  used  the 
Gryzinski  method,  and  the  fact  that  the  method  appears  to  be  as  accvirate 
as  other  techniques. 

Gryzinski’ s  method  is  based  on  a  classical  analysis  describing  the 
collision  as  a  two-body  coulomb  interaction  between  the  orbital  and 
incident  electrons.  In  earlier  work  (Ref.  34)  he  assumed  that  the  orbital 
electrons  had  a  certain  fixed  energy.  In  later  work  (Refs.  33,  35,  and 
36)  this  was  extended  to  the  specific  case  of  atomic  collisions  in  which 
the  bound  electron  has  a  velocity  distribution 

*v  =  ( <-''/%> 
where  Vjjj  is  the  mean  electronic  velocity. 

The  choice  for  this  distribution  appears  to  be  arbitrary,  being  based 
on  phenomenological  observation  of  residts  rather  than  developed  from 
first  principles.  However,  the  relative  success  of  the  technique  is  suffi¬ 
cient  to  justify  its  use. 


The  cross  sections  obtained  by  the  Gryzinski  technique  yield  a  result 
for  a  transfer  of  energy  equal  to  or  greater  than  some  value  Upq  from  ' 
the  incident  electron  to  the  bound  electron.  Expressed  in  terms  of  a  para¬ 
meter  X,  the  Gryzinski  cross  section  is  written 


(X) 


TTe 


u 


_  A 
2  X 


pq 


X 

A+x 


3/2 


x-1 

X 


2A+1 

A+1 


(10) 


29 


AEDC-TR-76-5 


where 


X 

A 


E/U. 
p'  'pq 


pq' 
E^/U. 


In  the  above,  E  is  the  kinetic  energy  of  incident  electron,  Ep  is  the 
kinetic  energy  of  bound  electron,  e  is  the  electronic  charge,  and  e  is  the 
natural  logarithm  base. 

Identifying  Upq  with  the  threshold  energy  for  excitation  to  the 
state  from  the  state,  the  expression  above  gives  the  cross  section 
for  an  energy  transfer  which  not  only  is  sufficient  to  excite  the  bound 
electron  to  the  state  Eq,  but  may  be  large  enough  to  excite  to  the  state 
Eq+i.  Hence,  in  the  determination  of  Qpq(x),  it  is  necessary  to  form 
the  difference: 


Q  (x)  =  Q 

pq  E 


1  E  I  _ 

E 

Iv)  Wi 

1  ^pq+1  . 

(11) 


Note  that  Q pq+ltE/Up^^j)  is  vanishingly  small  for  E  <  Upq+j- 

Using  the  parametric  form  of  Q  (x),  the  rate  coefficient  K(p;q)  is 

fjq 


K(p;q)  = 
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U  xdx  (12) 


Thus  far  there  have  been  no  applications  of  the  results  of  quantum 
mechanics.  The  above  discussion  centered  on  a  straightforward  energy 
transfer  from  an  incident  classical  electron  impinging  upon  an  orbital 
electron,  the  orbital  electron  receiving  enough  energy  from  the  incident 
electron  to  raise  it  to  some  higher  energy  state.  Now  from  quantum 
mechanics  it  is  known  that  some  transitions  will  be  allowed  and  that 
other  transitions  are  forbidden.  Hence,  the  dipole  selection  rules  AL  = 
±1,  AS  =  0  must  be  applied  so  that  Eqs.  (9)  and  (11)  are  used  only  for 
those  transitions  between  energy  states  satisfying  these  rules. 

Another  process  in  these  excitation  collisions  which  can  occur  is 
that  of  an  exchange  collision  in  which  the  orbital  electron  is  ejected  and 
the  impinging  electron  is  captured  by  the  nucleus  in  some  higher  energy 
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state.  For  these  collisions,  Gryzinski's  development  assumes  that  the 
orbital  electron  has  zero  velocity  and  results  in  the  parametric  form 


pq 
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ire 
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pq 


X“1  ^  „ 

A (x+A)  (x+A-1) '  ^  ^ 

B-1  ..  .  „ 
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(13) 


•where 


and 
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A  = 


B  =  U  ^,/U 

pq+1  pq 


U  .,  =  E  -  E  . . 
pq+1  p  q+1 


The  factor  of  1/2  in  the  expression  is  the  probability  that  an  imping¬ 
ing  electron  has  the  proper  spin  orientation  to  effect  the  transition.  The 
rate  coefficient  will  have  the  sanie  functional  form  as  Eq.  (12),  the  only 
real  difference  being  the  functional  form  of  the  cross  section  (Eq.  (13) 
as  opposed  to  Eq.  (10))  and  the  restriction  of  its  application  to  those 
transitions  whose  initial  and  final  states  violate  the  dipole  selection  rules 
AL  =  ±1,  AS  =  0. 

There  is  some  judgment  required  in  determining  precisely  what 
state  to  use  for  Up^+i.  For  those  transitions  which  are  dipole  allowed, 
it  is  obvious  that  the  state  to  use  for  Upn+i  is  the  next  higher  state 
which  will  also  satisfy  the  selection  rules.  For  dipole -forbidden  transi¬ 
tions,  however,  the  choice  is  not  quite  so  obvious.  The  exchange  cross 
section  is  applied  to  any  transition  which  is  dipole  forbidden.  This  means 
any  transition  which  has  AL  ^  ±1  or  AS  0.  For  a  particular  tran¬ 
sition,  there  is  a  specific  AL  or  AS.  This  implies  specific  conditions  on 
the  impinging  electron.  The  next  higher  excitable  state  will  be  that  one 
which  is  compatible  with  the  AL  or  AS  involved  in  the  specific  excitation  of 
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•interest.  Hence,  the  q  +  1st  state  will  be  the  next  higher  state  with 
the  same  AL  or  AS  value  as  the  state.  Dugan  and  Sovie  (Ref.  37) 
examine  the  behavior  of  the  Gryz inski  cross  sections  with  various  choices 
of  Upq+i  and  this  choice  is  in  accord  with  their  findings. 

Because  of  the  large  bulk  of  the  data  involved  in  these  calculations, 
a  detailed  accounting  of  the  values  for  the  rate  coefficients  will  not  be 
given  here  but  rather  the  results  of  using  the  cross  section  expression 
applied  to  selected  transitions  will  be  discussed. 

Figures  1  and  2  illustrate  the  results  of  the  Gryz  inski  cross  section 
calculations  as  well  as  calculations  based  on  the  Ochkur  approximation 
compiled  in  Ref.  31  and  experimental  data  (Refs.  38  and  39)  for  the 
llS  -  3lP  and  the  l^S  -  3^8  transitions,  respectively.  Results  of 
ground  state  excitation  are  assumed  to  be  representative  for  the  other 
excitations.  A  scarcity  of  data  in  the  literature  prevents  comparisons, 
of  excitations  among  the  higher  states.  As  can  be  seen,  the  calculation 
based  upon  the  Ochkur  approximation  gives  more  satisfying  comparison 
with  experimental  data  for  both  the  direct  excitation  and  the  exchange 
excitation.  However,  use  of  the  Ochkur  approximation  is  in  itself  an 
adaptation  of  the  Born-Oppenheimer  method  and  is  thus  theoretically 
applicable  only  for  the  high-energy  collisions.  The  electronic  kinetic 


Figure  1.  P  S  -  S'*  P  excitation  cross  section. 
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energies  with  which  this  study  is  concerned  will  peak  at  about  1  ev  or  less 
so  that  the  low  energy  region  of  the  cross  section  is  the  most  important. 

In  examining  the  lower  energy  region  of  Fig.  1  it  can  be  seen  that 
the  Gryzinski  cross  section  will  overestimate  the  experimental  work  of 
St.  John,  et  al.  (Ref.  38)  by  at  most  about  a  factor  of  two  in  the  energies 
close  to  threshold  (x  =  1).  It  is  difficult  to  be  more  precise  in  these 
comparisons  because  of  the  scarcity  of  experimental  work  in  these  re¬ 
gions  close  to  threshold.  Calculations  of  the  rate  coefficient  for  the 
2^8  -  transition  using  the  Gryzinski  technique  and  the  Ochkur  approx¬ 
imation  for  cross  sections  yield  results  which  agree  well  within  a  factor 
of  two.  It  is  thus  reasonable  to  expect  the  calculated  rate  coefficient  to 
overestimate  the  true  rate  coefficient  by  no  more  than  a  factor  of  two 
and  in  fact  probably  to  be  better  than  this  because  of  the  underestimation 
of  the  cross  section  at  higher  energies  by  the  Gryzinski  method. 

In  examining  the  exchange  cross  section  (Fig.  2)  one  can  see  that  the 
comparison  between  calculation  and  experiment  is  not  as  good  as  that 
for  direct  excitation,  although  there  does  appear  to  be  a  convergence  of 
sorts  between  the  Ochkur  approximation  and  the  later  experimental  work 
of  Moussa,  et  al.  (Ref.  39).  The  Gryzinski  cross  section  overestimates 
the  other  work  by  several  orders  of  magnitude  at  the  sharp  peak  and 
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underestimates  the  work  of  Moussa  by  typically  a  factor  of  two  at  the  high 
energies.  Hence,  because  of  the  sharp  peak  in  the  region  of  lowest  energy, 
it  is  anticipated  that  the  calculated  rate  coefficient  will  underestimate  the 
true  rate  coefficient  by  somewhat  less  than  the  factor  of  two.  At  any  rate 
the  exchange  cross  sections  will  be  somewhat  less  important  than  the 
cross  sections  for  direct  excitation  because  the  rate  coefficients  for  direct 
excitation  are  much  larger  than  those  for  exchange  excitation.  For  example, 
at  an  electronic  kinetic  temperature  of  10,000®K,  the  rate  coefficient  for 
the  2lS  -  sip  direct  excitation  is  about  10"  1®  cm^/sec,  whereas  for  the 
exchange  excitation  to  the  adjacent  S^D  level  (2ls  -  3^D)  the  rate  coeffi¬ 
cient  is  about  10"  13  cm3 /sec.  This  is  not  necessarily  true  for  transitions 
among  sublevels  of  the  higher  excited  states  where  the  energy  levels  are 
quite  dose  together.  However,  these  upper  level  rate  coefficients  will 
be  large  and  the  popdation  densities  will  rapidly  establish  an  equilibrium 
configviration  among  the  sublevels.  Requiring  a  detailed  accounting  of 
each  sublevel  is  unnecessary  for  all  but  the  earliest  stages  of  the 
computation. 

The  preceding  discussion  has  been  concerned  only  with  inelastic  ex¬ 
citation  between  bound  states  of  the  atom.  The  inverse  process,  a  super¬ 
elastic  collision  resulting  in  de -excitation  of  the  atom,  is  also  included 
in  this  study.  The  rate  coefficients  for  these  processes  come  quite  easi¬ 
ly  from  equilibrium  considerations. 

At  collisional  equilibrium  between  two  states,  the  rates  of  change  of 
the  two  states  will  be  equal.  Further,  detailed  balance  requires  that 
rates  due  to  individual  processes  be  equal.  Hence,  one  can  set 

ne^pK(p;q)  -  ngngK(q;p)  =  0 

Thus, 

n 

^  K(p;q) 

q 

and,  since  this  is  an  equilibrium  condition, 

^  exp  [-(e  -E  VkT  1 
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Hence, 
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3.2.2  Collisional  Ionization  and  Three-Body  Recombination 
Coefficients,  K(p;c}  and  K(c;p} 


The  collisional  ionization  cross  section  is  very  readily  obtained  from 
the  two -body  internal  transition  cross  section.  This  is  accomplished  by 
noting  that  the  Gryzinski  formulation  for  the  cross  section,  Eq.  (9),  is 
inversely  dependent  upon  the  square  of  energy  difference  between  the 
initial  and  final  states,  Upq.  For  the  final  state  being  the  continuum, 
this  Upq  just  becomes  the  ionization  potential  for  state  p.  Hence,  the 
parametric  cross  section  for  collisional  ionization  Qpc(x)  is 

1"  <«  +  ''5c=T)j  ' 


TTe  1 


Sil'I 


>•1 


X  =  E/1 


and  Ip^  is  the  ionization  potential  of  the  p'*^^  quantum  level.  Hence  the 
collisional  ionization  rate  coefficient  is  written 


K(p;c)  = 
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(14) 


In  this  case  again,  M  =  2  for  the  ground  state  atoms  since  there  are 
two  electrons  available  for  ionization  from  the  ground  state  and  M  =  1 
for  atoms  in  higher  energy  states. 


Figures  3  and  4  show  results  of  the  Gryzinski  calculations  for  the 
ionization  of  the  llS  and  2^S  levels,  respectively,  compared  to  experi¬ 
mental  results  reported  in  the  literature  (Refs.  40,41,42,  and  43).  The 
Gryzinski  cross  section  overestimates  the  ionization  cross  section  for 
the  groimd  state  atom  by  about  25  percent  compared  to  the  experimental 
data.  Further,  the  peak  of  the  calculated  cross  section  occurs  some¬ 
what  earlier  than  the  experimental  peak  but  for  the  purpose  of  this  study 
the  calculations  can  be  considered  to  give  satisfactory  results. 

The  calculated  cross  section  for  the  2^8  ionization  does  not  agree  as 
satisfactorily  with  experiment  as  is  shown  in  Fig.  4.  For  this  the  data 
of  Long  and  Geballe  (Ref.  42)  are  to  be  considered  superioirto  those  of 
Fite  and  Brackmann  (Ref.  43)  since  the  former  data  were  taken  from  a 
source  known  to  be  the  23s  atom,  whereas  Fite  and  Brackmann’s  source 
was  an  indeterminant  mixture  of  the  2^8  and  2^8  atoms. 
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For  the  purposes  of  this  study,  the  agreement  is  considered  satis¬ 
factory.  It  is  expected  that  the  characteristics  of  the  ionization  cross 
section  calcvQations  will  follow  the  same  pattern  for  the  higher  excited 
states. 

The  inverse  of  coUisional  ionization,  three-body  recombination,  is 
described  by  the  rate  coefficient  K(c;p).  This  coefficient  is  also  obtained 
from  equilibrium  considerations.  When  equilibrium  between  state  p  and 
the  continuum  is  established  by  these  two  processes,  the  rate  of  filling 
state  p  due  to  three -body  recombination  exactly  balances  the  rate  the  state 
is  depopulated  via  coUisional  ionization,  or 

2 

ngnpK(p,c)  =  n^  n^K(c,p) 


36 


AEDC-TR-76-5 


Figure  4.  2^8  ionization  cross  section. 


where  m  is  the  population  density  of  positive  ions.  Hence, 

K(c;p)  =  K{p;c) 

e  i 


Since  only  single  ionization  is  considered  here,  n^  =  nj.  Further, 
this  is  an  equilibrium  relationship  between  state  p  and  the  continuum  so 
that  Up  and  Ug  are  related  the  modified  Saha  equilibrium  relationship 
so  that 


=7:; - 7572  5^  K(p;c)  (15) 
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3.2.3  Radiative  Recombination  Rate  Coefficients,  j3(p) 

The  radiative  recombination  rate  coefficient  describes  the  filling 
of  state  p  due  to  two-body  encounters  between  electrons  and  ions  in  which 
the  ion  captures  the  electron  into  energy  state  p  with  the  attendant  emis¬ 
sion  of  radiation.  It  is  obtained  in  much  the  same  way  as  the  other  rate 
coefficients  although  there  is  some  difference  in  the  physical  analysis. 
The  basic  functional  form  for  the  rate  coefficient  is  written  as  before. 
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where  j3(p)  is  the  radiative  recombination  rate  coefficient  for  state  p, 
and  Q2.(p)  is  the  cross  section  for  the  capture. 

The  evaluation  of  Q^Cp)  does  not  proceed  as  straightforwardly  here, 
however.  To  obtain  this  cross  section  it  is  necessary  to  use  the  prin¬ 
ciple  of  detailed  balance  for  the  relationship  between  the  radiative 
recombination  cross  section  and  the  photoionization  cross  section  (or 
absorption  coefficient)  thus  (Ref.  44), 


Q^tP) 


(hv)^ 

2m  E 
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a^(p) 
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where  gp  and  g^on  are  the  statistical  weights  of  the  atom  and  ion,  res¬ 
pectively,  E  is  the  kinetic  energy  of  the  electron,  hi/  is  the  photon 
energy  emitted  (or  absorbed)  in  a  photoelectric  transition  between  the 
continuum  and  the  state  p,  and  a^  (p)  is  the  absorption  coefficient  for 
state  p. 

Noting  that  the  energy  of  an  emitted  photon  in  the  capture  will  include 
both  the  ionization  potential  of  state  p  and  the  kinetic  energy  of  the 
electron  before  capture. 


hv  =  I  +  E 

poo 
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Substituting,  one  obtains  with  this  change  in  variable 

00 

e(p)  =  ^  |l) ^ 

I 

P 

exp  (-hv/kTg)  a^(p)  d{hv)  (18) 


The  photoionization  cross  section  is  given  for  hydrogenic  atoms  by  Sea¬ 
ton  (Ref.  45); 


a^Cp) 


64a7Ta 

o 


2 


3/3 


gii(p/u) 


(19) 


where  a  is  the  fine  structure  constant,  ao  is  the  Bohr  radius,  Z  is 
the  nuclear  charge,  u  is  the  ratio  E/Ip®,  p  is  the  principal  quantum 
number,  and  gjj  is  the  Gaunt  factor; 


’II 


=  1  +  0.1728 


-  0.0496 


u  +4/3  u+1 


<20) 


Noting  that 

hv  =  I  (1  +  u) 

poo 


and  substituting. 


B(p)  =  B(p) 


9ll 


where 


B(p) 


64/TTa _ ^o 

3/3  (mekTe)^/^  ^ion 


E_ 


(21) 


(22) 
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The  expression  is  based  upon  a  hydrogenic  approximation  in  -which  p 
is  the  principal  quantum  number  of  the  energy  level.  Because  of  this 
hydrogenic  approximation,  one  would  expect  the  calculations  to  give  their 
worst  results  for  the  lower  quantum  levels  and  to  converge  to  correct 
values  for  the  upper  levels.  Figure  5  shows  the  results  of  the  calcvila- 
tion  of  the  absorption  cross  section  for  the  ground  state  helium  atom 
compared  to  the  experimental  results  of  Lowry,  Tomboulian,  and  Ederer 
(Ref.  46).  In  this  figure,  in  which  the  abcissa  is  wavelength,  energy 
increases  to  the  left,  threshold  energy  (electronic  velocity  of  0)  occurring 
at  about  505  A.  In  this  case,  a  1-ev  electron  wovild  be  ejected  at  an 
abcissa  very  near  this  threshold  wavelength.  Thus,  for  the  rate  coeffi¬ 
cient  calculation,  the  rightmost  portion  of  the  curve  is  most  important, 
and  as  can  be  seen,  the  calculations  reproduce  the  experiment  in  this 
region  to  within  a  factor  of  two  or  better.  Thus,  one  would  expect  the 
radiative  recombination  rate  coefficient  to  be  correct  to  within  a  factor 
of  two  or  better  for  the  ground  state  and  become  progressively  better 
for  the  higher  quantum  levels. 


Figure  5.  1^  S  absorption  cross  section. 
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4.0  TECHNIQUES 


This  chapter  summarizes  the  various  approximations  and  techniques 
which  were  utaized  in  obtaining  the  solution  to  the  ERE.  A  detailed 
listing  of  the  computer  programs  will  not  be  included  since  the  techniques 
are  quite  straightforward.  Rather,  a  short  discussion  of  programming 
considerations  will  be  given;  since  the  solution  to  the  ERE  is  obtained 
numerically,  certain  approximations  and  compromises  were  necessary 
and  these  are  discussed  in  detail, 

4.1  ENERGY  LEVELS 

The  helium  energy  levels  used  in  this  study  were  obtained  from 
Moore  (Ref.  47),  who  gives  all  the  energy  levels  from  which  spectral 
lines  have  been  observed.  The  precise  energy  values  of  all  the  detailed 
substates  considered  in  the  program  are  not  known,  or  are  so  closely 
adjacent  to  other  substates  of  the  same  quantum  level  that  they  have  not 
been  resolved.  The  energy  differences  between  these  states  will  be  so 
small  that  they  are  negligible. 

Evaluation  of  the  rate  coefficients  required  different  substates  to 
have  different  energies,  however.  This  is  because  in  the  course  of 
evaluation,  differences  in  energies  of  the  states  are  taken  and  a  zero 
energy  difference  causes  machine  errors  to  occur.  To  prevent  this 
condition  simple  linear  interpolation  was  used  to  generate  different 
numerical  energies  for  the  substates  not  given  explicitly  in  Moore  (Ref. 
47).  A  listing  of  the  energy  levels  used  in  this  study  is  given  in  Table  1, 

The  size  of  the  computer  core  dictated  that  not  all  energy  states 
listed  in  Table  1  be  used.  For, the  higher  principal  quantum  states  the 
distribution  of  densities  among  the  substates  will  be  maintained  in  a 
Boltzmann  configuration.  The  total  density  of  the  level  immediately 
gives  the  density  of  any  desired  substate. 

To  meet  the  computer  core  requirement,  the  a  priori  assumption  of 
the  Boltzmann  configuration  among  the  substates  of  quantum  levels  7  and 
8  was  made.  These  quantum  levels  were  described  by  a  single  energy 
state  chosen  to  be  the  lowest  energy  exhibited  by  any  of  the  substates  of 
that  level.  This  assumption  was  used  only  for  the  final  determination  of 
states  considered  by  the  ERE.  The  rate  coefficients  were  determined 
using  detailed  energies  of  the  substates  of  these  quantum  levels,  and  av¬ 
erage  rate  coefficients  for  the  total  quantum  level  were  then  determined 
from  these  values. 
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Table  1.  Helium  Energy  Levels 


state 

g 

E  (1/cm) 

State 

g 

E  (1/cm) 

l^S 

1 

0.0 

5^G 

27 

193921.0 

2^S 

3 

159850.3 

5^p 

3 

193936.8 

2^S 

1 

166271.7 

6^S 

3 

194930.5 

2^P 

9 

169081.1 

6^S 

1 

195109.2 

2^P 

3 

171129.2 

6^P 

9 

195187.2 

3^S 

3 

183231.1 

6^D 

15 

195254.4 

3^3 

1 

184859.1 

6^D 

5 

195255.0 

3^P 

9 

185558.9 

6^ 

7 

195256.7 

3^D 

15 

186095.9 

6^F 

21 

195256.3 

3^D 

5 

186099.2 

6^G 

9 

195258.0 

3^P 

3 

186203.6 

6^G 

27 

195258.5 

4^S 

3 

190292.5 

6^H 

11 

195260.0 

4^S 

1 

190934.5 

6^H 

33 

195260.5 

4^ 

9 

191211.4 

6^P 

3 

195269.2 

4^D 

15 

191438.8 

7^S 

3 

195862.6 

4^D 

5 

191440.7 

7^S 

1 

195973.2 

4^F 

21 

191446.6 

7^P 

9 

196021.7 

4^F 

7 

191447.2 

7^D 

15 

196064.0 

4^P 

3 

191487.0 

7^D 

5 

196064.3 

5^S 

3 

193341.3 

7^F 

7 

196065.4 

5^S 

1 

193657.8 

7^F 

21 

196065.5 

5^P 

9 

193795.1 

7^G 

9 

196066.5 

5^D 

15 

193911.5 

7^G 

27 

196066.6 

5^D 

5 

193912.5 

7^H 

11 

196067.5 

5^F 

7 

193914.3 

7^H 

33 

196067.6 

5^F 

21 

193915.8 

7^1 

13 

196068.5 

5^G 

9 

193920.0 

7^1 

39 

196068.6 

7^P 

3 

196073.4 

B^H 

11 

196591.6 

B^S 

3 

196455.8 

B^H 

33 

196591.7 

8^S 

1 

196529.0 

8^1 

13 

196592.2 

8^P 

9 

196561.1 

8^1 

39 

196592.3 

B^D 

15 

196589.4 

8^J 

IS 

196592.8 

B^D 

5 

196589.7 

B^J 

45 

196592.9 

B^F 

7 

196590.3 

B^P 

3 

196595.6 

B^F 

21 

196590.4 

9 

324 

196856.4 

B^G 

9 

196591.0 

10 

400 

197139.8 

B^G 

27 

196591.1 

11-15 

3420 

197347.0 

/ 

42 


AEDC-TR-76-5 


Quantum  levels  9  through  15  have  very  small  energy  widths  and  a 
progressively  larger  number  of  substates  with  each  principal  quantum 
level.  Consequently  each  was  considered  as  a  single  energy  level  for  all 
rate  coefficient  determinations. 


The  criterion  for  truncating  the  system  of  equations  comprising  the 
ERE  is  that  some  upper  level  will  be  in  equilibrium  with  the  free  electron 
density.  In  order  to  maintain  a  correct  description  of  the  rate  of  this 
upper  level,  states  above  it  must  also  be  considered  in  order  to  have 
approximately  correct  values  for  collisional  and  radiative  processes. 
Quantum  level  10  was  chosen  to  be  the  maximum  quantum  level  for  the 
critical  level  for  this  study.  Since  quantum  levels  11  through  15  will 
therefore  be  in  an  equilibrium  configuration,  they  were  averaged  together 
to  form  a  reservoir  for  the  various  processes  for  quantum  level  10. 


An  additional  phenomenon  which  relates  to  the  energy  levels  in 
partially  ionized  plasmas  and  must  be  taken  into  account  is  that  of  the 
lowering  of  the  ionization  potential.  This  arises  because  when  a  bound 
electron  is  in  a  sufficiently  high  energy  state  it  can  exceed  the  Debye 
shielding  length  for  the  plasma  and  thus  becomes  indistinguishable  from 
the  free  electrons.  However,  at  typical  conditions  encountered  in  this 
study,  ng  =  10^®  l/cm^,  T^  =  10^  °K,  the  expression  for  the  lowering 
of  the  ionization  potential  given  by  J.  Richter  {Ref.  48)  gives  a  value 
of  74  1/cm,  which  is  insignificant. 


In  summary,  plasma  conditions  are  restricted  to  such  that  at  least 
quantum  level  10  is  in  Saha  equilibrium  with  the  free  electrons.  A 
detailed  accounting  of  all  singlet  and  triplet  states  is  maintained  through 
quantum  level  6.  Principal  quantum  levels  7  through  10  are  maintained 
separate  from  each  other  although  a  Maxwell -Boltzmann  density  distri¬ 
bution  of  the  substates  of  each  of  these  principal  quantum  levels  is 
assumed. 


4.2  PROGRAMMING  CONSIDERATIONS 

As  indicated  above,  one  important  consideration  in  the  development 
of  the  computer  program  is  available  computer  core  size.  Arrays  for 
each  of  the  rate  coefficients  must  be  maintained  as  well  as  a  work  area. 
Utilizing  a  two-dimensional  array  for  the  K(p;q)  and  A(p;q)  with  the 
other  parameters  necessary  in  the  solution,  47  sepEirate  energy  states 
could  be  considered  in  approximately  116,000  bytes  of  core  in  the  IBM 
S360/50  computer. at  the  Arnold  Engineering  Development  Center. 
Although  more  core  was  available  at  the  time  of  this  work,  scheduling  qf 
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Figure  7.  n(2''  P)/n(Saha)  versus  time  for  two  error  criteria. 

error  bounds.  These  calculations  were  performed  for  precisely  the 
same  plasma.  Case  C  to  be  presented  later,  with  the  error  bound  the 
only  difference.  The  first  of  these,  10"®,  is  typical  of  the  error  bounds 
used  in  the  other  studies.  The  other  error  bound,  10“12  for  t  <  10"^®  sec 
and  10"^^  for  t  >  10"^®  sec,  is  an  extremely  tight  bound.  Some  experi¬ 
mentation  was  conducted  at  significantly  lower  error  bounds  and  it  was 
found  that  the  solution  would  begin  to  oscillate  if  the  bound  was  relaxed 
too  much,  although  a  subsequeni  tightening  of  the  error  bound  would  cause 
the  solution  to  become  stable  again. 

The  repeatability  of  the  solution  under  the  two  bounds  indicated  above 
illustrates  the  stability  of  the  solution  at  the  10"®  bound  which  was  used 
dtiring  the  early  part  of  the  solution.  It  was  foimd  that  after  the  early 
transients  had  died  out  the  error  bound  could  be  relaxed  substantially 
with  no  loss  of  stability.  This  stability  along  with  numerical  agree¬ 
ment  with  other  independent  calculations  of  the  CRR  coefficients,  as  will 
be  indicated  in  Section  5. 0,  shows  that  the  numerical  technique  used  in 
these  studies  yields  valid  results. 
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Quantum  levels  9  through  15  have  very  small  energy  widths  and  a 
progressively  larger  number  of  substates  with  each  principal  quantum 
level.  Consequently  each  was  considered  as  a  single  energy  level  for  all 
rate  coefficient  determinations. 

The  criterion  for  truncating  the  system  of  equations  comprising  the 
ERE  is  that  some  upper  level  will  be  in  equilibrium  with  the  free  electron 
density.  In  order  to  maintain  a  correct  description  of  the  rate  of  this 
upper  level,  states  above  it  must  also  be  considered  in  order  to  have 
approximately  correct  values  for  collisional  and  radiative  processes. 
Quantum  level  10  was  chosen  to  be  the  maximum  quantum  level  for  the 
critical  level  for  this  study.  Since  quantum  levels  11  through  15  will 
therefore  be  in  an  equilibrium  configuration,  they  were  averaged  together 
to  form  a  reservoir  for  the  various  processes  for  quantum  level  10. 

An  additional  phenomenon  which  relates  to  the  energy  levels  in 
partially  ionized  plasmas  and  must  be  taken  into  account  is  that  of  the 
lowering  of  the  ionization  potential.  This  arises  because  when  a  bound 
electron  is  in  a  sufficiently  high  energy  state  it  can  exceed  the  Debye 
shielding  length  for  the  plasma  and  thus  becomes  indistinguishable  from 
the  free  electrons.  However,  at  typical  conditions  encountered  in  this 
study,  Ug  =  10^^  1/cm^,  Tg  =  10^  ®K,  the  expression  for  the  lowering 
of  the  ionization  potential  given  by  J.  Richter  (Ref.  48)  gives  a  value 
of  74  1/cm,  which  is  insignificant. 

In  summary,  plasma  conditions  are  restricted  to  such  that  at  least 
quantum  level  10  is  in  Saha  equilibrium  with  the  free  electrons.  A 
detailed  accounting  of  all  singlet  and  triplet  states  is  maintained  through 
quantum  level  6.  Principal  quantum  levels  7  through  10  are  maintained 
separate  from  each  other  although  a  Maxwell-Boltzmann  density  distri¬ 
bution  of  the  substates  of  each  of  these  principal  quantum  levels  is 
assumed. 


4.2  PROGRAMMING  CONSIDERATIONS 

As  indicated  above,  one  important  consideration  in  the  development 
of  the  computer  program  is  available  computer  core  size.  Arrays  for 
each  of  the  rate  coefficients  must  be  maintained  as  well  as  a  work  area. 
Utilizing  a  two-dimensional  array  for  the  K(p;q)  and  A(p;q)  with  the 
other  parameters  necessary  in  the  solution,  47  separate  energy  states 
could  be  considered  in  approximately  116,  000  bytes  of  core  in  the  IBM 
S360/50  computer. at  the  Arnold  Engineering  Development  Center. 
Although  more  core  was  available  at  the  time  of  this  work,  scheduling  of 
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core  partitions  dictated  that  the  system  requirements  be  kept  within  133K 
bytes  to  expedite  the  work.  Extension  of  the  program  to  include  signifi¬ 
cantly  more  detail  would  have  caused  the  core  requirements  to  exceed 
the  optimum  13  3K  b3d:es. 

Also,  because  of  available  core,  it  was  impractical  to  calculate  the 
rate  coefficients  in  the  same  program  that  obtained  the  solution  to  the 
ERE.  Consequently,  these  calculations  were  accomplished  in  a  separate 
program  and  rate  coefficient  data  were  preserved  on  magnetic  tape  in 
binary  form  to  expedite  I/O  overhead. 

The  rate  coefficients  requiring  integration  are  evaluated  using  32- 
point  Laguerre -Gauss  quadrature  and  double  precision  arithmetic.  This 
is  a  standard  technique  which  is  readily  available;  details  of  the  integra¬ 
tion  are  included  in  Appendix  A. 

The  solution  to  the  ERE  was  obtained  using  the  modified  Euler's 
technique,  which  is  an  implicit  method.  Details  of  the  application  of  this 
technique  to  the  ERE  are  also  included  in  Appendix  A.  A  result  of  the 
technique  is  control  of  the  increment  of  calculation  based  upon  the  desired 
relative  error  bound.  The  calculations  presented  here  used  various 
bounds,  usually  as  high  as  10~®  to  10"^  in  the  portion  of  the  solution  in 
which  adjustments  were  taking  place  rapidly.  After  the  rapid  transients 
had  died  out,  error  bounds  as  low  as  10"^  were  used. 

All  calculations  were  performed  with  double  precision  arithmetic  and 
with  the  program  compiled  with  the  IBM  Fortran  H-level  compiler  using 
fiill  optimization  so  as  to  generate  the  most  efficient  machine  code.  In 
addition,  variables  and  calciolations  were  arranged  in  storage  so  as  to 
take  as  full  an  advantage  as  possible  of  the  optimization.  To  conserve 
compiler  time,  which  is  a  minor  consideration  except  for  programs  which 
are  run  many  times,  as  this  one  was,  the  compiled  program  was  preserv¬ 
ed  in  machine  language  on  disk  pack  for  immediate  recall  at  execution 
time.  These  efforts  to  reduce  execution  time  realized  an  approximate 
reduction  of  25  to  35  percent  compared  to  an  earlier  version  of  the  pro¬ 
gram.  Even  with  the  efforts  to  reduce  program  execution  time,  the 
run  times  can  become  quite  long.  Typical  run  times  for  each  case  stud¬ 
ied  here  were  45  minutes  to  one  hour.  Subsequent  computer  system 
improvements  have  reduced  these  computer  run  times  typically  an  or¬ 
der  of  magnitude. 

Consequently,  to  conserve  total  machine  time  and  insure  that  a 
solution  that  was  invalid  for  some  reason  would  not  waste  an  inordinate 
amount  of  machine  time,  the  program  was  setup  to  preserve  important 
solution  parameters.  Thus  one  can  restart  the  calculations  at  the  point 
in  time  the  parameters  were  preserved. 
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The  preservation  of  these  parameters  occurred  at  selected  clock 
time  intervals  and  plasma  time  intervals  so  that  an  inadvertent  machine 
or  program  failure  wovild  not  cause  the  loss  of  a  large  amount  of  compu¬ 
ter  time.  In  addition  the  parameters  were  preserved  at  the  last  calcu¬ 
lation  each  time  the  program  was  run  so  that  if  more  calculations  were 
necessary,  the  program  could  be  restarted  at  the  last  CEdcvilation  with  a 
minimum  of  overhead.  A  desirable  consequence  of  this  is  that  external 
judgment  control  can  be  exercised  over  the  execution  of  the  program. 
Hence,  the  error  bound  can  be  changed  as  the  solution  indicates  and  the 
number  of  calculations  judged  necessary  to  effect  the  solution  are 
under  external  control  so  that  unnecessary  calculations  are  minimized. 

4.3  PROGRAM  STABILITY 

As  indicated  above  and  shown  in  Appendix  A,  the  time  step  size  in 
the  solution  is  controlled  by  the  desired  error  bound.  Figures  6  and  7 
show  the  results  of  obtaining  a  solution  to  the  ERE  under  two  desired 


Time  (sec) 


Figure  6.  lie/nd)  versus  time  for  two  error  criteria. 
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na*  (sec) 


Figure  7.  n(2^  P)/n(Saha)  versus  time  for  two  error  criteria. 

error  bounds.  These  calculations  were  performed  for  precisely  the 
same  plasma.  Case  C  to  be  presented  later,  with  the  error  bound  the 
only  difference.  The  first  of  these,  10"^,  is  typical  of  the  error  bounds 
used  in  the  other  studies.  The  other  error  bound,  10 "^2  for  t  <  10“^^  sec 
and  10"^^  for  t  >  10"^®  sec,  is  an  extremely  tight  bound.  Some  experi¬ 
mentation  was  conducted  at  significantly  lower  error  bounds  and  it  was 
found  that  the  solution  would  begin  to  oscillate  if  the  boimd  was  relaxed 
too  much,  although  a  subsequent  tightening  of  the  error  bound  would  cause 
the  solution  to  become  stable  again. 


The  repeatability  of  the  solution  xmder  the  two  bounds  indicated  above 
illustrates  the  stability  of  the  solution  at  the  10"®  bound  which  was  used 
during  the  early  part  of  the  solution.  It  was  found  that  after  the  early 
transients  had  died  out  the  error  bound  could  be  relaxed  substantially 
with  no  loss  of  stability.  This  stability  along  with  numerical  agree¬ 
ment  with  other  independent  calculations  of  the  CRR  coefficients,  as  will 
be  indicated  in  Section  5.  0,  shows  that  the  numerical  technique  used  in 
these  studies  yields  valid  results. 
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5.0  PURE  AFTERGLOWS 

The  transient  solution  to  the  !ERK  under  the  restrictions  and  assump¬ 
tions  discussed  above  was  obtained  for  several  plasma  conditions.  This 
chapter  is  devoted  to  the  study  of  the  pure  afterglows  in  which  the  plasma 
was  allowed  to  relax  from  an  initial  distribution. 

5.1  REPORTING  OF  DATA 

Because  of  the  large  bulk  of  the  data  from  the  calculations,  the 
majority  of  the  data  will  be  presented  in  graphical  form  as  continuous 
functions.  Actual  points  from  the  calculations  will  not  be  shown  since  the 
lines  are  drawn  to  connect  the  points.  Additionally,  detailed  plots  from 
^  calculations  will  not  be  presented  but  rather  cases  will  be  selected  to 
illustrate  trends  of  the  plasma  decay.  The  cases  not  shown  in  the  detailed 
studies  wni  be  used  for  completeness  in  illustrating  certain  gross  trends 
of  the  decay  characteristics.  Tabular  information  and  precise  numbers 
from  the  calculations  will  be  used  only  as  needed  for  detailed  comparisons 

or  when  description  of  some  phenomenon  requires  resolution  beyond  that 
of  a  graph. 

5.2  PLASMA  CONDITIONS 

Transient  solutions  to  the  ERE  under  the  restrictions  and  assump¬ 
tions  discussed  above  were  obtained  for  several  afterglow  plasma  condi¬ 
tions  summarized  in  Table  2;  The  lower  limits  of  electron  temperature, 
degree  of  ionization,  or  total  number  density  represent  approximate 
lower  bounds  to  the  plasmas  which  could  be  solved  with  the  present  com¬ 
puter  program  because  of  the  critical  level  limits  described  earlier. 
Although  conditions  above  the  upper  limits  of  the  parameters  have  not 
been  investigated,  it  is  believed  that  extension  will  not  yield  notably 
different  data  on  the  mechanism  by  which  the  QSS  is  established. 

From  Table  2  one  can  see  that  cases  A  through  G  compare  plasmas 
at  different  total  densities  with  the  same  degree  of  ionization  and  electron 
temi^rature  and  cases  H  through  K  compare  plasmas  at  similar  selected 
densities  but  at  a  higher  electron  temperature,  14,000“K.  Cases  L 
through  O  compare  plasmas  at  constant  density  and  electron  temperature 
but  different  degrees  of  ionization.  Cases  C,  I,  P,  and  Q  compare 
plasmas  with  the  same  total  density  and  degree  of  ionization  but  different 
electron  temperatures.  By  quantitative  comparisons  of  the  results  for 
these  conditions,  insight  into  the  trends  of  the  processes  by  which  QSS 
is  established  will  be  obtained. 


47 


AEDC-TR-76-5 


Table  2.  Summary  of  Plasma  Conditions  for  Afterglows 
Studied  with  ERE 


Case 

"o 

1/cm^ 

"e/^o 

T 

e 

“K 

A 

17® 

1.63^' 

0.01 

1.0^ 

B 

4.4l« 

0.01 

1.0^ 

c 

I.5I® 

0.01 

1.0* 

D 

6.59^5 

0.01 

1.0* 

E 

5.0^^ 

0.01 

1.0* 

F 

1.5^^ 

0.01 

1.0* 

G 

5.0^^ 

0.01 

1.0* 

H 

1.63^’ 

0.01 

1.4* 

I 

i.s^e 

0.01 

1.4* 

J 

5.0^" 

0.01 

1.4* 

K 

5.0^^ 

0.01 

1.4* 

L 

2.69^^ 

0.2 

1.0* 

M 

2.69^^ 

0.01 

1.0* 

N 

2.69^^ 

0.001 

1.0* 

O 

2.69^^ 

0.0001 

1.0* 

P 

I.5I® 

0.01 

8.0^ 

Q 

1.5l« 

0.01 

6.0^ 

*Super scripts  denote  powers  of  10  by  which  the 
numbers  are  to  be  multiplied. 


The  initial  distribution  for  each  of  these  cases  was  chosen  such  that 
all  excited  state  densities  were  in  Saha  equilibrium  with  the  free  electron 
density.  Consequently,  the  excited  state  densities  exhibit  a  Boltzmann 
distribution.  In  most  laboratory  plasmas,  the  excited  state  distribution 
does  not  exhibit  this  Boltzmann  distribution  in  its  entirety.  Rather,  there 
can  be  significant  deviations  from  this  equilibrium  distribution  for  the 
first  few  excited  states.  In  this  sense  then,  the  initial  conditions  chosen 
for  the  afterglow  studies  are  at  variance  with  nature.  However,  the 
intent  is  to  examine  the  mechanisms  by  which  the  QSS  distribution  is 
established  and  thus  to  determine  the  important  phenomena  required  to 
maintain  the  QSS  in  an  actual  physical  situation.  In  a  physical  plasma 
various  of  the  plasma  parameters  will  be  changing  and  the  QSS  is  not 
a  static  condition.  By  studying  the  pure  afterglow  phenomena  with  the 
perturbational  studies  the  effect  upon  the  QSS  because  of  a  change  in 
the  plasma  parameters  can  be  examined  for  valid  application  to  physical 
plasmas. 
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5.2.1  Quasi-Steady-State  Recombination  Coefficients 


As  a  final  check  of  the  numerical  validity  of  the  computer  program 
the  collisional -radiative  recombination  coefficients  ce  were  obtained. 
These  results  were  obtained  by  allowing  the  transient  solution  to 
develop  until  rig  and  ri{l)  agreed  to  within  five  percent  and  then  cal¬ 
culating  or  with  the  instantaneous  values  of  rig  and  rig^,  a  =  -rig/ng^. 


Table  3  shows  the  value  of  a  obtained  from  the  ERE  and  results 
reported  by  Chen  (Ref.  21)  using  the  QSS  approximation.  As  is  seen 
in  Table  3,  the  present  resirits  compare  favorably  with  those  of  Chen, 
being  consistently  slightly  lower  except  for  cases  P  and  Q.  No  signifi¬ 
cance  is  placed  upon  the  variance  of  case  P,  being  less  than  ten  per 
cent.  Case  Q  will  be  discussed  in  more  detail  below.  Chen  also  used 
Gryzinski's  technique  for  obtaining  collision  cross  sections  but  used 
a  different  source  for  the  radiative  transition  probabilities.  Chen  did 
not  include  all  collisional  transitions,  as  this  work  does,  and  he  used 


Table  3.  Quasi-Steady-State  Collisional-Radiative  Recombination  Coefficients 
from  the  Eigenstate  Rate  Equations  Solution 


Case 

■^e 

"K 

1/cm^ 

a (ERE) 
1/cml/sec 

l/cm5/sec 

A 

4*’ 

1.0 

1.6315 

5.9-12 

7.0-12 

B 

1.0^ 

4.4I* 

3.04-12 

3.4-12 

C 

1.0^ 

I.5I* 

2.16-12 

2.25-12 

D 

1.0^ 

6.59II 

1.62-12 

1.8-12 

E 

l.o'* 

5.0I2 

8.04-15 

9.5-15 

P 

1.0^ 

1.5!^ 

6.77-15 

8.1-15 

G 

1.0* 

5. oil 

4.58-15 

6.5-15 

H 

1.4* 

I.63I5 

1.91-12 

2.1-12 

I 

1.4* 

I.5I* 

7.98-15 

1.05-12 

J 

1.4* 

5.0l=* 

3.69-15 

5.6-15 

K 

1.4* 

5.0II 

2.48-15 

3.7-15 

L 

1.0* 

5.381* 

3.32-12 

3.7-12 

ri 

1.0* 

2.69II 

1.26-12 

1.35-12 

N 

1.0* 

2.69!^ 

6.86-15 

8.5-15 

0 

1.0* 

2.69II 

3.97-15 

6.0-15 

p 

8.0^ 

I.5I* 

4.43-12 

4.1-52 

Q 

6.o3 

I.5I* 

1.25-11 

8.5-12 

®Chen  (Ref.  21) 

^^Superscripts  denote  powers  of  10  by  which  the 
numbers  are  to  be  multiplied. 
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the  simplified  technique  of  obtaining  the  minimum  in  the  total  rate  of 
de-excitation  as  a  function  of  the  energy  levels  of  excited  states  (Refs. 

2  and  23).  These  differences  can  easily  account  for  the  slight  differ¬ 
ences  in  Table  3.  It  is  believed  that  the  inconsistent  behavior  of  case 
Q  is  due  to  explicit  inclusion  of  exchange  excitations  from  the  metas¬ 
table  2^S  level  to  the  adjacent  metastable  2^S  level.  This  has  the  effect 
of  maintaining  a  higher  collisional  rate  to  the  adjacent  2^P  level  (rapid 
direct  excitation)  so  that  a  higher  rate  of  radiation  to  the  ground  state  is 
maintained.  This  rate  at  QSS  is  communicated  back  through  the  bound 
states  to  the  free  electrons,  thus  resulting  in  a  high  rate  of  recombina¬ 
tion.  Other  workers  (Ref.  49)  have  also  calculated  recombination  co¬ 
efficients  but  their  conditions  and  resolution  of  their  data  make  compar¬ 
isons  difficult. 

It  is  notable  that  these  recombination  coefficients  are  typically  a 
factor  of  two  less  than  those  reported  for  hydrogen  by  Bates,  et  al. 

(Ref.  3).  This  corresponds  to  a  general  trend  in  reported  experimental 
values  for  helium  recombination  rates  to  be  less  than  hydrogen  recom¬ 
bination  rates  (Refs.  12,  22,  50,  and  51).  The  true  value  will  be  some¬ 
what  in  doubt  because  of  the  inability  of  the  Gryzinski  cross  sections  to 
predict  accurately  the  true  collision  rates.  However,  the  use  of  a 
consistent  set  of  cross  sections  for  both  hydrogen  and  helium  will 
indicate  relative  differences. 

The  favorable  comparison  of  the  results  of  the  transient  ERE 
solution  and  the  QSS  calculations  shown  in  Table  3  indicates  that  the 
numerical  technique  used  for  the  ERE  will  yield  physically  correct 
results.  This  observation  coupled  with  the  stability  considerations 
described  in  Section  4.  0  constitutes  proof  that  the  numerical  results  to 
be  presented  in  the  following  pages  do  represent  physical  processes 
to  the  limit  of  the  accuracy  of  the  parameters  used  in  the  calculation, 

5.2.2  Location  of  the  Critical  Level 

The  initial  distribution  for  each  of  the  plasmas  was  a  Boltzmann 
distribution  in  equilibrium  with  the  free  electrons  from  the  2^S  level 
through  the  upper  states  considered.  As  the  ERE  solution  develops  in 
time,  the  levels  below  the  critical  level  adjust  to  the  final  QSS  distri¬ 
bution  while  those  levels  above  and  including  the  critical  level  remain 
in  Saha  equilibrium  with  the  free  electrons. 

Although  the  concept  of  a  critical  level  is  not  used  in  the  present 
computations,  the  observation  of  the  development  of  such  a  level  in  the 
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transient  solution  and  the  location  of  it  is  of  interest  as  it  pertains  to  the 
internal  consistency  of  the  calculations.  The  location  of  the  critical 
level  is  useful  for  experimental  determinations  becuase  it  indicates 
that  lowest  level  whereby  spectral  intensity  measurements  can  be  used 
to  ascertain  the  kinetic  temperature  of  the  free  electron  distribution. 

Table  4  lists  the  critical  levels  found  for  each  of  the  plasmas.  Two 
criteria  were  used  for  identification  of  the  critical  level:  3  percent  and 
10  percent,  referring  to  the  maximum  deviation  from  Saha  equilibrium 
with  the  free  electrons  for  any  of  the  substates  in  a  principal  quantum 
level  for  that  principal  quantum  level  to  be  considered  as  the  critical 
level.  The  critical  levels  found  for  the  high  electron  density,  high 
temperature  cases  with  the  10 -percent  criteria  compare  favorably  with 
predictions  by  Hinnov  and  Hirschberg  {Ref.  11)  who  indicate  that  at  those 
conditions  quantum  number  4  will  be  the  lowest  level  in  Saha  equilibrium 
with  the  free  electrons.  No  comparisons  were  made  at  the  low  densi¬ 
ties  and  temperature  other  than  that  they  qualitatively  show  expected 
behavior. 


Table  4.  Total  Number  Density,  Electron  Density,  Electron 
Temperature,  and  Location  of  Critical  Level  for 
Each  of  the  Plasmas  Studied 


Case 

"o  , 

"e 

Critical  Level 

^  3 
1/cm 

l/cm^ 

“K 

(Principal 

3% 

Quantum  Numbei 
10% 

A 

1.63I’* 

1.63^^ 

1.0^ 

4 

4 

B 

4.41® 

4-.414 

1.0^ 

5 

4 

C 

1.5I& 

I.5I4 

1.0^ 

5 

4 

D 

6.59^5 

6.59^^ 

1.0« 

6 

5 

E 

5.0I2 

1.0^ 

8 

6 

F 

I.5I4 

1.5^2 

1.0^ 

9 

7 

G 

5.0^^ 

5.0^^ 

l.o'‘ 

10 

7 

H 

1.63^ 

1.63^^ 

1.4^ 

4 

3 

I 

l.S^® 

1.5“ 

1.4* 

5 

4 

J 

5.0^4 

5.0^2 

1.4* 

B 

6 

K 

5.0“ 

5.0^^ 

1.4* 

9 

7 

L 

2.69^^ 

5.38^^ 

1.0* 

5 

4 

M 

2.69^® 

2.69^2 

1.0* 

6 

5 

N 

2.69^^ 

2.69^2 

1.0* 

8 

6 

0 

2.69^ 

2.69^^ 

1.0* 

10 

8 

P 

i.sie 

1.5!^- 

B.O^ 

5 

4 

Q 

i.s^fi 

i.si^ 

6.0^ 

6 

5 

^Superscripts  denote  powers  of  10  by  which  the 
numbers  are  to  be  multiplied. 
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5.3  TIME  DEVELOPMENT  OF  THE  ERE  SOLUTION 


The  time  development  of  selected  quantum  level  densities  for  eight 
of  the  cases  listed  in  Table  2  are  presented  on  the  following  pages. 

The  ordinate  for  each  plot  is  the  ratio  of  the  instantaneous  population 
density  to  the  density  required  to  be  in  Saha  equilibrium  with  the  free 
electron  density.  Thus,  at  t  =  0.0  sec,  each  quantum  level  would  have 
an  ordinate  of  1.0.  The  substates  of  only  quantum  levels  2  and  3  are 
shown.  Those  quantum  levels  between  quantum  level  3  and  the  critical 
level  will  show  characteristics  similar  to  the  triplet  and  singlet  states 
of  quantum  level  3,  and  any  quantum  level  above  and  including  the 
critical  level  will  show  a  constant  ratio  of  1.0. 

The  two  quantum  levels  are  displayed  together  because  substates 
of  quantum  level  3  are  characteristic  of  the  other  upper  levels  while 
substates  of  quantum  level  2  are  in  general  quite  different  from  each 
other  and  any  of  the  other  levels.  The  ground  state  density  is  not 
included  since  it  shows  nothing  more  than  a  slow,  steady  increase 
throughout  time.  The  study  of  the  rate  of  change  of  the  ground  state 
density  is  more  pertinent  and  studies  of  certain  of  these  will  be 
presented  later. 


Figure  8  shows  the  time  development  for  Case  A,  a  high  density 
plasma  which  is  generally  collision  dominated.  Figure  8  does  not 
include  any  of  the  substates  of  quantum  level  3  since  their  values  fall 
into  the  same  region  as  the  2^S  and  2^P  levels  and  would  just  clutter 
the  plot. 


This  plasma  is  sufficiently  collision  dominated  that  the  only  signifi¬ 
cant  change  is  due  to  the  radiative  depopulation  of  the  2lp  level.  The 
metastable  2^S  and  the  2^P  levels  are  coupled  together  strongly  by  two- 
body  collisional  processes  so  that  as  the  2^P  level  decays  radiatively, 
it  is  depressed  below  the  equilibriimi  configuration.  Consequently,  the 
collisional  excitation  rate  from  the  2^S  level  is  greater  than  the  super¬ 
elastic  de- excitation  from  the  2^P  to  the  2^S  level,  and  there  is  a  net 
transfer  of  atoms  from  the  lower  to  the  upper  state.  The  2^P  level  con¬ 
tinues  to  decay  radiatively  and  thus  causes  the  2^S  level  to  decay  also. 

At  about  1  X  10"®  sec,  the  2^P  density  has  depressed  sufficiently  that 
the  radiative  depopulation  rate  has  decreased  to  the  point  that  the 
collisional  population  rate  will- balance  the  total  rate  of  change  towards 
zero,  or  the  QSS.  There  are  other  collisional  processes  possible, 
namely  exchange  collisions,  so  that  as  the  2^P  density  decreases,  the 
other  quantum  levels  in  the  vicinty  of  the  2^P  level  are  also  depressed  by 
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Tlae  (see) 

Figure  8.  n(p)/n(Saha)  versus  time  for  quantum  level  two,  case  A. 


collisional  transfers  to  the  2^P  and  2^S  levels,  and  they  all  show  a 
slight  depression  below  the  equilibriiim  value.  All  of  the  singlet  states 
of  quantum  level  3  show  a  tendenQr  to  be  populated  slightly  less  than 
the  triplet  states,  reflecting  the  more  efficient  collision  processes  to 
the  2lP  state  from  the  other  singlet  states.  However,  the  difference 
between  the  singlet  and  triplet  states  is  only  about  10  percent,  and 
graphical  distinction  is  not  justified. 

Figure  9  shows  the  time  development  of  the  substate  of  quantum 
levels  2  and  3  for  case  C.  The  3®P  and  3^D  as  well  as  the  S^S,  3^D, 
and  3^P  levels  are  close  enough  together  energetically  that  at  this 
density  they  are  coupled  together  strongly  and  decay  together.  Thus, 
the  two  trivet  states  are  plotted  together  as  a  single  curve  as  are  the 
three  singlet  states  of  quantum  level  3.  Case  C  is  at  a  lower  total 
density  than  case  A  and  Fig.  9  illustrates  this  through  the  larger  spread 
in  the  various  densities.  Again  the  2^P  level  tends  to  depopvilate  radi- 
atively  from  t  =  0.  0  sec.  The  depression  of  the  2^P  density  tends  to 
depress  the  density  of  the  other  singlet  states  adjacent  to  it  because  of 
collisional  transfer  to  the  2^P  state.  The  largest  effect  is  on  the  2ls 
level  because  the  2^S  and  2^P  levels  are  closer  together  energetically 
than  the  2lp  and  3^S  levels. 
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Figure  9.  n(p)/n(Saha)  versus  time  for  quantum  levels  two  and  three,  case  C. 


The  metastable  2^8  level  has  a  large  radiative  population  rate  so 
that  from  t  =  0.  0,  it  is  increasing  in  density.  Even  though  the  2^8  level 
is  collision  dominated,  the  radiative  rate  is  of  comparable  magnitude 
at  this  density.  Since  the  initial  condition  for  this  plasma  was  a  Boltz¬ 
mann  distribution,  the  coUisional  population  and  depopulation  rates 
balance  each  other  so  the  net  effect  on  the  total  rate  for  the  23s  level 
is  radiative.  At  about  3  x  10"  sec,  the  23s  and  2^P  densities  {the 
2^P  increases  because  of  strong  coUisional  coupling  to  the  2^S  level) 
have  increased  in  density  enough  that  exchange  transitions  to  the  2^8 
and  2^P  levels  balance  the  radiative  rate  from  the  2lP  level  and  then, 
as  the  triplet  states  continue  to  increase  in  density,  the  coUisional 
rate  dominates  and  the  densities  of  the  singlet  states  increase  until 
the  plasma  achieves  QSS  and  computations  are  terminated. 


Figure  10  shows  the  time  development  of  the  2^P  radiative  rate, 
two-body  coUisional  rate,  and  total  rate  for  case  C.  Each  quantity  is 
the  algebraic  sum  of  the  populating  and  depopulating  processes.  As 
can  be  seen,  the  major  contributor  to  the  total  rate  for  short  times  is 
radiative  depopulation. 
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Figure  10.  2^  P  collisional,  radiative,  and  total  rates 
versus  time,  case  C. 


Throughout  this  early  time  period,  the  radiative  rate  dominates 
sufficiently  that  the  2^P  density  decreases.  Thus,  the  radiative  rate 
also  decreases.  This  decrease  in  density  causes  a  net  depopulation  of 
adjacent  states  so  the  coUisional  rate  decreases.  In  the  region  1  x  10~° 
sec  5  t  <  3  X  10"®  sec  the  2^S  and  2®P  levels  attain  a  high  enough 
density  that  the  coUisional  rate  begins  to  increase.  Although  the  scale 
of  Fig.  10  is  too  large  to  provide  visual  resolution,  at  about  3  x  10"® 
sec  the  total  rate  becomes  slightly  positive,  reflecting  the  balance  be¬ 
tween  the  radiative  and  coUisional  rates.  Subsequently,  the  coUisional 
rate  continues  to  increase  because  of  the  continuing  increase  in  density 
of  the  adjacent  levels.  Reflecting  the  increase  in  the  2lp  density,  the 
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radiative  depopulation  rate  also  tends  to  increase.  Finally,  at  about 
1  X  10“6  sec  the  coUisional  rate  begins  to  flatten  because  QSS  is  being 
attained  by  the  other  levels  and  the  collisional-radiative  processes 
strike  a  final  balance  for  the  2lp  level. 

Figure  11  shows  the  two-body  coUisional  depopulation  rate,  the 
two-body  coUisional  population  rate,  and  the  radiative  rate  for  the  2^8 
level  of  case  C.  The  radiative  rate  has  been  increased  an  order  of 
magnitude  for  plotting  purposes,  and  the  coUisional  depopulation  is 
plotted  as  a  positive  quantity. 


Figure  11.  2'  S  radiative  population,  two-body  coUisional 

population  and  two-body  coUisional  depopulation 
rates  versus  time,  case  C. 
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The  decay  of  the  2^S  level  is  dominated  strongly  by  the  collisional 
depopulation  rate  until  about  about  3  x  10“8  gee.  The  collisional  popu¬ 
lation  rate  shows  a  marked  decrease  through  the  early  portion  of  the 
decay,  primarily  because  of  the  decrease  in  the  2lp  density.  By  about 
2  X  10 "8  sec,  however,  the  2^'P  density  has  decreased  enough,  and  the 
23s  and  2^P  densities  increased  enough,  that  the  exchange  collisions 
from  the  23s  and  23p  levels  are  an  appreciable  factor  in  the  2^S  colli¬ 
sional  population  rate.  Subsequently,  as  these  triplet  states  increase, 
so  does  the  collisional  population  rate  of  the  2^S  level.  The  collisional 


depopulation  rate  shows  a  very  steep  decrease  until  about  3  x  10"®  sec 
reflecting  the  rapid  decrease  of  the  2lS  level  because  of  this  process. 
At  3  X  10 “3  sec,  however,  the  radiative  rate  and  the  collisional  popula¬ 
tion  rate  are  sufficient  to  dominate  the  total  rate,  so  the  2^S  density 
increases.  With  the  increase  in  density,  the  collisional  depopulation 
rate  also  increases.  The  noticeable  difference  between  the  collisional 


population  and  depopulation  rates  reflects  the  fact  that  the  2^P  density 
is  below  the  equilibrium  configuration  with  the  2  density,  and  colli¬ 
sional  processes  are  attempting  to  bring  the  2lp  density  into  this  con¬ 


figuration.  The  flattening  of  the  rates  at  about  10-6  gee  shows  that  QSS 
for  the  entire  distribution  is  being  achieved.  At  10"®  sec,  the  2ls  total 
rate  is  6  x  lO^  atoms/cm^/sec,  significantly  less  than  any  of  those 
rates  plotted,  and  this  shows  the  collisional  radiative  balance  which  is 


maintained  at  QSS. 


Figure  12  shows  the  time  development  of  the  substates  of  quantum 
levels  2  and  3  for  case  E.  This  plasma  is  at  the  same  electron  tempera¬ 
ture  and  degree  of  ionization  as  cases  A  and  C  but  the  total  density  is 
lower.  Hence,  the  plasma  is  more  radiation  dominated  and  the  spread 
of  the  various  quantum  levels  is  more  dramatic.  Because  the  plasma 
is  less  collision  dominated,  the  2^S  density  tends  not  to  decrease  as 
far  as  in  the  previous  cases.  The  triplet  states  of  quantum  level  3 
show  some  separation  at  the  early  and  the  later  times  and  the  2^1*  level 
shows  a  very  dramatic  decrease  in  density. 

The  separation  of  the  singlet  and  triplet  states  of  quantum  level  3 
is  much  more  evident  here.  The  P  and  D  states  are  quite  close  together 
energetically  and  thus  are  cblllsionally  coupled  together  quite  strongly. 
The  3 1 P  state  is  influenced  by  the  large  radiative  rate  to  the  ground 
state  and  hence  the  3^P  and  3^D  levels  show  a  noticeable  reduction  in 
density  compared  to  the  3^8  state.  In  the  early  stages  of  decay  the  3^? 
rate  is  determined  principally  the  radiation  to  the  23s  level  so  that 
the  36 P  and  3^0  levels  decay  rapidly.  The  3^8  level  rate  is  also 
largely  determined  by  radiation  to  the  2  P  level  but  this  proceeds  at  a 
slower  rate  than  the  33p,  3®D  decay.  At  about  2  x  lO""^  sec,  the 
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Tlae  (sec) 


Figure  12.  n(p)/n(Saha}  versus  time  for  quantum  levels 
two  and  three,  case  E. 

density  of  the  23s  level  has  increased  enough  that  collisional  transfers 
to  the  2^F  cause  its  density  to  begin  to  increase.  During  the  period 
of  balancing  of  collisional-radiative  rates  between  the  triplet  states 
of  quantum  levels  2  and  3,  6  x  10"®  sec  £  t  <  3  x  10"^  sec,  collisional 
processes  between  the  triplet  states  of  quantiim  level  3  can  bring  these 
states  into  a  Boltzmann  configuration  with  each  other.  After  about 
3  X  10"*^  sec,  the  2^8  and  2^P  densities  have  increased  to  the  point  that 
the  triplet  states  also  start  to  increase  in  density.  Because  of  the 
larger  radiative  rate  from  the  3^ P  level,  and  the  fact  that  as  it  adjusts 
it  must  also  maintain  the  state  in  a  Boltzmann  configuration,  the 
33p  and  33d  states  do  not  show  as  marked  an  increase  in  density  as 
the  3^8  state. 


S8 


AEDC-TR-76-5 


Figure  13  shows  the  same  quantum  levels  for  case  J  at  the  same 
conditions  as  case  E  except  the  electron  temperature  is  14.000°K.  As  a 
result  the  collisional  processes  for  this  plasma  are  stronger  than  case 
E  and  the  curves  do  not  show  qmte  as  much  spread.  The  3^8  density  is 
separated  somewhat  more  from  the  3^?  and  3^D  densities  because  of 
the  stronger  coupling  between  the  2^?  and  3^8  state  at  the  higher  temper¬ 
ature.  The  2^8  level  tends  to  rise  to  a  higher  density  for  case  J  than 
for  case  E.  This  is  because  collisional  processes  are  more  efficient  at 
the  higher  temperature  so  that  the  collisional  population  and  depopula¬ 
tion  rates  remain  closer  together,  tending  to  cancel  each  other.  Thus, 
the  2^8  level  must  rise  to  a  higher  density  in  order  for  collisional 
depopulation  to  dominate  collisional  population  sufficiently  to  also 
cancel  the  radiative  population  rate. 


Figure  13.  n(p)/n(Saha)  versus  time  for  quantum  levels 
two  and  three,  case  J. 
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Figure  14  shows  the  time  development  of  substates  of  quantum  levels 
2  and  3  for  case  N.  This  plasma  is  at  a  different  total  density  than  case 
E  (5.  0  X  atoms/cm^  for  case  E  and  2.  69  x  10^^  atoms/cm^  for 


Figure  14.  n(p)/n(Saha)  versus  time  for  quantum  levels  two  and  three,  case  N. 

case  N),  but  the  electron  densities  are  nearly  the  same  (5  x  10^^  for 
case  E  compared  with  2.  69  x  10^^  for  case  N),  and  the  electron  tempera¬ 
tures  are  the  same.  The  character  of  the  curves  is  about  the  same  in 
Figs.  14  and  12,  with  the  slight  differences  reflecting  the  slight 
decrease  in  the  magnitude  of  collisional  processes  due  to  the  lower 
electron  density.  This  is  well  illustrated  by  the  fact  that  the  2^8  density 
tends  to  be  slightly  greater  than  the  2^ P  level  until  about  3  x  sec. 
This  is  due  to  the  small  decrease  in  the  collisional  coupling  of  the 
2 Is  to  the  2^P  level  while  the  radiative  transitions  from  the  to  the 
2^8  level  play  a  proportionately  greater  role  for  the  early  portion  of  the 
decay.  The  2^P  and  2^8  levels  are  not  included  in  their  entirety  in 
Fig.  14  since  their  characteristics  are  now  well  documented  by  the 
previous  figures.  The  2lp  level  reaches  a  minimum  ordinate  of  0.0023 
at  about  5  x  10 sec  and  returns  to  about  0.  01  at  QSS.  The  2^8  level 
has  a  peak  ordinate  of  96.2  at  QSS,  t«4  x  10“^  sec. 

Figures  15,  16,  and  17  show  the  time  development  of  the  substates  of 
quantum  levels  2  and  3  for  cases  I,  Q,  and  P,  respectively.  These 
plasmas  are  all  at  the  same  densities  and  degrees  of  ionization  but  have 
electron  temperatures  of  14, 000°K,  8,000“K,  and  6,000“K,  respectively. 
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Figure  15.  n{p)/n(Saha)  versus  time  for  quantum  levels  two  and  three,  case  I. 


n(p) 

n(Saha) 


Figure  16.  n(p)/n(Saha)  versus  time  for  quantum  levels  two  and  three,  case  P. 
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Figure  17.  n(p)/n(Saha)  versus  time  for  quantum  levels  two  and  three,  case  Q. 


Figure  9  (case  C)  showed  the  development  of  a  plasma  at  the  same  den¬ 
sity  and  degree  of  ionization  with  an  electron  temperature  of  10, 000“K. 
Comparison  of  Figs.  15,  9,  16,  and  17  shows  the  steady  decrease  of 
importance  of  collisional  processes  and  the  steady  increase  of  radia¬ 
tive  processes  as  temperatvire  is  decreased. 

The  interplay  of  these  collisional -radiative  processes  is  quite  com¬ 
plex  and  subtle.  Detailed  study  of  the  change  in  the  decay  characteris¬ 
tics  will  be  taken  up  later  in  this  study.  At  this  point,  it  will  siiffice  to 
say  that  at  the  low  temperature.  Fig.  17,  collisional  processes  are 
sufficiently  inefficient  that  the  radiative  transitions  of  the  2^P  level 
depress  it  far  enough  that  adjacent  levels  are  also  depressed  in  attempt¬ 
ing  to  bring  the  2^P  density  back  to  the  equilibrium  configuration.  This 
depression  of  adjacent  levels  in  turn  lowers  particularly  the  radiative 
transition  rate  into  the  2^S  level  so  that  it  does  not  show  a  marked  in¬ 
crease  in  density.  Without  this  increase  in  the  2^S  density,  there  is  no 
subsequent  increase  in  the  other  quantum  levels  because  of  collisional 
transfers  from  the  2^8  level.  Thus,  the  entire  distribution  in  the  low- 
lying  states  is  drawn  below  the  equilibrium  with  the  free  electrons. 
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5.3.1  Acquisition  of  Quasi-Steady  State 

A  consequence  of  the  QSS  for  the  optically  thin  plasma  is  the  condi¬ 
tion  that  the  rate  of  change  of  the  ground  state  density  and  the  continuum 
density  are  equal  in  magnitude  but  opposite  in  sign  {Section  2. 0).  Hence, 
observation  of  the  time  development  of  the  ratio  rie/h(l>  will  provide 
graphic  illustration  of  the  time  development  of  the  plasnia  as  a  whole  to 
the  QSS.  At  QSS,  this  ratio  has  a  value  of  -1. 

The  initial  conditions  used  for  the  csdcvdations  will  certainly  have  a 
significant  bearing  upon  the  time  development  of  ng/n{l).  However,  as 
described  on  the  preceding  pages,  the  initial  conditions  were  suph  that 
the  metastable  2%  level  density  had  to  increase  to  reach  the  QSS  and 
the  other  levels  were  both  directly  and  indirectly  affected  the  metas¬ 
table  levels.  Thus,  if  a  physical  plasma  suffers  a  parametric  change 
such  that  the  densities  of  these  levels  are  below  their  QSS  value,  the 
initial  conditions  chosen  fpr  this  study  dp  give  realistic  resiilts. 

Figure  18  shows  the  time  development  of  ng/ri(l)  for  cases  A  through 
G  and  case  M.  These  plasmas  all  have  an  electron  temperature  of 
10, 000®K,  one  percent  ippization,  and  progressively  lower  total  densities. 
As  one  would  expect,  as  conditions  progress  such  that  the  plasma  is  less 
collision  dominated,  the  49  §99^41^9  becomes  progressively 


Figure  18.  he/n(1)  versus  time,  cases  A,  B,  C,  D,  E,  F,  G,  and  M. 


/ 
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longer.  These  times  range  from  3  x  10"®  sec  for  case  A  to  3.  5  x  10“^ 
sec  for  case  G.  The  criterion  used  for  these  and  subsequent  cases  is  that 
when  rig  and  h(l)  agree  within  5  percent  (-1.  05  <  hg/ri  (1)  <  -0,  95)  the 
plasma  is  considered  to  be  at  the  QSS.  A  consistent  (case  A  the  only 
exception)  characteristic  is  the  overshoot  of  the  QSS  condition  early  in 
the  decay  of  each  plasma  and  then  the  subsequent  slow  recovery  back  to 
the  QSS  conditions.  This  phenomenon  will  be  studied  in  detail  below. 

Figure  19  shows  the  time  development  of  he/n(l)  for  cases  H  through 
K.  These  plasmas  are  at  similar  densities  and  degrees  of  ionization  as 
cases  A,  C,  E,  and  G,  respectively,  except  that  the  electron  tempera¬ 
ture  is  elevated  to  14,000®K.  The  character  of  the  curves  in  Fig.  19 
compared  to  those  in  Fig.  18  is  quite  similar,  except  that  the  QSS  is 
achieved  slightly  sooner  (2.  5  x  10“®  sec  for  case  H  compared  to  2.  5  x 
10"®  sec  for  case  A,  2.  3  x  10"^  sec  for  case  K  compared  to  3.  5  x  10 
sec  for  case  G)  and  the  magnitude  of  overshoot  is  slightly  greater. 


Time  (sec ) 


Figure  19.  he/h(1)  versus  time,  cases  H,  I,  J,  and  K. 


Figure  20  shows  the  time  development  of  ne/n(l)  for  cases  L 
through  O.  These  cases  are  all  at  the  same  total  density  and  electron 
temperature  but  at  different  degrees  of  ionization.  The  character  of 
these  curves  is  also  consistent  with  those  in  Figs.  18  and  19. 

Figure  21  shows  the  time  development  of  ng/n(l)  for  cases  C,  I,  P, 
and  Q.  These  plasmas  are  all  at  the  same  density  and  ionization  but 
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Figure  20.  lie/nd)  versus  time,  cases  L,  M,  N,  and  O. 


10"®  10“®  lO"''  10"® 

Time  (sec ) 

Figure  21.  he/h(1)  versus  time,  cases  C,  I,  P,  and  Q. 
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each  is  at  a  different  electron  temperature.  This  again  shows  the  con¬ 
sistent  behavior  that  as  coUisional  rates  get  larger  (increasing  tempera¬ 
ture)  the  time  for  the  plasma  to  achieve  QSS  gets  progressively  smaller. 
One  can  also  see  the  tendency  for  the  overshoot  to  become  progressively 
larger  as  the  electron  temperature  becomes  larger.  Case  Q,  the  lowest 
temperature  case,  does  not  show  the  overshoot  and  will  be  discussed  in 
more  detail  below. 

Regardless  of  the  character  of  the  decay,  the  time  to  achieve  QSS 
was  not  extremely  sensitive  to  the  electron  temperature,  ranging  from 
2  X  lO""^  sec  for  case  I  (Tg  =  14, 000“K)  to  1.  2  x  10"®  sec  for  case  Q 
(T^  =  6,000»K). 

5.3.2  Overshoot  of  Quasi -Steady  State 

A  consistent  character  of  the  ng/n(l)  curves  is  the  overshoot  of 
the  QSS  early  in  the  decay  of  most  of  the  plasmas  (cases  A,  H,  and  Q 
the  only  exceptions)  and  then  the  slow  recovery  back  to  the  QSS.  This 
could  be  due  to  an  overshoot  of  QSS  conditions  of  either  ng  or  n(l)  but 
in  fact  both  quantities  show  the  character  to  a  greater  or  lesser  degree. 

It  is  also  tempting  to  ascribe  the  overshoot  to  numerical  inaccuracies 
in  the  program  but  it  is  in  fact  a  physical  phenomenon  depending  upon 
the  initial  conditions. 

The  precise  balancing  between  the  various  processes  included  in 
the  ERE.  induces  large  time  gradients  in  the  excited  state  densities 
during  the  early  stages  of  computation.  However,  these  numerical 
requirements  are  satisfied  very  early  (t  <  10"^®  sec)  in  the  decay  with 
no  appreciable  changes  in  the  excited-state  densities.  The  subsequent 
plasma  decay  is  attributable  to  the  physical  imbalance  of  the  terms  in 
the  ERE  because  of  the  initial  distribution.  Thus  physical  plasmas 
which  suffer  parametric  changes  which  require  a  significant  2^S  density 
increase  to  achieve  QSS  will  also  be  subject  to  the  overshoot  phenomenon. 

Case  C  was  chosen  to  examine  the  overshoot  phenomenon  in  more 
detail.  This  case  shows  the  phenomenon  quite  well  yet  does  not  require 
such  large  times  to  re-establish  QSS  that  study  can  become  imwieldy. 

Figure  22  shows  the  detailed  time  development  of  the  rates  of 
change  of  the  continuum  as  well  as  the  l^S,  2^S,  and  2^?  bound  states 
for  case  C.  This  shows  that  both  the  ground  state  and  free  electron 
rates  overshoot  their  QSS  values  and  converge  back  toward  their  QSS 
value  slowly.  Of  specific  interest  in  this  context  is  the  2^P  level.  By 
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Time  (sec ) 


Figure  22.  Total  rate  of  change  of  VS,  2^8,  2^P  quantum  levels  and  free 
electron  density  versus  time,  case  C. 


comparing  the  shape  of  the  curves  for  the  time  development  of  the  2^P 
and  continuum  rates,  it  is  seen  that  the  continuum  overshoot  follows 
the  sudden  growth  of  the  2^P  rate  from  about  5  x  10"^  sec  through  about 
3  X  lO”"^  sec  and  then  slowly  falls  to  the  QSS  value  with  the  fall  of  the 
2^S  and  23p  rates  to  the  QSS.  In  this  region  of  time  some  of  the  other 
quantum  levels  will  also  show  a  similar  variation  in  their  rates  for  the 
described  below.  The  early  part  of  the  continuum  rate  curve, 
sec,  follows  the  adjustment  of  the  other  quantum  levels  toward 


reason^ 
t  <  lO’” 
their  QSS  value. 


Figure  23  examines  the  time,  developnient  of  the  2^P  rate  in  more 
detail  illustrating  both  the  collisional  and  radiative  contributions  to  the 
net  rate.  Initial  conditions  for  this  plasma  were  such  that  the  23p 
density  is  less  than  its  QSS  value  so  there  is  a  tendency  for  the  2^P 
level  to  populate  radiatively  in  the  early  portion  of  the  decay.  As  the 
level  continues  to  populate  radiatively,  the  density  increases  so  that 
radiative  population  consequently  increases  and  collisional  depopulation 
decreases. 


Hence,  as  time  progresses,  the  radiative  rate  becomes  progressive¬ 
ly  smaller  and  the  collisional  depopulation  rate  becomes  larger  for 
t  <  5  X  10“^  sec.  At  this  point,  however,  the  collisional  rate  reverses 
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Figure  23.  2^P  collisional  rate,  radiative  rate,  and 
total  rate  versus  time,  case  C. 

directions  from  a  depopulating  effect  to  a  populating  effect.  With  the 
subsequent  increase  in  the  population  density,  the  radiative  rate  serves 
to  increasingly  depopulate  the  state  until  about  10"^  sec  when  the  colli¬ 
sional  and  radiative  rates  balance  each  other  and  the  state  approaches 
QSS.  By  referring  to  Fig.  22  one  sees  that  throughout  this  oscillation  of 
the  23p  rate  the  metastable  2^S  level  has  been  populating  radiatively  at 
a  high  rate.  At  about  5  x  10“®  sec  the  2^S  level  has  become  sufficiently 
populous  that  direct  excitation  2^8  -  2^P  collisions  with  the  continuum 
can  begin  to  have  a  discernible  effect  upon  the  23p  rate.  Collisional 
transfers  out  of  the  23$  level  to  all  levels,  however,  are  not  sufficient 
to  balance  the  radiative  rate  into  it  so  that  the  23s  level  continues  to 


68 


AEDC-TR-76-5 


populate,  and  thus  the  transfer  rate  to  the  2^?  level  becomes  larger 
until  QSS  is  reached.  As  the  collisional  rate  from  the  2^S  to  the  2^? 
rate  becomes  large  enough  to  begin  to  balance  the  radiative  rate,  the 
2^S  level  rate  decreases  to  its  QSS  value.  As  its  total  rate  decreases 
toward  zero,  the  collisional  rate  to  the  other  quantum  level  flattens 
out  so  that  the  other  collisional  radiative  processes  can  balance  the 
total  rate  for  each  individual  quantum  level  to  the  QSS  value.  The  slow 
equilibration  of  the  collisional  and  radiative  rates  for  the  2^S  and  2^P 
levels  is  communicated  by  collisional  processes  back  through  the  bound 
states  to  the  electron  continuum.  Thus  as  these  rates  fall  to  their  QSS 
value,  they  have  a  decreasingly  important  contribution  to  the  continuum 
rate,  and  the  continuum  rate  can  also  slowly  subside  to  the  QSS  value. 


Figure  24  illustrates  the  time  development  of  the  population  densi¬ 
ties  of  the  2^S,  2^P,  2^S,  and  2^P  levels  for  case  C.  This  shows  the 
slow  increase  of  the  triplet  states  to  their  QSS  value  while  examination 
of  the  singlet  states  reveals  the  cause  of  the  ground  state  rate  overshoot 
shown  in  Fig,  22.  The  2^P  level  has  a  large  radiative  transition  proba¬ 
bility  to  the  ground  state.  Hence,  its  population  density  shows  a  very 
rapid  decrease  through  the  early  part  of  the  decay.  As  its  density  falls 
below  the  equilibrium  configuration,  collisions  have  a  tendency  to  fill  it 
and  especially  collisional  transfers  from  the  2^S  state.  Hence,  the  2^S 
state  also  shows  a  reduction  in  density  but  lags  slightly  behind  the  2^P 
level.  At  about  1  x  10"®  sec,  the  2^P  level  has  lowered  in  density 


Figure  24.  2^8,  2^P,  218,  and  2^P  population  densities  versus  time,  case  C. 
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sufficiently  that  coUisional  rates  can  begin  to  balance  radiative  rates, 
so  the  density  curve  begins  to  flatten.  However,  the  triplet  states  are 
continuing  to  increase  in  density  so  that  exchange  transfers  between  the 
triplet  and  singlet  states  become  more  numerous  and  serve  to  elevate 
the  population  densities  of  the  singlet  states.  The  largest  single  contri¬ 
butor  to  the  ground  state  rate  is  radiative  transitions  from  the  2  ^  P 
state.  By  comparing  the  ground  state  rate  in  Fig,  22  to  the  2^P  popula¬ 
tion  density  in  Fig.  24,  one  sees  that  as  the  density  decreases,  goes 
through  a  minimum,  and  slightly  increases,  so  does  the  ground  state 
rate.  This  effect  of  the  2^?  level  upon  the  ground  state  rate  is  preva¬ 
lent  for  all  cases  which  show  overshoot. 


In  summary,  those  cases  that  show  an  overshoot  of  ng/n(l)  do  so 
basically  because  of  the  need  of  the  2^S  density  to  increase  to  its  QSS 
value.  This  increase  in  density  causes  the  continuum  rate  overshoot 
because  the  coUisional  transfers  between  the  2^S  level  and  adjacent 
quantum  levels  induce  time  rates  of  changes  in  these  levels  which  are 
significant  with  respect  to  the  continuum  and  ground  state  rates.  The 
ground  state  rate  overshoot  is  caused  specifically  because  the  increase 
in  the  2^S  density  induces  an  increase  in  the  2^P  density  with  an  atten¬ 
dant  increase  in  the  1^S-2^P  radiation. 


5.3.3  No  Quasi-Steady-State  Overshoot 

As  indicated  previously,  three  of  the  plasmas  studied  showed  no 
overshoot  of  the  QSS  condition  but  rather  decayed  directly  to  QSS.  Two 
of  these  cases,  A  and  H,  are  high  density  cases  which  are  strongly 
coUision  dominated.  Little  needs  to  be  said  for  these  two  as  to  why  they 
do  not  show  the  overshoot  characteristics.  Since  they  are  strongly  coUi¬ 
sion  dominated,  slight  changes  in  any  quantum  level's  density  is  quickly 
communicated  to  the  other  levels  so  that  the  QSS  distribution  is  quickly 
reached.  No  single  process,  e.  g. ,  2^P  radiation,  can  cause  one  quan¬ 
tum  level  to  over -depopulate  but  rather  all  levels  approach  QSS  together. 

The  other  case  which  shows  no  overshoot,  case  Q,  is  somewhat 
different  in  nature,  however.  In  this  low  temperature  situation  (Tg  = 
6,000“K)  the  coUisional  processes  are  sufficiently  inefficient  that  the 
radiative  decay  of  the  2^P  density  acts  as  a  sink  for  the  rest  of  the 
distribution. 

Table  5  shows  the  radiative  popvilation  rate  for  the  2^S  level  for  each 
of  the  cases  C,  I,  P,  and  Q  at  initial  conditions  and  at  QSS.  Inspection 
of  this  table  shows  that  of  these  cases,  case  Q  is  the  only  one  which 
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Table  5.  Comparison  of  23$  Radiative  Rates  at  Initial  Conditions 
and  Quasi-Steady  State  for  Casas,  C.  I,  P,  and  Q 


Case 

T.. 

2^S  Radiative  Rate 

e 

atoms/cm^/seo 

Initial 

QSS 

I 

1.4’ 

7.98^® 

1.19^® 

C 

1.0^ 

2.91^* 

4.11^® 

P 

8.0^ 

8.97^® 

I.O4I7 

Q 

6.0^ 

6.14^^ 

3.77I’ 

^Superscripts  denote  powers  of  10  by  which  the 
numbers  are  to  be  multiplied. 


shows  a  decrease  in  the  rate  with  which  the  2^8  level  increases  because 
of  radiative  transitions.  This  implies  that  particularly  the  2^P  level 
has  suffered  a  significant  decrease  in  density  (at  QSS  for  case  I,  the 
23p  contributes  55  percent  of  the  2^8  radiative  rate).  To  understand, 
how  this  comes  about,  it  is  informative  to  examine  the  rate  coefficients 
for  collisional  transitions  between  the  substates  of  quantum  level  2. 
These  are  shown  in  Table  6  for  electron  temperatures  of  6,  000®K  and 
10,  000“K  for  comparison.  Included  in  the  table  for  discussion  purposes 
is  the  minimum  temperature  the  projectile  electron  must  exhibit  to 
effect  the  transition. 


Table  6.  Comparison  of  Collisional  Rate  Coefficients  for 
Quantum  Level  Two  Transitions  at  10,000”  K 
and  6,000”  K 


Transition 

Rate  Coefficients  omVsec 

P  =  1.0®*K  =  6.0^"K 

e  e 

Minimum 

Temperature 

•K 

2^S  -  2^S 

-9® 

7.03  ^ 

3.10"® 

9.2® 

2^S  -  2^P 

1.58“^ 

3.98"® 

1.3" 

2^S  -  2^P 

3.76"® 

1.17"® 

1.6® 

2^S  -  2^P 

1.42-® 

8.98"® 

4.0® 

2^S  -  2S 

1.21-® 

4.79"’ 

6.9® 

2^P  -  2^P 

1.77"® 

1.24"® 

3.0® 

^Superscripts  denote  powers  of  10  by  which  the 
numbers  are  to  be  multiplied. 
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Examination  of  the  rate  coefficients  at  10,000^K  shows  that,  as 
would  be  expected,  the  direct  transitions,  23S  -  2^P  and  2^S  -  2^P, 
are  the  most  important  by  at  least  an  order  of  magnitude.  However,  at 
6,000°K  the  2^F  -  2^P  transition  is  competitive  with  the  2^S  -  2^P 
transition.  This  is  because  at  6, 000"K  most  of  the  electron  distribution 
is  still  available  to  effect  the  transition  (a  3,000°K  electron  has  suffi¬ 
cient  energy)  and  even  though  it  is  an  exchange  transition,  the  magnitude 
of  the  rate  coefficient  is  not  significantly  altered.  For  the  direct  2^8  - 
2^P  transition,  a  13, 000®K  electron  is  required  and  at  6,  000“K  an  appreci¬ 
able  portion  of  the  electron  distribution  has  insufficient  energy  to  effect 
the  transition.  Hence,  the  rate  coefficient  is  markedly  reduced.  Simi¬ 
lar  arguments  hold  for  the  2ls  -  23p  transition,  which  requires  a 
4,000®K  electron  for  the  transition. 

Thus  at  6,000°K,  as  the  2 ^P  level  depresses  radiatively,  the  other 
levels  attempt  to  bring  it  back  to  the  equilibrium  configuration  and  are 
also  depressed.  The  details  of  this  depend  upon  the  following  coupling 
processes. 

The  2 Is  level  and  the  23p  level  are  the  major  contributors  to 
re  populating  the  2^P  level.  The  2^F  radiation' dominates  this  repopula¬ 
tion  so  that  these  two  contributors  are  also  depressed.  At  the  high 
temperatures,  the  2^S  -  2^P  rate  coefficient  dominates  the  2^P  -  2^P 
rate  coefficient  by  an  order  of  magnitude.  Thus,  the  2^P  level  is 
prevented  from  being  depopulated  by  23p  -  2^P  transitions.  At  the  low 
temperature,  the  2^8  -  23p  rate  coefficient  is  of  the  same  order  of 
magnitude  as  the  2^P  -  2^P  rate  coefficient  and  as  a  result,  the  2'^P  - 
2^P  transition  forces  the  23p  density  significantly  below  the  equilibrium 
configuration.  The  2^8  level  is  consequently  depleted  in  attempting  to 
fill  the  23 P  level.  This  depletion  of  the  2^8  density  in  turn  affects  all 
the  other  level  densities  and  thus,  at  6,000®K,  the  entire  distribution 
decays  below  equilibrium. 

The  foregoing  discussion  of  overshoot  has  shown  the  specific  mecha¬ 
nisms  by  which  overshoot  is  established  for  the  various  plasmas.  It 
has  been  shown  that  the  overshoot  characteristics  are  determined  princi¬ 
pally  by  the  2^S  state  and  the  2^P  state.  The  overshoot  of  the  quantity 
iie/n{l)  is  caused  by  the  overshoot  of  both  the  continuum  rate  and  the 
ground  state  rate  with  the  overshoot  of  the  ground  state  rate  showing  a 
somewhat  more  consistent  trend  with  varying  plasma  conditions.  This 
is  illustrated  in  Table  7,  which  has  reordered  the  various  plasmas 
studied  in  order  of  increasing  overshoot  of  n^/n  (1).  Included  in  the 
table  are  the  ng(max)/ne(QSS),  n(l)(min)/h(l)(QSS),  the  ratio  of  the 
2^8  collisional  depopulation  rate  to  the  radiative  population  rate  at 
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Table  7.  Maximum  Value  of  he/h(1),  he/n  (QSS),  Minimum  Value  of  n(1)/n(1)  (QSS), 
and  Ratio  of  Collisional  Rate  to  Radiative  Rate  at  Initial  Conditions  for  the 
2^S  and  2^  P  Quantum  Levels  for  Each  of  the  Afterglow  Plasmas  Studied 


Cue 

"a 

1/ca^ 

l/c»^ 

•K 

n^/nll> 

nasiauB 

overshoot 

n^(Bax) 

A^(oss) 

n(l)  (nln) 
A(l) (056) 

Initial  ' 

Collisional  Rate 
2^S 

Condxtione 

/  Radiative  Rate 
2^P 

A 

1.63“ 

1.0^ 

1.0 

1.0 

1.0 

17.12 

0.9 

H 

l.«3” 

1.63“ 

1.4« 

1.0 

1.0 

1.0 

15.04 

1.21 

L 

2.69« 

5.3*“ 

l.0« 

1.10 

1.06 

0.98 

5.67 

0.29 

■ 

4.4“ 

1.4“ 

i.o" 

1.1$ 

l.OB 

0.95 

4.64 

0.24 

P 

1.5“ 

1.5“ 

8.0^ 

1.30 

1.11 

0.84 

1.67 

0.06 

C 

1.5“ 

1.5“ 

l.0« 

1.50 

1.23 

0.82 

1.57 

0.08 

I 

1.5“ 

1.5“ 

1.4« 

1.65 

1.43 

0.87 

1.44 

0.10 

D 

6.59^® 

6.59“ 

i.o" 

1.75 

1.26 

0.73 

0.69 

0.036 

H 

2.69“ 

2.69“ 

1.0« 

1.95 

1.22 

0.63 

0.28 

0.015 

B 

5.0^* 

5.0“ 

1.0« 

2.10 

1.0» 

0.52 

0.05 

0.002 

P 

1.5“ 

1.5“ 

1.0« 

2.15 

1.06 

0.48 

0.016 

0.0008 

N 

2.69“ 

2.69“ 

1.0* 

2.15 

1.07 

0.50 

0.028 

0.0015 

G 

5.0“ 

s.o“ 

1.0* 

3.35 

1.04 

0.46 

0.005 

0.0002 

0 

2.69*® 

2.69*'' 

1.0* 

3.35 

1.04 

0.44 

0.003 

0.0002 

J 

5.0“ 

5.0“ 

1.4* 

3.45 

1.27 

0.53 

0.046 

0.0046 

X 

5.0“ 

5.0“ 

I.V 

3.26 

1.15 

0.35 

0.0046 

0.0013 

0 

1.5“ 

1.5“ 

6,0^ 

I.O 

1.0 

1.0 

1.64 

0.05 

^Superscripts  denote  powers  of  10  by  which  the  numbers 
are  to  be  multiplied. 

the  initial  conditions,  and  the  2^P  collisional  population  rate  to  the 
radiative  depopulation  rate  at  the  initial  conditions.  Case  Q  has  been 
included  at  the  bottom  of  the  table  for  completeness.  Table  7  shows 
that  for  those  cases  (A  and  H)  which  have  strong  collision  dominance 
in  the  2^8  level  and  competition  between  collisions  and  radiation  in 
the  2lp  level  there  is  no  overshoot  for  any  of  the  parameters.  In 
general,  as  radiation  becomes  more  important,  the  overshoot  of  lig/iid) 
becomes  greater.  The  ground  state  overshoot  tends  to  follow  a  more  or 
less  consistent  trend  with  changing  plasma  conditions  whereas  the  con¬ 
tinuum  rate  overshoot  appears  to  go  through  a  maximum  in  the  region  of 
conditions  in  which  there  is  no  dearcut  domination  of  the  2^8  collisional 
or  radiative  rates.  Cases  J  and  K  appear  to  be  the  principal  inconsistent 
cases  in  the  table  but  the  better  communication  between  the  lower  level 
rates  and  the  continuum  rates  have  already  been  discussed  for  these 
higher  temperature  cases.  It  should  also  be  noted  that  cases  J  and  K 
are  consistent  with  cases  H  and  I,  which  are  eilso  at  an  electron  tempera¬ 
ture  of  14,  000°K.  The  reason  for  no  overshoot  for  case  Q  has  also  been 
discussed  previously. 
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5.4  VALIDITY  OF  CRR 

The  preceding  discussion  of  the  mechanisms  by  which  QSS  is  estab¬ 
lished  showed  that,  except  for  a  limited  range  of  plasma  conditions,  the 
continuum  rate  does  not  deviate  drastically  from  its  QSS  values.  Thus 
one  would  expect  the  results  of  CRR.  theory  to  be  apjiLicable  to  the  free 
electron  density  even  though  the  QSS  is  not  a  physical  reality.  This  in 
turn  implies  that  CRR  can  occur  physically  independent  of  the  QSS  and 
indeed,  the  QSS  may  be  unnecessary  except  as  a  mathematical  conven¬ 
ience  or  in  studying  the  quantum  levels  below  the  critical  level. 

This  is  exactly  the  situation  as  is  shown  in  Figs.  25,  26,  27,  and 
28.  Each  of  these  figures  displays  the  time -dependent  CRR  coefficient, 
a<t),  based  upon  the  electron  density  decay,  a  (t)  =  -rie/ne^,  referred 
to  its  value  at  QSS.  Each  of  the  Figs.  25  through  28  shows  results  for 
the  same  plasmas  shown  in  Figs.  18  through  21.  As  expected,  only 
those  cases  which  showed  significant  continuiim  overshoot  (Table  7,  cases 
C,  D,  I,  J,  and  M)  show  significant  overshoot  of  the  QSS.  Even  so  the 
departure  from  the  QSS  is  less  severe  than  one  would  expect  from  exam¬ 
ination  of  ng/n(l)  shown  in  Figs.  18  through  21. 


Figure  25.  a(t)/a(QSS)  versus  time,  cases  A,  B,  C,  D,  E,  F,  and  G. 


74 


AEDC-TR-76-5 


Figure  26.  a(t)/a(QSS)  versus  time,  cases  H,  I,  J.  and  K. 


The  conclusion  from  this,  that  CRR  can  be  applied  to  the  continuum 
rate  for  low  density  (hq  <  10  1/cm^),  low  excitation  temperature 

(Tg  <  10, 000®K)  plasmas  before  QSS  is  acquired  is  inescapable.  The 
reason  for  this  applicability  of  CRR  to  the  continuum  decay  is  fairly 
straightforward.  If  one  writes  the  EBE  equation  for  the  rate  of  change 
of  the  ground  state  and  factors  out  Ug^  to  put  it  into  the  proper  form  for 
identifying  the  CRR  coefficient,  one  sees  that  the  dominant  contributors 
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Time  (sec) 


Figure  28.  a(t)/a(QSS)  versus  time,  cases  C,  I,  P,  and  Q. 

to  the  ground  state  rate  are  those  terms  involving  the  low-lying  energy 
levels,  particularly  radiation  from  the  2^P  level.  Because  of  the  ener¬ 
getically  close  proximity  of  the  2 level  to  the  2^S,  2^P,  and  2^S  levels 
and  the  significant  times  required  for  them  to  reach  their  QSS  value,  it 
will  require  significant  times  for  the  QSS  to  be  reached,  so  far  as  the 
ground  state  is  concerned.  However,  for  the  "less  collision  dominated" 
plasmas,  the  upper  quantum  levels  are  the  most  important  so  far  as 
the  transfer  of  electrons  between  the  continuum  and  the  bound  states  are 
concerned.  The  lower  levels  populate  radiation  from  the  upper  levels. 
Hence,  the  continuum  develops  equilibrium  with  the  upper  states  very 
rapidly  and  decays  with  them  according  to  CRR  as  they  radiatively  fill 
the  low  levels  until  QSS  for  the  entire  plasma  is  obtained. 

A  dominant  feature  of  the  analysis  of  these  plasmas  has  been  the 
effect  of  the  metastable  levels  upon  the  time  required  to  reach  QSS. 

The  other  levels  adjust  to  some  value  near  or  below  their  QSS  value 
early  in  the  decay,  typically  at  times  of  order  10"®  sec  which  is 
characteristic  of  radiative  decay  times.  Then,  as  the  metastable  levels 
continue  to  populate  radiatively,  coUisional  processes  from  the  metas¬ 
table  levels  to  the  adjacent  levels  and  subsequent  excitation  to  higher 
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levels  serve  to  raise  the  densities  of  these  levels  to  their  final  QSS 
values.  The  amount  of  adjustment  is  dictated  by  the  energy  proximity 
to  the  metastable  levels  and  whether  the  collisional  excitations  were 
governed  by  direct  or  exchange  calculations.  The  importance  of  these 
metastable  levels  upon  the  time  for  a  plasma  to  achieve  QSS  should  not 
be  neglected. 


6.0  PERTUBATIONS 


The  study  of  the  QSS  development  for  the  pure  afterglow  cases  has 
provided  a  great  deal  of  insight  into  the  important  mechanisms  in  a 
decaying  helium  plasma.  Actual  physical  phenomena  will  more  realis¬ 
tically  involve  the  plasma  changing  from  the  QSS  at  one  set  of  plasma 
conditions  to  the  QSS  at  another  set  of  plasma  conditions.  The  study  of 
the  reacquisition  of  the  QSS  after  a  perturbation  of  existing  QSS  condi¬ 
tions  will  enhance  the  insight  gained  in  the  study  of  the  pure  afterglow 
plasmas  for  extension  to  the  physical  situation.  For  this,  the  QSS 
distribution  was  subjected  to  various  perturbations,  and  these  are 
summarized  in  Table  8. 


Table  8:  Summary  of  t^^urbatioris  of  Case  C 


Case 

Initial  Conditions 

CA 

All  excited-state  densities  doubled 

CB 

Ail  excited-state  densities  halved 

CC 

Principal  qbahtiini  level  3  density 
increased  one  order  of  magnitude 

CD 

Principal  quantum  level  3  density 
doubled 

CE 

Principal  quantum  level  6  density 
increased  one  order  of  magnitude 

CF 

3^P  density  increased  one  order  of 
magnitude 

CG 

6^P  density  increased  one  order  of 
magnitude 

CH 

Tg  from  10,0d0°K  to  9,500“K  in  one  step 

Cl 

T  from  lOiOOO^K  to  9,500“K  in  100®K 
increments 

CJ 

Te  from  10,000®K  to  9,9000k  in  10“K 
increments 

CK 

All  K(p;q)  halved 
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6.1  PERTURBATIONS  OF  ENTIRE  DISTRIBUTION,  CASES  CA  AND  CB 

Cases  CA  and  CB,  in  which  the  entire  excited-state  QSS  density 
distribution  from  case  C  was  doubled  and  halved,  respectively,  is  the 
largest  perturbation  since  it  affects  the  entire  distribution.  The  quan¬ 
tities  ne/n(l)  and  a{t)/a{QSS)  for  these  plasmas  are  illustrated  in 
Fig.  29  along  with  the  same  parameters  from  case  C  for  comparison. 
Figures  30  and  31  show  the  time  development  of  quantum  levels  2  and  3 
population  densities  normalized  to  their  Saha  values  for  cases  CA  and 
CB,  respectively.  Included  in  Figures  30  and  31  is  the  time  develop¬ 
ment  of  quantum  level  6,  which  lies  above  the  critical  level. 


Figure  29.  he/n(1)  and  a(t)/a(QSS)  versus  time,  cases  C,  CA,  and  CB. 


The  comparison  shows  that  the  time  to  acquire  QSS  is  quite  similar 
for  both  conditions  in  spite  of  the  fact  that  the  initial  conditions  are 
quite  diverse.  Case  C  started  with  all  excited  states  in  equilibrium  with 
the  continuum,  and  the  approach  to  QSS  is  primarily  determined  by  the 
lower  states.  Cases  CA  and  CB  have  as  their  initial  conditions  the  den¬ 
sities  of  all  excited  states  either  double  or  half  their  QSS  value,  re¬ 
spectively,  and  again  QSS  is  determined  by  the  lower  levels. 
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Time  (sec) 


Figure  30.  n(p)/n(Saha)  versus  time  for  quantum  levels  two,  three,  and  six,  case  CA. 
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For  case  C  the  continuum  rate  and  thus  a  are  determined  by  recom- ' 
bination  into  the  upper  levels  which  then  depopulate  by  radiation  to  the 
lower  levels.  For  case  CA  these  levels  are  overpopulated  by  a  factor  of 
two  and  because  of  ionization  collisions  the  upper  levels  reach  equilibrium 
with  the  free  electrons  early  (t  <  10"®  sec)  in  the  return  to  QSS.  The 
lower  levels  are  overpopulated  and  depopulate  principally  through  coUi- 
sional  excitation  to  a  higher  level  and  then  subsequent  ionization.  Thus, 
in  this  case  the  continuum  rate  is  held  towards  a  more  positive  value 
until  QSS  is  obtained  and  there  is  no  overshoot. 


For  case  CB  the  levels  are  all  underpopulated  by  a  factor  of  two. 
Again  the  upper  levels  establish  equilibrium  with  the  free  electron  den¬ 
sity  earlier  (t  <  10"®  sec)  than  the  time  required  for  the  lower  levels  but 
the  actual  QSS  is  not  obtained  imtil  considerably  later  (t  *  5  x  lO""^  sec). 
The  similarity  in  the  decay  characteristics  of  cases  C  and  CB  for  1  x 
10"®  sec  <  t  <  5  X  10"^  sec  is  worth  noting.  For  case  C  the  energy 
states  above  the  critical  level  were  in  equilibrium  with  the  free  electron 
density  from  t  =  0.0  sec.  For  case  CB  the  initial  conditions  were 
just  half  the  initial  conditions  of  case  C  for  the  states  above  the  critical 
level  and  near  the  same  for  the  quantum  states  below  the  critical  level. 
Consequently,  the  decay  characteristics  of  the  two  cases  are  quite 
different  during  the  early  part  of  the  decay  when  the  upper  states  can 
appreciably  affect  the  characteristics.  However,  by  t  =  1  x  10"®  sec, 
all  levels  with  the  exception  of  the  metastable  levels  have  established  a 
collisional- radiative  balance  among  themselves  and  subsequent  develop¬ 
ments  are  a  result  of  the  coupling  witn  the  metastabie  levels.  Thus, 
because  of  the  similarity  of  the  initial  density  of  the  2®S  level,  the  decay 
characteristics  for  t  >  1  x  10“8  sec  ai’e  quite  similar. 


6.2  PERTURBATIONS  OF  PRINCIPAL  QUANTUM  LEVEL  DENSITY, 

CASES  CC,  CD,  AND  CE 

Cases  CC,  CD,  .and  CE  each  suffered  somewhat  smaller  perturba¬ 
tions,  having  but  one  principal  quantum  level  affected.  Cases  CC  and 
CD  have  the  QSS  density  of  quantum  level  3  (below  the  critical  level) 
increased  a  factor  of  10  and  2,  respectively.  Case  CE  has  the  QSS  den¬ 
sity  of  principal  quantum  level  6  (above  the  critical  level)  increased  by  a 
factor  of  10.  Figure  32  summarizes  the  results  for  these  cases  with 
the  ordinate  the  values  of  the  ratios  n^/n  (1)  and  a  (t)  /  a  (QSS)  as 
before.  The  plots  show  the  not  unexpected  trends  that  case  CE,  in 
which  quantum  level  6  was  increased  an  order  of  magnitude,  comes  to 
QSS  very  rapidly  in  both  derivatives  and  CRR  coefficients.  Cases  CC 
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Figure  32.  rie/h(1)  and  a(t)/a(QSS)  versus  time,  cases  CC,  CD,  and  CE. 


and  CD,  which  had  quantum  level  3  increased  factors  of  10  and  2,  res¬ 
pectively,  take  longer  -  only  slightly  for  case  CD  but  an  order  of  mag¬ 
nitude  for  case  CC.  The  nQ/n(l)  curve  for  case  CD  shows  an  interest¬ 
ing  character  illustrating  the  role  of  bound  states  above  the  critical  level 
as  a  communicator  between  the  lower  level  and  the  continuum.  Early 
in  the  decay  {t  <  10"^®  sec)  few  collisions  have  occurred  which  transfer 
electrons  from  quantum  level  3  to  some  quantum  level  adjacent  to  the 
continuum.  Thus,  there  is  not  sufficient  ionization  from  the  upper  levels 
to  appreciably  affect  the  .continuum  rate.  However,  for  10"^®  <  t  <  3  x 
10"®  sec  there  are  sufficient  numbers  transferred  to  these  upper  levels 
that  ionization  overbalances  recombination  and  the  continuum  rate  act¬ 
ually  goes  positive  for  a  short  period  of  time  before  driving  rapidly  to 
the  QSS  by  t  <•  10”®  sec.  Case  CC  exhibited  this  same  character  but 
because  of  the  higher  initial  overpopulation  of  quantum  level  3,  the 
ratios  could  not  be  included  on  the  scale  of  the  plot.  Case  CE  differed 
in  this  respect  in  that  the  continuiun  rate  showed  a  very  large  positive 
derivative  from  t  =  0. 


Figures  33  and  34  explore  cases  CC  and  CE  a  little  farther,  showing 
the  time  development  of  the  normalized  population  densities  of  the  per- 
timbed  and  two  adjacent  levels.  In  each  the  ordinate  is  the  ratio  of  the 
time -dependent  population  density  to  the  density  in  equilibrium  with  the 
free  electrons.  Figure  33  shows  that  the  singlet  and  triplet  states  have 
slightly  different  decay -characteristics  for  quantum  level  3.  This  is 
primarily  because  of  an  initial  difference  in  population  densities.  The 
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Figure  33.  n(p)/n(Saha}  versus  time  for  33p,  3^  P,  and  4^  P  quantum  levels, 
case  CC. 

transient  coupling  to  other  states  is  illustrated  by  the  development  of 
the  4^P  density  which  shows  an  increase  in  density  peaking  at  t  10"® 
sec  before  falling  back  to  the  QSS  value.  All  higher  levels  in  fact  show 
the  same  general  natiire  as  the  4lP  level,  and  the  QSS  values  are  not 
reached  until  times  of  order  10"^  sec. 

Figure  34  illustrates  the  transient  coupling  between  quantum  level 
6  and  the  adjacent  levels.  Because  quantum  level  7  is  closer  to  6 
energetically  than  quantum  level  5,  it  shows  a  rise  in  population  density 
before  quantum  level  5  and  then,  because  ionization  is  more  efficient 
from  quantum  level  7,  it  has  a  tendency  to  return  to  QSS  before  the 
others.  This  is  also  why  quantum  level  5  tends  to  lag  behind  both  quan¬ 
tum  levels  6  and  7  in  returning  to  QSS. 

6.3  PERTURBATIONS  OF  DENSITY  OF  ONE  EXCITED  STATE,  CASES  CF  AND  CG 

Figure  35  illustrates  the  character  of  the  decay  of  cases  CF  and  CG 
in  which  the  densities  of  the  S^P  and  e^P  states,  respectively,  were 
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Figure  34.  n(p)/n{Saha)  versus  time  for  quantum  levels 
five,  six,  and  seven,  case  CE. 

perturbed  upwards  by  a  factor  of  10.  The  ordinate  again  gives  the  values 
of  the  ratios  ne/n(l)  and  a  (t)/a(QSS). 

These  plasmas  .show  characteristics  similar  to  cases  CC  and  CE  hy 
coming  to  QSS  quite  rapidly  with  the  CRR  showing  the  QSS  value  slightly 
before  the  rate  ratio.  Comparison  with  Fig.  32  shows  that  these  plasmas 
achieve  the  QSS  typically  an  order  of  magnitude  before  cases  CC  and  CE 
which  had  the  entire  quantum  level  elevated  by  a  factor  of  10.  Figures 
36  and  37  illustrate  the  time  development  of  selected  population  densities 
in  the  same  manner  as  Figs.  33  and  34.  Again  the  transient  coupling  is 
evident  between  the  affe,cted  states  and  adjacent  states  with  those  lying 
energetically  below  the  states  lagging  slightly.  Although  exchange 
collisions  are  included  in  the  calculations,  these  rate  coefficients  are 
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Figure  35.  h«/h(1)  and  a(t)/a(QSS)  versus  time, 
cases  CF  and  CG. 


very  small  and  the  effect  of  these  perturbations  upon  singlet  states, 
though  resolvable  in  the  calcvilations,  is  insignificant  compared  to  the 
effect  on  the  triplet  states.  Figure  37  does  not  show  the  effects  upon  a 
higher  level  because,  although  a  slight  increase  is  observable  in  these 
higher  levels,  it  is  fractional  and  would  only  tend  to  clutter  the  figure. 
This  figure  illustrates  the  strong  coupling  between  the  sublevels  of  the 
high  quantum  levels  in  that  the  63  P  level  and  the  other  triplet  states  of 
quantum  level  6  converge  to  the  same  relationship  with  the  continuum  at 
about  3  X  10"  gee  and  then  decay  to  QSS  together. 
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Figure  36.  n(p)/n(Saha)  versus  time  for  S^P,  S^S,  and  43$  quantum 
levels,  case  CF. 
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Figure  37.  n(p)/n(Saha)  versus  time  for  quantum  level  and  triplet  states  of 
quantum  levels  five  and  six,  case  CG. 
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i6.4  ELECTRON  TEMPERATURE  DECAY,  CASES  CH,  Cl,  AND  CJ 

Most  experimental  studies  of  plasmas  are  subject  to  a  decaying 
electron  temperature.  If  the  acquisition  of  the  QSS  under  these 
conditions  requires  times  in  excess  of  the  observed  gradient  times,  then 
application  of  CRR  theory  may  be  inappropriate.  The  effects  of  a 
relaxing  electron  temperature  were  initially  studied  with  the  ERE  with 
two  cases,  CH  and  Cl.  Each  case  used  the  QSS  distribution  from  case 
C  for  the  initial  distribution  but  a  different  lower  electron  temperature. 
Although  the  present  computer  program  will  not  support  a  variable 
electron  temperature,  the  effect  of  an  electron  temperature  change  can 
be  approximated  replacing  the  rate  coefficients  at  the  QSS  tempera¬ 
ture  by  rate  coefficients  calculated  at  some  other  temperature.  Case 
Cl  investigated  the  approach  to  QSS  by  stepping  the  electron  temperature 
from  the  10,  000“K  of  case  C  to  9,  500“K  in  lOO^K  steps.  The  calculation 
was  continued  at  each  temperature  until  ng/h(l)  indicated  the  plasma 
was  within  five  percent  of  QSS.  Case  CH  repeated  the  same  conditions 
except  the  500*’ K  drop  was  made  in  one  step. 

Comparison  of  the  two  cases  is  made  in  Table  9.  The  first  four 
colunms  refer  to  case  Cl  with  the  first  coliunn  indicating  the  electron 
temperature  of  each  succeeding  calciiLation.  The  second  column  is  the 
time  at  which  the  ratio  ng/nil)  came  within  the  five  percent  criterion 
indicated  above.  Since  provisions  were  not  made  to  halt  computations 

Table  9.  Times  for  QSS  to  be  Established  for  Cases  CH  and  Cl 


t( 

*K 

-n^/nCl)  <  1.05) 

sec 

t  (computations 
stopped) 
sec 

Ag/n(l)  t(-ng/n(l)  <  1.05) 

t  stop 

9.9^^ 

4.2“® 

8.1"® 

-1.03 

9.8^ 

8.8"® 

8.8-® 

-1.05 

9.7^ 

1.4-7 

1.8-7 

-1.04 

9.6^ 

1.0-7 

3.3-7 

-1.02 

9.5^ 

6.0“® 

3.0-7 

-1.02 

Sum 

=  3.92-7 

3.1-7 

^Superscripts 

indicate  powers  of  10 

by  which  the  numbers  are  to  be 

multiplied. 
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as  soon  as  the  five  percent  criterion  was  met,  the  third  column  indicates 
the  time  at  which  computations  were  stopped  at  each  temperature  with 
the  fourth  column  showing  the  ratio  ng/n{l)  at  the  time  computations 
were  stopped. 

The  QSS  distribution  at  each  of  these  temperatures  was  used  for 
the  initial  distribution  at  the  next  lower  temperature.  The  last  entry 
in  the  second  column  is  the  sum  of  times  to  reach  QSS  at  the  five  percent 
criterion  and  is  to  be  compared  to  the  single  entry  in  column  five,  which 
is  the  time  required  for  case  CH,  the  single  SOO^K  drop,  to  reach  QSS 
at  the  five  percent  criterion.  Comparison  of  the  second  and  fifth  columns 
of  Table  9  shows  that  approximating  the  Te  change  a  series  of  small 
changes  will  result  in  the  plasma  requiring  longer  to  achieve  QSS  at  the 
final  Tg  than  if  the  change  is  made  in  one  large  jump. 

An  interesting  feature  of  the  decay,  comparing  columns  four  and 
two  of  Table  9,  is  that  the  farther  away  from  QSS  the  initial  conditions 
are,  the  longer  the  time  required  for  the  plasma  to  reachieve  the  QSS. 
Thus,  for  Tg  =  9,  900®K,  it  only  required  4.  2  x  10“^  sec  for  ng/n<l) 
to  come  within  five  percent  of  QSS  (with  ng/n(l)  =  -1.000  at  Te  = 

10,  000“K)  but  it  required  1.4  x  lO""^  for  Tg  =  9,  700®K  (rie/n(l)  =  -1.05 
at  Te  =  9,800®K).  This  suggests  that  the  instantaneous  rates  are  quite 
sensitive  to  the  value  of  the  rate  coefficients  even  though  final  QSS 
rates  may  not  be  appreciably  different.  This  is  supported  by  noting  that 
a  QSS,  the  instantaneous  total  rate  of  change  of  a  quEuitum  level  will  be 
several  orders  of  magnitude  below  the  instantaneous  rate  of  change  of 
either  collisional  population  or  collisional  depopulation.  Hence,  a 
very  slight  change  in  the  rate  coefficient  wiU  cause  appreciable  changes 
in  the  ne/n(l)  characteristic.  Because  of  this,  small  changes  in  Tg 
will  contribute  significantly  to  the  ratio  ng/n(l). 

With  the  electron  temperature  decaying  continuously,  as  in  a  phys¬ 
ical  plasma,  in  the  time  following  a  Tg  change  and  with  the  plasma 
approaching  the  QSS,  there  will  be  another  Tg  change,  albeit  slight, 
with  the  attendant  deleterious  effect  upon  the  rates  and  the  QSS.  These 
observations  lead  to  the  hypothesis  that  in  an  environment  in  which 
the  electron  temperature  decays  continuously,  the  QSS  will  continually 
"chase"  the  Tg  change  and  the  plasma  will  never  satisfactorily  achieve 
QSS  until  the  Tg  time  gradient  becomes  "sufficiently"  flat. 

These  observations  prompted  the  last  temperature  perturbational 
study,  case  CJ.  The  initial  conditions  were  again  the  same  as  for 
cases  CH  and  Cl.  This  time,  however,  the  plasma  was  subjected  to  a 
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constant  temperature  gradient,  -9  x  10^°K/sec,  and  the  temperature, 
increments  were  10®K  drops  from  10,000“K  to  9,  900"K.  The  rate 
coefficients  were  assumed  linear  over  the  1U0°K  temperature  range. 
Note  that  this  gradient  is  only  about  one  third  that  indicated  in  Table  9 
for  the  first  100°  K  drop. 

Figure  38  displays  the  ratio  ne/ri(l)  at  the  point  the  computations 
were  terminated  at  each  temperature.  The  points  at  9,  990°  K  and 
9,  980°K  are  plotted  at  their  respective  plasma  times  of  1. 1  x  10"^  sec 
each  since  those  decays  were  inadvertently  allowed  to  continue  to  a 
plasma  time  of  1  x  10"^  sec.  The  terminal  distribution  at  each  temper¬ 
ature  provided  the  initial  conditions  for  the  next  lower  temperature. 


Figure  38.  ne/n(1)  at  times  computations  terminated,  case  CJ. 


The  figure  shows  that  if  the  QSS  plasma  of  case  C  is  subjected  to 
the  constant  temperature  gradient  of  -9  x  10^°K/sec,  it  will  still 
not  be  to  the  QSS  after  a  one  percent  temperature  drop  at  a  plasma 
time  of  1.  15  X  10"®  sec.  This  time  is  almost  three  times  that  indicated 
in  Table  9  for  the  same  temperature  drop. 

This  case  CJ  reinforces  the  hypothesis  just  presented  that  in  an 
environment  in  which  T©  is  continuously  changing,  the  QSS  will 
"chase"  the  Tg  change.  Although  the  conclusions  must  be  based  upon 
phenomenological  observations,  it  is  apparent  that  in  the  10,000°K  to 
9,  900°K  temperature  range,  the  maximum  temperature  gradient  which 
will  allow  application  of  the  QSS  is  less  than  -9  x  10^®  °K/sec.  Exten¬ 
sion  of  these  observations  to  the  cases  in  which  the  plasma  suffers  a 
500°K  electron  temperature  decay  indicates  that  the  time  for  the  plasma 
to  return  to  the  QSS  would  be  considerably  in  excess  of  the  3.  92  x  10"'^ 
sec  reported  in  Table  9. 
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6.5  PERTURBATION  OF  RATE  COEFFICIENTS,  CASE  CK 

The  discussion  of  the  rate  coefficients  in  Section  3.  0  indicated  that 
they  should  be  overestimated  by  no  more  than  a  factor  of  two.  In 
order  to  investigate  the  effect  of  this  possibility,  the  transient  decay 
of  case  C  was  recomputed  with  all  rate  coefficients,  K(p;q),  reduced 
by  a  factor  of  two.  The  time  development  of  the  population  density  of 
quantum  levels  two  and  three  ratioed  to  the  Saha  equilibrium  value  for 
this  case  CK  is  shown  in  Fig.  39. 

Comparison  of  the  plasma  decay  illustrated  in  Fig.  39  to  case  C 
illustrated  in  Fig.  9  shows  that  the  population  densities  are  spread  out 
more  for  case  CK  than  for  case  C.  For  example,  the  2^8  level  has 
a  QSS  ordinate  of  4.  9  for  case  CK  and  only  2.  7  for  case  C.  This  would 
be  expected  since  lower  rate  coefficients  would  physically  mean  poorer 
coUisional  efficiency,  and  thus  radiation  will  have  a  greater  effect.  An 
interesting  consequence  of  this  is  that  the  2lS  density  has  a  larger  QSS 
value  for  case  CK  than  for  case  C  while  the  2^P  level  has  the  same  QSS 
density  for  both  cases.  Thus  the  CRR  coefficient  will  have  approximate¬ 
ly  the  same  value  for  both  case  CK  and  case  C.  Quantum  level  3  states 


Figure  39.  n(p)/n(Saha)  versus  time  for  quantum  levels  two  and  three,  case  CK. 
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generally  show  lower  densities  for  case  CK  than  for  case  C,  as  one  would 
expect,  with  collisions  being  less  important.  The  other  quantum  levels 
below  the  critical  level  will  also  show  slightly  lower  densities  for  case 
CK  than  for  case  C.  Those  quantum  levels  above  the  critical  level  will 
show  the  same  densities  for  both  cases.  The  critical  level  for  case  CK 
with  the  10  percent  criterion  was  five;  that  for  case  C  was  four. 

Figure  40  shows  the  ng/nd)  and  a(t)/a(QSS)  curves  for  case  CK  and 
case  C.  This  shows  that  case  CK,  with  the  collisional  coefficients  cut 
in  half,  requires  approximately  twice  as  long  to  acquire  QSS  as  does 
case  C.  Fxirther,  the  overshoot  of  case  CK  is  larger  than  that  of  case  C. 

This  is  all  because,  with  lower  collisional  coefficients,  particularly 
the  23s  level  must  achieve  a  higher  density  at  QSS.  This  causes  the 
larger  overshoot  because  the  2^?  density  is  allowed  to  depopulate 
further  before  collisional  processes  balance  radiative  processes. 

A  consequence  of  these  observations  is  that  the  results  for  the 
other  cases  studied  will  generally  be  conservative.  If  the  Gryzinski 
cross  sections  used  in  this  study  are  too  large,  as  anticipated,  the  true 
times  for  the  acquisition  of  the  QSS  will  be  larger. 


Time  (sec ) 


Figure  40.  ha/h{1)  and  a{t}/a(QSS)  versus  time,  case  CK. 
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6.6  ABSORPTION 

Although  absorption  phenomena  require  precise  knowledge  of  geo¬ 
metry  for  explicit  study,  the  effects  can  be  accounted  for  by  simply 
reducing  the  appropriate  Einstein  transition  probability.  Thus,  if  the 
plasma  absorbs  half  of  the  radiation  of  a  specific  wavelength,  this  is 
the  same  as  if  only  half  the  radiation  were  emitted  and  this  is  described 
by  simply  reducing  the  transition  probability  for  the  emission  by  50 
percent. 

The  effects  of  absorption  were  studied  in  this  way  for  a  jdasma  at 
Te  =  10,000‘’K,  no  =  1,5  X  10^®  l/cm^,  one  percent  ionization,  and 
an  equilibrium  distribution  from  the  first  excited  state,  the  same  initial 
conditions  as  case  C.  Four  plasmas  were  studied  in  which  all  radiative 
transition  probabilities  from  an  elevated  bound  state  to  the  ground  state 
were  reduced.  The  plasmas  were  assumed  optically  thin  to  continuum 
radiation,  a  minor  constituent  at  these  conditions.  The  identification 
and  amount  of  reduction  for  each  case  is  given  below: 

1.  Case  C25:  25  percent. 

2.  Case  C50:  50  percent. 

3.  Case  C75:  75  percent. 

4.  Case  ClOO:  100  percent. 

Figure  41  summarizes  the  decay  characteristics  of  these  plasmas. 
The  ordinate  is  the  ratio  ne/n(l).  With  increasing  reduction  in  tran¬ 
sition  probability  the  overshoot  of  QSS  of  both  rig  and  n  (1)  increases 
slightly  through  the  75  percent  reduction.  These  plasmas  all  reach 
QSS  in  very  near  the  same  times, _  again  progressively  longer  with 
increasing  reduction.  The  transition  frpna  the  75  percent  to  the  100 
percent  reduction  shows  a  dramatic  increase  in  the  overshoot  of  the 
ne/n{l)  curve.  The  scale  of  Fig.  41  does  not  permit  showing  the  point 
at  which  ne/n(l)  reaches  its.  ininimuin,  which  is  -30  att  =  1  x  10"8 
sec.  Even  with  such  a,n  increase  in  the  overshoot,  however,  the  plasma 
achieves  QSS  within  five  percent  at  t  =  1.  6  x  10"®  sec,  just  slightly 
greater  than  the  other  conditions, 

The  sudden  change  in  the  ng/n(l)  characteristic  is  easily  explained 
by  observing  the  change  in  the  transition  probabilities.  The  spontan¬ 
eous  transition  probability  for  the  2lp  -  lls  transition  for  example  is 
1.  8  X  10®  1/sec.  In  going  from  optically  thin  to  the  75  percent  reduc¬ 
tion,  the  transition  probability  is  reduced  to  4.  5  x  10®  or  approximately 
one  order  of  magnitude.  In  going  from  the  75  percent  to  the  100  percent 
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reduction,  the  probability  is  reduced  to  zero,  or  eight  orders  of  magni¬ 
tude,  a  much  greater  change.  Hence,  with  absorption  less  that  100 
percent,  there  will  still  be  appreciable  depopulation  of  the  2 level 
and  then  a  subsequent  increase  in  the  2^P  level  as  the  metastable  levels 
increase  in  density.  For  case  ClOO,  however,  the  2 ^P  level  does  not 
have  this  mechanism  to  depopulate  and  hence  there  is  no  ground  state 
rate  overshoot. 

Illustrated  in  Figs.  42,  43,  44,  and  45  is  the  time  development  of 
substates  of  quantum  levels  2  and  3  for  cases  C25,  C50,  C75,  and  ClOO, 
respectively,  in  the  same  manner  that  Fig.  9  illustrates  case  C.  These 
show  the  lessening  2^P  depopulation  with  increasing  absorption  with  the 
somewhat  dramatic  development  of  no  2^P  depopulation  for  case  ClOO 
in  Fig.  45.  The  substates  of  quantum  level  3  show  an  initial  decrease  in 
density  in  Fig.  45  because  of  radiative  transitions  to  substates  of 
quantum  level  2.  Another  interesting  observation  is  the  gradual  increase 
of  the  2^3  QSS  density.  This  is  because  the  2^P  level  provides  a  pro¬ 
gressively  smaller  sink  to  the  veirious  substates  of  quantum  level  2. 
Hence,  coUisional  depopulation  of  the  2^S  level  is  smaller  and  the  state 
must  reach  a  higher  QSS  density  to  balance  radiative  population  with 
collision  depopulation. 
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Figure  42.  n(p)/n(Saha)  versus  time  for  quantum  ievels  two  and  three,  case  C25. 


Figure  43.  n(p)/n<Saha)  versus  time  for  quanturh  levels  two  and  three,  case  C50. 
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Figure  44.  n(p}/n(Saha)  versus  time  for  quantum  levels  two  and  three,  case  C75. 


Time  (sec ) 

Figure  45.  n<p}/n(Saha)  versus  time  for  quantum  levels  two  and  three,  case  C100. 


The  overshoots  of  ng/nd)  for  case  C,  C25,  C50,  and  C75  were 
all  due  to  overshoot  of  both  he  and  n(l).  However,  for  case  ClOO,  the 
overshoot  is  due  entirely  to  overshoot  of  the  continuum  density.  This 
is  because  for  case  ClOO,  the  ground  state  rate  is  determined  princi¬ 
pally  by  super-elastic  collisional  populating  from  all  the  higher  levels. 
There  is  no  way  for  these  levels  below  the  critical  level  to  depopulate 
rapidly  to  a  density  less  than  their  QSS  value  and  then  to  fill  collision- 
ally  with  the  growth  of  the  metastable  level  as  in  case  C.  Instead  the 
ground  state  rate  follows  the  collisional  rate  from  the  upper  levels  as 
they  fill  to  their  QSS  value.  The  large  value  of  the  overshoot  of  case 
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ClOO  is  in  part  due  to  the  value  of  the  ground  state  rate  which  is  an  order 
of  magnitude  less  than  the  cases  which  allow  some  radiative  trauisitions. 
For  those  cases  which  allow  at  least  some  radiation  to  the  ground  state, 
this  radiative  rate  contributes  to  the  ground  state  rate  as  well  as  even¬ 
tually  determining  the  electronic  rate.  For  case  ClOO,  this  radiation  is 
not  allowed,  so  that  the  QSS  rate  is  determined  by  collisions  and  for 
these  conditions  is  approximately  one  order  of  magnitude  slower.  Hence, 
through  the  development  of  the  QSS,  the  electronic  rate  is  principally 
determined  by  radiative  transitions  downward  to  the  metastable  levels. 
This  proceeds  at  a  rapid  pace  until  QSS  is  approached  yielding  an  inordi¬ 
nately  high  continuum  rate.  As  the  QSS  is  approached,  the  low-lying 
levels  become  sufficiently  populated  radiatively  that  they  collisionally 
feed  back  to  the  continuum  and  slowly  reduce  its  rate.  The  low  levels 
continue  to  fill  and  reduce  the  continuum  rate  until  approaches  the  rate 
with  which  the  ground  state  is  filling,  which  is  significantly  smaller  than 
when  Ip  transitions  are  allowed.  Thus,  the  large  overshoot  of  QSS  for 
case  ClOO  is  due  to  the  comparatively  low  QSS  rate. 


The  effect  of  partial  absorption  has  little  effect  upon  a  (t)/a{QSS) 
characteristics  of  the  plasma  showing  the  same  general  character  as 
the  rig/n(l)  curves.  Case  C25  showed  a  rather  significant  overshoot, 
as  would  be  expected  from  the  overshoot  of  the  electronic  rate.  The  QSS 
values  of  the  CRR  coefficient  for  each  of  the  cases  C,  C25,  C50,  C75, 
and  ClOO  is  given  below. 


Case  C: 

a 

Case  C25; 

a 

Case  C50j 

a 

Case  C75: 

a 

Case  ClOO; 

a 

2.  16  X  10"^^  cm^/sec 
1.  90  X  10-J2  cm^/sec 
1.68  X  10", 2  cm^/sec 
1.43  X  10  cm^/sec 
1.01  X  10"  cm^/sec 


The  effect  of  absorprtion  thus  has  noticeable  effects  upon  the  decay 
characteristics  of  an  afterglow  plasma.  A  consequence  of  these  results 
is  that  a  plasma  which  is  100  percent  absorbing  to  these  transitions 
will  have  an  afterglow  lasting  more  than  20  times  as  long  as  one  which 
is  completely  optically  thin. 


7.0  EFFECTS  ON  MEASUREMENTS 

As  has  been  demonstrated  in  the  previous  chapters,  the  QSS  can 
require  physically  significant  times  to  be  established.  This  was 
particularly  illustrated  by  cases  CA,  CB,  CH,  Cl,  and  CJ.^ 


95 


AEDC-TR-76-5 


The  cases  CH,  Cl,  and  CJ  are  perhaps  the  most  meaningful  of  the 
perturbation  studies  made  since  these  examine  changes  in  the  electron 
temperature  which  affects  the  entire  distribution.  Experimental  after¬ 
glow  studies  are  invariably  subjected  to  changes  in  the  electron  temper¬ 
ature.  The  present  studies  have  shown  that  if  the  plasma  conditions 
are  such  that  the  metastable  levels  are  increasing  in  density,  the 
recombination  rates  based  upon  the  electron  density  will  be  higher  than 
recombination  rates  based  upon  the  excited  state  density  distribution. 
Further,  it  was  shown  that,  depending  upon  the  temperature  gradient 
and  the  amount  of  temperattire  change  studied,  the  QSS  might  never  be 
reached  even  from  a  pragmatic  standpoint. 

Although  there  has  been  considerable  investigation  of  helium  recom¬ 
bination,  little  of  that  reported  contains  sufficient  detail  to  examine  the 
results  critically  with  the  ERE.  Also,  much  of  the  reported  work  is  at 
high  enough  pressures  that  the  molecular  helium  ion  is  present  and  must 
be  accounted  for.  An  exception  is  the  work  of  Johnson  and  Hinnov  (Ref. 
22)  in  which  the  excited  state  densities  and  electron  densities  are  studied 
as  functions  of  time  for  low  density  (5  x  10“^  torr  -1  x  10"^  torr)  helimn 
plasma  produced  in  the  C-Stellerator.  Diiring  the  early  time  periods 
of  their  studies  the  temperatm*e  gradient  was  quite  steep,  typically  of 
the  order  of  -1  x  10^“K/sec,  but  flattened  to  only  about  -2  x  104®K/sec 
near  the  end  of  their  studies. 

The  goal  of  their  work  was  to  determine  a  semi-empirical  formula 
for  collisional  cross  sections  for  helium  that  wovild  reproduce  experi¬ 
mental  observations  so  as  to  allow  comparison  of  experimental  awH 
theoretical  predictions.  Their  technique  was  simply  to  use  an  assumed 
cross  section  functional  form  with  certain  adjustable  parameters  for 
the  collisional  transitions  and  then  to  use  the  QSS  assumption  to  calcu¬ 
late  excited  state  densities  and  compare  them  to  measured  densities. 
They  include  exchange  collisional  transitions  between  singlet  and  trip¬ 
let  states  with  the  same  principal  quantum  number  and  orbital  angular 
momentum  for  L  >  2  and  also  have  introduced  cross  sections  for 
optically  forbidden  (AL  ^  1)  transitions  for  "some"  of  their  calcula¬ 
tions.  This  is  in  contrast  to  the  present  work  which  simply  assumed 
Gryzinski's  form- of  the  exchange  calculation  for  all  transitions  which 
are  optically  forbidden.  In  either  case  these  transitions  are  not  partic¬ 
ularly  important  except  for  those  substates  of  quantum  level  2. 

Their  work  was  successful  in  that  they  were  able  to  obtain  fairly 
good  agreement  between  measured  and  calculated  population  densities 
at  low  electron  temperatures  (Tg  <  ~  3,  500 “K).  They  make  little 
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reference  to  the  higher  temperature  data  except  to  note  that  it  was 
necessary  to  include  the  2^S  -  2ls  collisional  transition  to  obtain  agree¬ 
ment  for  the  low-lying  levels.  They  do  not  report  any  population  densi¬ 
ties  below  quantum  level  3. 

Although  the  bulk  of  their  reported  analysis  is  carried  out  on  data 
whose  conditions  lie  below  those  which  can  be  conveniently  examined  by 
the  present  ERE  program,  they  do  report  one  set  of  data  which  can  be 
examined.  The  electron  temperature  versus  time  for  that  set  of  data 
is  shown  in  Fig.  46. 
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In  an  effort  to  examine  this  work  with  the  ERE  solution,  a  QSS 
distribution  was  obtained  at  an  electron  temperature  of  7,  500“K. 
to  compare  to  their  measurements  at  7,  543° K.  They  do  not  report  the 
accuracy  of  their  population  density  measurements  except  through 
reference  to  some  earlier  work  (Ref.  52)  which  reports  an  accuracy  of 
within  20  percent  for  "reasonably  strong  signals"  for  the  absolute 
intensity  measurements.  Electron  density  measurements  were  estimat¬ 
ed  accurate  within  15  percent.  Considering  the  uncertainty  in  their 
electron  density  measurements  and  since  the  states  above  the  critical 
level  can  be  described  by  the  Saha  relation,  which  has  quadratic  electron 
density  dependence,  there  is  an  immediate  uncertainty  of  30  percent  in. 
comparing  their  results  to  the  ERE.  Table  10  shows  the  ratio  of  the 
computed  or  measured  population  density  to  the  corresponding  density 
which  would  be  in  Saha  equilibrium  with  that  state  at  7, 500°K. 

Table  10.  Ratio  of  Excited  State  Population  Densities  to  Saha  Equilibriu.n 
Values  for  ERE  Calculations  and  Measurements,  ne  =  2.5  x 
1/cm3,  Ta  =  7.500°K. 


state 

n  (p)  /n  (8aha) 

ERE 

n  (p)  /n  (Saha) 
measurement^ 

3^P 

O.IB 

0.24 

3^D 

0.18 

0.26 

3^8 

0.50 

0.42 

3^P 

0.40 

0.63 

3^D 

0.40 

0.46 

4ls 

0.72 

1.00 

4^P 

0.74 

0.68 

4^D 

0.74 

0.92 

4^8 

0.83 

1.15 

4^P 

0.87 

1.00 

4^D 

0.87 

1.04 

5^8 

0.92 

1.19 

S^P 

0.93 

1.30 

5^D 

0.93 

0.96 

sh 

0.95 

1.10 

sh 

0.96 

1.10 

5^D 

0.96 

1.05 

6^8 

0.97 

0.87 

6^P 

0.94 

1.05 

6^D 

0.98 

0.74 

6^8 

0.98 

0.96 

6^P 

0.99 

0.92 

6^D 

0.99 

0.87 

9^D 

1.0 

1.22 

9^D 

1.0 

1.22 

^Measurements  from  Johnson  and  Hlnnov  (Ref.  22) 
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Their  data  and  the  ERE  calculations  show  good  agreement  for 
qucuitum  levels  5  and  6  if  allowance  is  made  for  the  uncertainties  in 
their  data  and  one  ignores  the  occasional  flyers  such  as  the  5^P  and 
6^D  levels.  For  these  calculations  the  critical  level  with  the  10  percent 
criterion  would  be  quantum  level  5.  Since  their  electron  density  was 
used  for  both  the  calculations  and  measurements  and  the  ratios  for  the 
high  levels  show  reasonable  agreement,  the  electron  density  measure¬ 
ment  is  satisfactorily  accurate  for  these  comparisons. 

Comparison  of  the  levels  below  the  critical  level  show  that  in 
general  the  measurements  and  calculations  are  still  within  experimen¬ 
tal  error  with  the  exception  of  the  S^P,  4^S,  4^S,  5^S,  and  5^P  states. 

Some  interesting  observations  can  be  made  about  the  distribution. 
The  ERE  calculations  show  the  results  of  very  strong  coupling  between 
the  ^P  and  states  as  well  as  the  ^P  and  states  as  evidenced  by 
their  relationship  to  the  free  electron  density.  This  is  not  necessarily 
evident  in  the  measurements  except  for  the  singlet  states  of  quantum 
level  3  and  the  triplet  states  of  quantum  levels,  4,  5,  and  6.  The 
coupling  is,  however,  evident  within  the  experimental  error  for  all  the 
levels.  The  calculations  show  a  generally  flatter  distribution  than  do 
the  measurements,  which  may  be  indicative  of  an  overestimation  of 
the  cross  sections  by  the  Gryzinski  technique.  Although  Johnson  and 
Hinnov  refer  to  entrapment  of  resonance  radiation  being  present,  they 
account  for  it  in  an  approximate  way  and  its  presence  is  not  well  indi¬ 
cated  by  the  measurements  except  for  the  5^P  and  S^P  levels.  If 
resonance  radiation  entrapment  were  present,  the  ^P  and  measure¬ 
ments  shotild  show  significantly  higher  relationships  to  the  continuum 
than  the  calculations,  which  assumed  no  absorption,  but  the  converse 
is  true. 


The  above  observations  have  been  directed  at  the  absolute  densities 
themselves  which,  because  of  uncertainties  in  both  rate  coefficients 
and  measurements,  are  quite  difficult  to  use  for  reliable  analysis  of 
plasma  conditions.  Of  a  somewhat  more  sensitive  nature  is  an  examina¬ 
tion  of  recombination  rates.  Johnson  and  Hinnov  also  investigated  this, 
using  the  cross  sectional  forms  they  deduced  from  their  measurements 
for  calculations  assuming  QSS  and  comparing  to  direct  measurement 
for  the  electron  density  rate.  A  characteristic  of  their  resialts  is  a 
lower  recombination  rate  from  calculation  than  from  measurement  for 
the  early  periods  of  measurements. 

One  feature  of  early  transitions  already  noted  is  that  when  the 
lower-lying  levels  are  increasing  in  density,  as  with  a  decaying  electron 
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temperature,  the  ground  state  rate  will  underestimate  QSS  conditions 
while  the  electron  density  rate  will  overestimate  the  QSS  condition.  As 
was  shown  earlier,  the  response  of  the  system  to  changes  in  the  electron 
temperature  is  quite  nonlinear  with  respect  to  the  time  scale  of  the 
temperature  gradient.  The  time  to  maintain  QSS  conditions  for  a 
series  of  small  temperature  changes  was  significantly  longer  than  when 
the  temperature  change  was  made  in  one  step. 

To  examine  the  possibility  that  the  QSS  had  not  been  established 
during  the  early  states  of  the  measurements  of  Ref.  22,  three  cal¬ 
culations  were  performed.  The  first  of  them  subjected  the  QSS  dis¬ 
tribution  of  case  M  reported  above  (n^  =  2.  69  x  10^^  1/ cm^)  to  a  sudden 
electron  temperature  drop  of  2,  500°K  to  7,  500°K.  The  second  calcula¬ 
tion  involved  suddenly  dropping  the  electron  temperature  from  7,  500°K 
to  7,  000°K  in  one  step  for  the  QSS  distribution  reported  in  Table  10. 

The  third  calculation  dropped  the  electron  temperature  from  7,  500°K 
to  7,000°K  in  10°K  increments  with  the  plasma  being  allowed  to  come 
to  the  QSS  before  each  succeeding  drop.  The  time  for  QSS  to  be  estab¬ 
lished  at  the  terminal  conditions  for  each  of  the  above  calculations 
was:  8  x  10“^  sec  (AT  =  2,  500°K,  one  step),  7  x  10"^  sec  (AT  =  500°K, 
one  step),  and  1.4  x  10"^  sec  (AT  =  500°K,  50  steps).  From  Fig.  46, 
page  173,  the  experimental  temperature  gradient  for  this  period  was 
about  -2.  1  X  10'  °K/sec.  Thus,  it  required  2.4  x  10"^  sec  for  the 
electron  temperature  to  decay  from  7,  500°K  to  7,  000°K  in  the  experi¬ 
ment. 

The  obvious  non-linearity  of  the  response  of  the  ERE  to  tempera¬ 
ture  changes,  the  close  proximity  of  the  time  calculated  with  the  10°K 
increments  compared  to  the  experimental  times,  and  the  qualitative 
agreement  of  the  predictions  of  this  study  with  their  results  for  the 
early  afterglow  makes  it  inescapable  that  the  QSS  had  not  been  establish¬ 
ed  during  the  early  portions  of  their  work.  Further,  their  observed 
discrepancies  between  the  measured  and  computed  recombination  rates 
are  in  direct  qualitative  agreement  with  the  predictions  of  this  study. 

To  be  more  quantitative  will  require  development  of  an  ERE  computer 
program  with  a  continuously  variable  electron  temperature. 

Since  the  data  (Ref.  22)  used  to  obtain  cross  sections  were  taken  at 
a  later  time  in  the  afterglow  when  the  electron  temperature  was  lower 
and  the  gradient  smaller,  their  work  will  probably  give  adequate  pre¬ 
dictions  for  most  cross  sections.  The  ERE  should  be  used  for  at 
least  the  early  portion  of  their  studies.  This  will  affect  their  reported 
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results  for  the  23s  -  2lS  cross  section  which  were  obtained  in  the 
early  portion  of  their  studies  with  the  QSS  approximation. 


8.0  SUMMARY 

The  system  of  rate  equations  describing  the  transient  behavior  of 
a  singly  ionized,  optically  thin,  monatomic  helium  plasma  has  been 
solved  numerically  for  various  plasma  conditions. 

The  objectives  of  this  study,  as  outlined  in  Section  1.2,  were 
achieved  by  studying  the  early  decay  characteristics  of  plasmas  which 
were  initially  in  a  Boltzmann  distribution  for  all  excited  states  and 
also  by  examining  the  time  to  reacquire  the  QSS  following  an  instan¬ 
taneous  perturbation  of  plasma  parameters.  To  the  author's  know¬ 
ledge,  this  is  the  first  study  of  these  early  plasma  decay  characteris¬ 
tics  or  of  the  perturbational  effects  upon  the  plasma  which  includes  the 
transient  coupling  between  the  various  energy  states. 


8.1  NUMERICAL  VALIDITY 

Because  these  studies  necessarily  involve  very  short  time  periods 
which  are  generally  beyond  experimental  technique,  the  validity  of  the 
solution  was  examined  by  comparing  the  collisional-radiative  recombi¬ 
nation  coefficients  determined  after  the  transient  solution  had  developed 
to  the  QSS  with  the  calculations  of  Chen  (Ref.  21),  who  determined  recom¬ 
bination  coefficients  based  upon  assuming  QSS  conditions.  These  com¬ 
parisons  (Table  3)  showed  favorable  agreement.  Stability  of  the  calcu¬ 
lations,  an  important  consideration  in  investigations  of  this  type,  were 
examined  under  two  criteria  on  the  time  increment  allowed  for  the 
calculation.  The  solution  was  found  to  adequately  reproduce  itself  under 
the  two  criteria.  Hence,  it  is  concluded  that  the  numerical  results 
presented  herein  do  represent  physical  processes  to  the  limit  of  the 
accuracy  of  the  parameters  used  in  the  calculation. 


8.2  PURE  AFTERGLOWS 

The  study  of  the  pure  afterglows,  in  which  the  plasma  relaxes  from 
an  initial  LTE  distribution,  given  in  detail  in  Section  5.0,  provided  a 
response  to  the  specific  objectives  outlined  in  Section  1.2.  These  may 
be  summarized  as  follows: 
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1.  There  is  no  single  common  mechanism  by  which  the 
QSS  is  established.  However,  the  establishment  of 
the  QSS  is  principally  controlled  by  the  time -depen¬ 
dent  characteristic  of  the  22s  density  and  its  effects 

on  the  densities  of  the  several  adjacent  quantum  levels, 

2.  Generally,  the  results  of  CRR  theory  can  be  applied 
to  helium  plasmas  before  the  QSS  condition  for  ex¬ 
cited  atomic  state  densities  is  in  fact  reached.  This 
is  becaused  the  continuum  rate  is  largely  determined 
by  recombination  into  the  high-lying  quantum  levels 
and  the  QSS  configuration  is  established  rapidly  for 
these.  The  ground  state  rate,  however,  is  largely 
determined  by  the  low-lying  quantum  levels  and  the 
QSS  distribution  for  these  requires  a  comparatively 
appreciable  time  to  be  established. 

3.  The  time  required  for  the  QSS  to  be  established  varied 
with  the  plasma  conditions,  being  shortest  for  those 
plasmas  which  were  strongly  collision  dominated  and 
surprisingly  long  for  those  which  were  least  collision 
dominated.  A  summary  of  these  times  for  each  plasma 
to  achieve  the  QSS  as  well  as  the  time  at  which  CRR 
could  be  applied  to  the  electron  density  decay  is  listed 
in  Table  11. 

4.  The  processes  by  which  plasmas  decay  are  indeed 
complex  and  subtle.  The  study  of  the  ERE  solutions 
has  illustrated  this  because  of  the  ability  to  examine 
each  individual  process  singly  as  well  as  its  coupling 
effects  with  other  mechanisms  in  the  plasma  decay. 
Very  rarely  in  low  density  helium  plasmas  such  as 
were  studied  here  can  one  mechanism  (e.  g.  ,  2^S  -  2^P 
collisional  transfers)  be  singled  out  as  a  dominant 
feature  of  the  decay  but  rather,  the  total  coupled  fea¬ 
tures  of  the  decay  must  be  considered  in  the  transient 
mode.  Steady-state  theories  largely  mask  the  inter¬ 
play  between  the  various  states  and  cannot  account  for 
such  phenomena  as,  say,  the  momentary  radiative 
overdepopulation  of  the  2lp  density  and  then  the  sub¬ 
sequent  rise  in  the  2^P  density  via  23s  -  2^P,  2^S  - 
23p  _  2^P,  2^S  -  2^S  -  2^P,  etc.  collisional  transfers. 
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Table  11.  Summary  of  Plasma  Conditions,  Time  to  Acquire  Quasi-Steady 
State,  and  Time  at  Which  Collisionai-Radiative  Recombination 
Theory  Applies  to  Free  Electron  Density 


Case 

T 

e 

n 

o 

n 

e 

Time  (QSS) 

Time  (CRR) 

°K 

l/cm^ 

l/cm^ 

sec 

sec 

A 

1.0^ 

1.63^’ 

1.63^^ 

2.0”® 

1.8-8 

B 

1.0^ 

4.40^® 

4.40^** 

7 . 0”® 

2 ,0”® 

C 

1.0^ 

1.50^® 

1.50^“’ 

8 .0"’ 

3.0"’ 

D 

1.0^ 

6.59^5 

6.59^3 

I.3-® 

1.2-6 

E 

i.c-* 

14 

5.00-^^ 

5.00^^ 

3.1”^ 

5.0-6 

F 

1.0^ 

1.50^^ 

I.50I2 

1.2"^ 

6.5-8 

G 

1.0^ 

5.00^3 

5.00^^ 

3.5"^ 

2.2-’ 

H 

1.4^ 

I.63I’ 

1.63^^ 

2.7"® 

1 . 

I 

1.4^* 

1.50^® 

1.50^'' 

2.0“^ 

2.0-’ 

J 

1.4^ 

S.OO^'’ 

5.00^^ 

1.8~^ 

1.0-6 

K 

1.4^ 

5.00^3 

5.00^1 

2.1-' 

1.0-' 

L 

1.0^ 

2.69^^ 

5.38^^ 

3.0”® 

2.5-8 

M 

i.o'^ 

2.69^^ 

2.69^3 

4.3-^ 

3.0-6 

N 

1.0^ 

2.69^^ 

2.69^2 

6.0"^ 

6.0-6 

0 

1.0^ 

2.69^^ 

2.69^^ 

6.6"' 

2.5-’ 

P 

8.0^ 

1.50^® 

1.50^' 

6.6"^ 

3.0"^ 

Q 

6.0^ 

1.50^^ 

1.50^^ 

1.2-6 

3.0-8 

In  addition  to  the  listing  above,  the  pure  afterglow  studies  showed 
an  additional  unexpected  phenomena.  This  was  overshoot,  in  which  the 
plasma  decay  would  overshoot  QSS  conditions  and  then  return  to  the  QSS 
at  a  later  time.  It  was  found  that  this  phenomena  was  physical  rather 
than  artificial.  The  study  of  overshoot  can  be  summarized  thus: 

1.  Plasmas  which  are  strongly  collision  dominated  at 
their  initial  conditions  will  not  show  overshoot  of  the 
QSS. 

2.  Plasmas  whose  initial  conditions  are  such  that  there  is 
no  clearcut  domination  of  radiation  or  collisions  show 
the  strongest  continuum  rate  overshoot. 

3.  As  excitation  temperature  increases  in  the  range  8,000°K 
<  Tg  <  14,  000°K,  overshoot  of  the  continuum  rate  becomes 
progressively  greater  because  of  improved  collisional 
efficiency. 
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4.  For  lower  electron  temperatures  Tg  <  8,000°K,  overshoot 
becomes  progressively  less,  and  there  is  no  overshoot 

of  the  continuum  rate  at  Tg  =  6,000°K  because  of  poor 
collisional  efficiency. 

5.  The  ground  state  rate  shows  generally  greater  over¬ 
shoot  as  the  collision  dominance  decreases  although 
there  are  inconsistencies  in  this  behavior  for  increas¬ 
ing  electron  temperature. 

The  plasma  conditions  and  the  magnitude  of  overshoot  encountered 
in  this  study  were  summarized  in  Table  7. 


8.3  PERTURBATIONS 

The  perturbational  studies.  Section  6.0,  in  which  the  plasma  relax¬ 
es  after  a  perturbation  of  some  plasma  parameter,  reinforced  the  con¬ 
clusions  drawn  from  the  pure  afterglow  studies.  They  additionally 
indicated  that  an  environment  supporting  a  transient  electron  tempera¬ 
ture  could  have  a  significant  effect  upon  the  QSS.  The  findings  of  the 
perturbational  studies  may  be  summarized  as  follows: 

1.  Parametric  changes  affecting  the  entire  distribution 
require  significantly  longer  times  to  achieve  the  QSS 
than  if  the  perturbation  affects  the  density  of  just  one 
quantum  level . 

2.  Perturbations  of  the  density  of  a  single  quantum  level 
will  induce  a  ripple  effect  through  the  rest  of  the  distri¬ 
bution  because  of  the  transient  coupling  to  the  density 
of  the  selected  level. 

3.  Errors  in  the  rate  coefficients  will  have  noticeable 
effects  upon  the  transient  decay  of  the  plasma.  There 
can  be  an  appreciable  effect  upon  the  QSS  distribution 
of  densities  for  those  levels  below  the  critical  level 
but  no  discernible  effect  upon  the  densities  of  levels 
above  the  critical  level.  The  location  of  the  critical 
level  can  be  affected.  It  was  estimated  that  the  rate 
coefficients  used  in  this  study  were  too  large;  hence, 
the  times  calcialated  in  this  study  will  be  short  when 
compared  to  the  physical  situation. 
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4.  The  mathematical  and  thus  the  physical  system  show 

a  strong  non-linear  response  to  changes  in  the  electron 
temperature.  The  effects  of  a  continuously  changing 
electron  temperature  were  approximated  by  a  series 
of  step  function  changes.  For  a  series  of  small  changes 
of  electron  temperature  the  plasma  required  appreci¬ 
ably  longer  times  to  reacquire  QSS  than  for  a  single 
step  with  the  same  total  change.  The  studies  led  to 
the  observation  that  in  an  environment  with  a  continuous¬ 
ly  changing  temperature  the  QSS  may  never  be  reached. 

5.  Absorption  of  resonance  radiation  has  a  calculable 
but  not  significant  effect  on  the  times  to  achieve  the 
QSS.  Absorption  of  resonance  radiation  can  have 
effects  on  the  QSS  density  distribution,  causing  the 

and  the  adjacent  states  to  maintain  higher  popvila- 
tion  densities. 


8.4  EFFECTS  ON  IVIEASU  REM  ENTS 

The  studies  of  the  various  plasma  conditions  indicated  that  when 
the  conditions  were  such  that  the  metastable  states  are  increasing  in 
density,  the  recombination  rates  based  upon  the  electron  density  wUl 
be  significantly  higher  than  recombination  rates  based  upon  population 
densities.  This  predicted  effect  was  observed  in  data  published  by 
Johnson  and  Hinnov  (Ref.  22).  In  calculations  with  the  ERE  at  their 
conditions  a  QSS  distribution  was  established  which  agreed  with  their 
measurements,  within  experimental  error.  Their  work  was  subjected 
to  severe  temperature  gradients  and  because  of  the  qualitative  agree¬ 
ment’  of  their  results  with  the  findings  of  this  work,  it  is  believed 
that  they  did  not  have  QSS  conditions,  at  least  in  the  early  stages  of  the 
decay. 

The  findings  of  this  study  have  shown  that  the  concept  of  the  plasma 
passing  from  one  QSS  to  another  QSS  as  the  plasma  decays  is  not  as 
"almost  instantaneous"  as  is  usually  assumed.  However,  the  results 
of  CRR  theory  can  generally  be  applied  to  electron  recombination  rates 
without  the  QSS  being  physically  necessary.  In  studies  in  which  plasma 
parameters  are  changing  rapidly,  great  care  must  be  exercised  in 
relating  density  measurements  of  quantum  states  below  the  critical 
level  to  the  electronic  recombination  rate  via  the  QSS  assumption.  This 


105 


AEDC-TR-76-5 


■will  be  particiiLarly  true  in  apparata  in  which  severe  electron  temperature 
gradients  are  present  or  when  there  are  short  time  scales  present. 


8.5  CONCLUDING  REMARKS 

The  study  has  provided  detailed  information  and  insight  into  plasma 
decays  and  the  mechanisms  by  which  the  nonequUibrium  eigenstate 
population  density  distribution  function  is  established.  These  mechanisms 
have  important  fundamental  implications  for  future  work  in  plasma 
decays  and  their  effects  cannot  be  fully  assessed  without  a  transient  solu¬ 
tion  to  the  ERE. 

Several  implications  for  future  applications  beyond  the  original 
intent  of  the  work  may  also  be  drawn.  Under  proper  experimental  con¬ 
ditions  and  with  an  analytic  approach  similiar  to  that  used  here,  a  very 
critical  detailed  examination  of  rate  coefficients  may  be  made.  This 
work  would  not  have  to  necessarily  be  restricted  to  plasmas  but  the 
approach  is  applicable  to  other  chemical  processes.  The  ripple  effect, 
discovered  in  the  calculations  concerned  with  perturbations  of  a  specific 
energy  level,  suggest  mechanisms  by  which  specific  rate  coefficients 
may  be  determined  through  a  selective  absorption  experiment.  Another 
application  which  suggests  itself  is  the  study  of  the  processes  involved 
in  popvilation  inversions  with  a  view  toward  prediction  of  possible  last¬ 
ing  transitions. 

Of  a  more  immediate  nature,  the  approach  of  this  study  should  be 
put  to  a  critical  test  for  comparison  of  calculations  and  experiment. 
Because  physical  plasmas  are  generally  characterized  by  a  spatial 
density  variation  and  invariably  by  electron  temperature  variations, 
the  work  reported  here  must  be  considered  as  a  first  approximation  to 
actual  physical  processes,  and  correct  modeling  will  need  to  account 
for  the  effects  of  these  additional  transient  phenomena.  The  need  to 
consider  continuous  variations  in  density  and  electron  temperature  as¬ 
sumes  added  importance  in  light  of  the  observations  of  the  effects  of 
step  function  electron  temperature  changes  and  the  physically  significant 
times  for  the  QSS  to  be  established.  The  development  of  these  capabili¬ 
ties  and  completion  of  the  critical  evaluation  will  provide  a  basis  for 
extension  to  more  complex  species,  such  as  diatomic  molecules,  so 
that  the  internal  energy  distribution  function  development  in  the  non- 
equRibrium  environment  may  be  studied. 
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APPENDIX  A 

NUMERICAL  TECHNIQUES 
1.  RATE  COEFFICIENTS 

The  rate  coefficient  integrations  indicated  in  Section  3,  0  were  per¬ 
formed  using  32 -point  Laguerre-Gauss  quadrature.  Quadrature  tech¬ 
niques  are  described  in  most  numerical  analysis  books  (Ref.  53,  for 
example)  and  derivation  will  not  be  given  here.  The  basis  of  the  tech¬ 
nique  is  to  use  a  weighting  function  so  that 


/  w(x)f(x)dx  =  E  w(x.)f(x.)  (A-1) 

^  J  i=l  ^  ^ 

a 

where  w(x)  is  the  weighting  function  and  w(xi)  is  the  vsilue  of  the  weight¬ 
ing  coefficient  at  the  point  x^^. 

It  turns  out  that  if  a  weighting  function 

w(x)  =  e  ^  (A-2) 

is  used,  the  weighting  coefficients  are  a  function  of  the  Laguerre  polyno¬ 
mials  (Ref.  53) 

(n!)^  X. 

w(x.  )  =  - ^  (A-3) 

where  the  Xj^  are  the  roots  of  the  n^^  Laguere  polynomial. 

Further,  use  of  Laguerre-Gauss  quadrature  is  valid  on  the  interval 
(0,  oo).  The  integrations  in  Section  3.  0  do  not  extend  to  oo,  but  the  lower 
limit  is  one.  Compatability  with  the  interval  (0,  °o)  is  accomplished  by 
a  simple  transformation  of  coordinate. 

Various  quadrature  techniques  are  generally  part  of  the  library  in 
numerical  analysis  laboratories.  The  particular  one  used  here  was 
generated  by  K.  R,  Kneile,  Analyst,  Central  Data  Processing,  Arnold 
Engineering  Development  Center  (AEDC),  and  is  in  subroutine  form  for 
ease  of  use.  The  32  abcissas  and  the  values  of  the  weighting  coefficients 
are  given  in  Table  A-1. 
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Table  A-1.  Abscissa  and  Weighting  Coefficients  for  32-Point  Laguerre-Gauss  Quadrature 


i 

X  . 

1 

w(x^) 

i 

X  . 

1 

w(x^ ) 

1 

0.04448936583326702 

0.1141871057681048 

17 

22.63088901319677 

2.884392092922042 

2 

0.2345261095196185 

0.2660652168976152 

18 

25.62863602245925 

3.113261327039586 

3 

0.5768846293018864 

0.4187931373248530 

19 

28.86210181632347 

3.356217692595803 

4 

1.072448753817818 

0.5725328464998047 

20 

32.34662915396474 

3.615869856484269 

5 

1.722408776444645 

0.7276487883809713 

21 

36.10049480575197 

3.895513044948550 

6 

2.528336706425795 

0.8845367193402497 

22 

40.14571977153944 

4.199394104711585 

7 

3.492213273021994 

1.043618875892077 

23 

44.50920799575494 

4.533114978534362 

8 

4.616456769749767 

1.205349274152353 

24 

49.22439498730864 

4.904270287611245 

9 

5.903958504174244 

1.370221338521781 

25 

54.33372133339691 

5.323500972023666 

10 

7.358126733186241 

1.538777256468645 

26 

59.89250916213402 

5.806333214233621 

11 

8.982940924212596 

1.711619352686457 

27 

65.97537728793505 

6.376614674159653 

12 

10.78301863253997 

1.889424063449484 

28 

72.6876280906627 

7.073526580707242 

13 

12.76369798674273 

2.072959340246534 

29 

80.18744697791352 

7.967693509295901 

14 

14.93113975552256 

2.263106633996964 

30 

88.73534041789240 

9.205040331278190 

15 

17.29245433671531 

2.460889072488236 

31 

98.82954286828397 

11.16301309076787 

16 

19.85586094033605 

2.667508126397117 

32 

111.7513980979377 

15.39018041526064 

2.  TRANSIENT  ERE 

The  numerical  technique  used  to  effect  the  solution  to  the  ERE  was 
the  modified  Euler's  method.  This  technique  has  enjoyed  some  success 
in  solution  of  reacting  flow  problems  at  AEDC  and,  because  of  the  simi¬ 
larity  in  the  mathematical  systems,  the  use  of  the  technique  was  sugges¬ 
ted. 


To  describe  the  application  of  the  method  to  the  problem  at  hand, 
it  is  convenient  to  attack  the  problem  in  more  general  form.  Thus, 
the  various  quantities  will  be  described  in  functional  form  and  as 
independent  or  dependent  variables  rather  than  closely  tied  to  physical 
interpretation.  The' system  of  Eos.  (7)  and  (8)  can  thus  be  written 

aj  =  f (y)  (A-4) 

where  the  single  bar  denotes  a  vectoral  quantity,  y  is  the  dependent 
variable,  x  is  the  independent  variable,  and  f(y)  is  the  functional  form  of 
the  ERE. 

The  modified  Euler's  method  proceeds  straightforwardly  from  dif¬ 
ference  formxila  considerations  and  is  written 

?i+i  =  ^  ^  . 
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whe^re  the  subscript  i  denotes  some  point,  Xj[,  at  which  the  solution  is 
known,  i  +  1  the  next  point,  at  which  the  solution  is  desired,  and 

Ax  is  the  difference  in  the  points,  -  Xj^.  This  is  obviously  an  im¬ 

plicit  method  since  a  function  of  the  ui^nown  solution,  y^+i,  resides  on 
the  right-hand  side  of  Eq.  (A-5). 

The  in  Eq.  (A-5)  can  be  approximated  by  expanding  f(y)  in 

Taylor's  Series  in  terms  of  the  solution  and  keeping  the  first  two  terms, 

f  (y)  =  f  (Yq)  +  f '  (Yq)  (y  -  Yq)  (A-6) 

where  f  (jq)  is  the  Jacobian 

\  ^  /y  =  Y 

Jf  Jf  o 

Thus,  f(y^-|-i)  can  be  written 


Defining  the  constant  matrices 

A  =  f(y.)  -  f  (y^)  7. 


and 


B  =  f ' (y^) 


f(yi+l)  is  written 


f(Yi+l)  =  A  + 


B  y 


i+1 


(A- 9) 


(A- 10) 


(A- 11) 


Note  that  with  the  definitions  of  A  and  B,  f(yj^)  can  be  written  in  the  same 
manner : 


f(y^)  =  A  +  B  Y^ 


(A-12) 
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Substituting  Eqs.  (A-12)  and  (A-11)  into  Eq,  (A-5)  and  solving  for 
one  has 


Ax 


-1  r- 


Ax 


=  (I  -  B)  "  [y.  +  (2A  +  B  y.)] 


(A-13) 


where  I  is  the  unit  matrix.  Thus,  the  unknown  solution,  yi+i,  is  known 
as  a  function  of  the  solution  of  the  previous  point,  y^,  and  the  indepen¬ 
dent  variable  increment.  Ax.  The  B  matrix  is  readily  evaluated  analy¬ 
tically  if  the  function  f(y)  is  known,  as  it  is  in  this  case. 


Relatively  simple  considerations  of  the  approximations  involved  in 
establishing  Eq.  (A-13)  yield  a  relationship  between  the  approximate 
error  and  Ax  and  has  been  worked  out  by  Mr.  F.  C.  Loper,  of  AEDC. 
The  error  associated  with  expressing  yi+i  in  the  manner  of  Eq.  (A-5) 
is  of  order  Ax^.  The  error  associated  with  expressing  f(^  in  the  man¬ 
ner  of  Eq.  (A-8)  is  of  order  (y^+j  -  Yi)^  which  is  of  order  Ax^.  This 
error  in  f(y^+]^)  substituted  into  Eq.  (A-5)  shows  the  subsequent  error 
in  y^+i  because  of  the  local  linearization  of  f(^+i)  is  thus  also  of  order 
Ax3,  Hence,  the  error  in  y^^.^  is  order  Ax^^and  the  error  can  be 
written 


^i+1 


^i+l 


3 

i+1 


(A- 14) 


where  g|-|-2  is  some  unknown  function.  Relating  this  to  the  error  in 
y^  and  solving  for  Ax^+i  yields 


1  V3 


(A-15) 


Assuming  the  g  function  does  not  vary  rapidly,  g^+j 

1/3 


gi;  then 


Ax. 


^i  +  l 


E. 

1 


(A-16) 


So,  with  a  desired  accuracy  in  y^+i  and  a  known  accuracy  in  y^,  which 
can  be  evaluated  from 


^  I  "  f^i-i  +  ^  (f(yi-i)  +  f(yi))]  (A-17) 
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where  f(yj^)  is  evaluated  exactly,  a  Ax  yielding  the  approximate 
accuracy  can  be  obtained.  It  should  be  noted  that  Ej  in  Eq.  (A- 17)  is  in 
fact  a  vector.  In  practice,  one  uses  the  largest  element  of  Ej^  from  Eq. 
(A- 17)  for  the  estimation  of  Axj^+j^  in  Eq.  (A-16). 

The  identification  of  the  elements  of  the  y  with  the  population  densi¬ 
ties  n(p),  X  with  time,  and  f(y)  with  the  right-hand  side  of  Eqs.  (7)  and  (8) 
in  Section  2.  0  complete  the  description  of  the  technique  in  terms  of  the 
physical  parameters. 


117 


